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| Introduction 


True PAPER DESCRIBES the design and test work con- 
ducted during a program to develop a three-stage, liquid-cooled 
gas turbine. It is a continuation of a paper presented to the 
ASME three years ago [1].!_ At that time a single-stage turbine 
had been operated successfully over 200 hours. 

The program, sponsored by the Bureau of Ships, U. 8. Navy, 
had as its objectives the development of a gas turbine that 
would utilize little or none of the so-called strategic materials: 
Nickel, chromium, molybdenum, tungsten, cobalt, columbium, 
beryllium. The rotor was to be cooled with a closed cooling 
system. Methods of stator cooling were optional. The gas tur- 
bine was to utilize the compressor developed by Solar Aircraft 
Company for its 500-hp Jupiter? gas turbine and over-all per- 
formance was not to be sacrificed to make cooling possible. 

Considerable development work had been done both in this 
country and abroad on the subject of liquid cooling. These 
efforts had been aimed primarily at solution of the theoretical 
problems associated with cooling. Data on internal and external 
heat-transfer coefficients had been obtained and detail studies of 
the mechanism of cooling had been conducted at the Lewis Labo- 
ratory of the NACA (now NASA). Almost all of this earlier 
work had been conducted with research equipment, and prac- 
tical designs had not been evolved. An excellent summary of gas- 
turbine cooling research was presented to the ASME [2]. 

The program conducted at Solar was intended to produce a 
practical turbine based on the basic data available from NACA 
and others. To this end a design was selected with considera- 
tions for data taking made secondary to economical manufacture. 
Initial design considerations ruled out single-piece turbine con- 
struction with blade, rotor, and shaft, all produced from a single 
billet. Instead, a more conventional design and method of 
fabrication were sought that would permit blades and rotor disks 
to be readily changed. 


It General Cycie and Cooling Considerations 

A Design Point. Though the use of an existing compressor was 
specified, thus limiting the freedom of design point selection, tur- 
bine inlet temperature was not specified. From combined cycle 


1 Numbers in brackets designate References at end of paper. 

2 Registered trade-mark. 

Contributed by the Gas Turbine Power Division and presented 
at the Gas Turbine Power Conference and Exhibit, Cincinnati, Ohio, 
March 8-11, 1959, of Tae American Society or MECHANICAL 
ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, Decem- 
ber 31, 1958. Paper No. 59—GTP-1. 
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Development of a Three-Stage 
Liquid-Cooled Gas Turbine 


A three-stage gas turbine was designed to operate at 1750-F inlet temperature with 
nonstrategic, low alloy components. 
system, and the stators are cooled with air. A practical design for the rotor was evolved 
and tested. Heat-transfer data correlated well with expected values. 


The rotor blades are cooled with water in a closed 


and cooling considerations, the turbine inlet temperature was set 
at 1750 F. Cycle design data are given in Table 1. 


Table 1 Turbine design data 


Ambient pressure 
Ambient temperature 
Air-flow rate 
Compressor inlet pressure 
drop 
Compressor pressure ratio 
Combustor pressure drop 
Combustion efficiency 
Turbine inlet temperature 
Turbine inlet pressure 
Turbine gas flow rate 
Turbine back pressure 


Se 


per cent of inlet pressure 


B Cooling Methods. [Extensive studies were conducted to es- 
tablish the cooling systems for the rotors and stators. For the 
rotors, a closed system was specified. To be established were the 
cooling fluid and type of system. Water was selected and a 
forced-convection system established. 

The advantages of water were obvious: Availability, economy, 
good heat-transfer medium, and no special handling precautions 
required. Its major disadvantages were its low boiling point and 
its exceedingly good heat-transfer characteristics. To eliminate 
boiling within the cooling passages where pressure was not availa- 
ble, fluid temperatures had to be kept below 200 F. At this 
temperature, the low turbine blade temperatures resulted in an 
excessive amount of heat being extracted by the coolant. This 
negative reheat effect lowered over-all performance by about 20 
per cent. 

A forced-convection cooling system was selected over the 
free-convection or mixed systems because the mechanism of forced 
convection was better understood and the blind holes used in the 
free-convection system were considered too vulnerable to plug- 
ging. 

The material limitations required the use of low alloy steels. 
Hence temperatures in highly stressed regions had to be main- 
tained below 800 F, «nd peak metal temperatures limited to about 
1100 F. The materials selected are discussed in reference [1]. 

The calculated variation of coolant temperature rise with flow 
for each of the three stages is given in Fig. 1. Also shown are the 
prevalent or mean blade temperatures. For the first-stage rotor 
blades, trailing edge tip temperatures of 1100 F were estimated. 

Air cooling was selected for the stators, as it appeared that heat- 
transfer losses could be minimized. Water, having too good a 
conductivity, would over-cool the stators. Inasmuch as stress 
levels are lower in the stators, higher temperatures could be 
tolerated. Air cooling was attractive if the coolant could be used 
at a pressure lower than that of the gas stream; with high pres- 
sure air the penalty involved in compressing the air was too great. 
Accordingly, the stators were designed to use air at a pressure 
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psia 
deg F 
Ib/eec 
in. H,O 
per cent 
deg F 
sia 
fb /sec 
in. 


PREVALENT BLADE TEMPERATURE dog F 
COOLANT TEMPERATURE RE dog 
| | 
4 


COOLANT FLOW RATE ib/sec COOLANT FLOW RATE COOLANT FLOW RATE sec 
&. STAGE! STAGED ©. STAGE 


Fig. 1 Variation of 


comparable to compressor second-stage bleed, The use of this 
low-pressure air required sealing to prevent leakage of the high- 
pressure hot gases into the cooling passages. 

Using NACA data, the stator blades were designed with in- 
ternal extended surfaces. 

C Heat-Transfer Rates and Over-All Performance. Heat-transfer 
calculations were used to establish coolant flow rates and heat re- 
jection rates. The cooling losses in the turbine were calculated 
as shown in Table 2. 


Table 2 Turbine heat rejection rates (Btu/sec) 


I I 

Stator shrouds................. 28.7 14.4 9.0 
147.8 98.4 80.2 


In computing turbine performance the following assumptions 
were used: 

1 The total temperatures entering the rotor would be reduced 
by the heat transfer in the preceding stator. 

2 The effective total temperature in the rotor would be further 
reduced by an amount equivalent to one half of the rotor-blade 
heat-transfer rate. 

3 The total temperature entering a stage would be reduced 
by the total of all preceding transferred heat. 


Over-all power plant performance is given in Table 3. 

All the design calculations were conducted on the basis of the 
gas turbine used in a shaft power output configuration. During 
the course of the program it was found expedient to use a larger 
compressor and to extract power as bleed air from the compressor. 
This eliminated the need for a new gear box, and a Solar T-300J 
air-bleed compressor and accessory drive gear box were used. In 
this configuration, the use of air bleed lowered the turbine pres- 
sure level and thus the heat-transfer rate. The reduction was 
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lant temperature rise and prevalent blade temperature with coolant flow rate 


Table 3 Design power-plant performance 


Turbine power output hp 1412 
Compressor power hp 864 
Cooling air power hp 17 
Power consumed in bearings hp 15 
Liquid coolant pumping hp 5 
Oil pumps hp 1 

902 
Shaft power available hp 510 
Gear loss hp 10 
Net power hp 500 
Fuel rate lb/hr 479 
Specific fuel consumption Ib/hp-hr 0.96 
Over-all thermal efficiency _ per cent 14.4 


calculated to be 10 per cent in the first stage and less in the 
second and third. 


A RotorBlades. For high efficiency a three-stage reaction 
turbine was designed. For effective cooling of the rotor blade 
trailing edges, fairly thick trailing edges were specified. The 
high reaction and the low velocities in the turbine allowed these 
relatively thick blades to be used without excessive losses. Other 
than this compromise, fairly conventional design practice was 
used, based on free vortex flow and resulting in twisted rotor 
blades. . 

The rotor blades were fabricated from precision castings. 
Steps in the fabrication of rotor blades include: 


1 Precision cast from low alloy steel. 

Anneal the casting. 

Drill coolant holes and manifold holes, 
Chromize for corrosion protection. 

Braze tip cap and manifold end plugs. 

Heat treat to tensile strength of 130,000 psi. 
Crush grind rotor-blade base. 

Braze O-ring studs into bottom of base. 


2 
3 
4 
5 
6 
7 
8 
9 Pressure and flow test. 
A 


sketch of the cooling passages in a roter blade is shown in 
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Fig. 3 Section through first-stage rotor blade 


Fig. 2. Six radial holes were used for the first-stage blades and 
four for the second and third stages. A photograph of a section 
through a blade is shown in Fig. 3. 

During early tests, several blades were fabricated with a copper 
insert in the trailing edge. These were designed to provide a 
high conductivity path so as to reduce trailing edge temperatures 
as suggested in reference [3]. Slots 0.060 in. wide were milled in 
the trailing edge, and copper inserts were cast into these slots by 
induction heating in a vacuum retort. A copper insert is shown 
in Fig. 3. 

B Rotor Assembly. As reported in reference [1], a single-stage 
turbine was built first. It utilized a rotor disk split in the radial 
plane with conventional axial fir-tree slots to retain the blade. 
The cooling fluid was supplied to each blade from a housing 
located in the center of the split wheel. To connect each blade 
and the housing, two '/s-in-diam tubes were brazed into the 
blade and sealed into the housing by O-rings. This design was 
tested for 100 hours at gas temperatures of 1700 F and was found 
to be adequate. The assembly operation, however, was difficult 
and a new design was sought for the three-stage turbine. 
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bly-simulated blade attachment 


Fig. 5 Test 


As finally developed, the three-stage turbine utilized separate, 
solid rotor disks. Blades were retained individually in the disk 
and the disks were stacked onto a shaft. The use of O-ring seals 
in the coolant passages between the blades and the disks was re- 
tained. These O-ring seals required that the blade be inserted 
radially into the rim of the wheel, eliminating the possible use of 
conventional axial fir-tree slots. In their place a simple Vee slot 
was used as shown in Fig. 4. The blade could be plugged in 
radially and retained by axial pins. Such pinned blade mounts 
were possible because of the effective cooling of the disk and blade 
base. Four pins, 0.083 inch in diameter, were used in each blade. 
Wheel slots were broached and blade bases crush ground, each 
providing one half of the cylindrical pin slots. Pull tests on a 
simulated blade mount indicated strength levels 3.5 times that 
required. A photo of the test pieces is shown in Fig. 5. 

A cross section of a typical rotor assembly is shown in Fig. 6. 
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Fig. 6 Cross section of rotor assembly 
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Fig. 7 Rotor shaft ts—before brazing 
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Fig. 9 Three-stage rotor assembly—radial coolant inlet Fig. 10 Section through stator vane 


4 / sANUARY 1960 Transactions of the ASME 


COOLANT 

oe" AXIAL COOLANT INLET 

j Fig. 8 Three-stage rotor assembly—axial coolant inlet 


Fig. 11 Stator, Stage ll—partial assembly 

The disks were turned from pancake forgings. A central hole 
was bored, and coolant supply and return cavities were machined. 
After finish machining, the blade slots were broached in the rim. 
From the bottom of each slot, supply and return holes were 
drilled into the cavities. These holes were '/s inch in diameter 
and were counterbored to a diameter of '/, inch to a depth of 0.312 
inch. Two studs were brazed into the base of eacn blade, each 
carrying a —1 size O-ring. The studs and counterbored portions 
of the coolant holes provided close dimensional control for good 
O-ring sealing. The fluid pressure built up in the wheel by cen- 
trifugal force is about 3000 psi in the first stage at the O-ring 
and about 4000 psi in the third stage. 

Cooling fluid was transferred to and from each stage by the 
turbine shaft. This shaft was fabricated in two pieces as shown 
in Fig. 7. The center solid member was fluted on its outside to 
provide 12 slots. This center member was copper-brazed into an 
outer cylindrical sleeve. The passages formed between the slots 
and the outer sleeve were used for coolant transfer. A common 
coolant supply was used with six passages, three on each side of 
the shaft. Coolant return from each stage was separated, two 
passages per stage, diametrically opposed. The flutes between 
the adjacent supply passages were locally milled out to ensure 
equal flow to each passage. A photograph of the assembled 
three-stage rotor is shown in Fig. 8. 

Initially, fluid entrance was through radial slots in the shaft at 
a station between the bearings, as shown in Fig. 9. The inlet was 
flooded and outlet throttling was used for control. Flow rates 
for each stage were determined from special tests in which the 
coolant from each stage was collected separately. In the final 
design, coolant flow was controlled by orifices separately sized in 
pairs for each stage. These can be seen in Figs. 8 and 9. 

The efficiency of the radial inlet left much to be desired, and 
necessary flow rates could not be obtained at design speed. Ac- 
cordingly, the radial inlet was abandoned and an axial inlet in- 
stalled on the downstream face of the shaft. The nut that re- 
tained the three rotors to the shaft was modified to serve as a seal 
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Fig. 12 Stator, Stage Il—complete assembly 


COOLING AIR 
DISCHARGE 


Fig. 13 Cross section of turbine—showing cooling air-flow path 


rub ring and coolant pickup. The rotor modified for axial inlet is 
shown in Fig. 8. 

C Stators. The stator vanes were fabricated of sheet metal 
0.025 inch thick. The trailing edges were welded. To provide 
corrosion protection to the lean alloy steel, chromizing was used. 
The internal corrugation was rolled from 0.008-inch mild steel 
sheet and also chromized. A typical stator vane after brazing is 
shown in Fig. 10. The individual vanes were brazed into inner 
and outer shroud rings and baffles added to guide the cooling air. 
Photographs of the second stage stator in process of fabrication 
are shown in Figs. 11 and 12. The cooling air-flow path for each 
stator can be seen in Fig. 13. Air was distributed through mani- 
folds surrounding the turbine casing. The casing was pierced to 
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Fig. 14 Turbine casing 


allow the air flow to penetrate as shown in Fig. 14. The holes 
matched corresponding holes in each stator assembly. Sealing 
between the hot gas and the coolant air was provided by convo- 


luted metal rings designed so that the higher pressure would ex- 
pand the convolution and accomplish the sealing. 


IV Test Program and Results 


The liquid-cooled gas turbine was subjected to extensive tests. 
The 100-hour test of the single-stage split disk configuration was 
reported in reference [1]. To prove the feasibility of the pin- 
mounted blade, a single-stage solid disk was tested. The rotor 
was straddle-mounted and coupled to a dynamometer. The 
rotor setup for these tests is shown in Fig. 15. The straddle- 
mounted configuration was advantageous for instrumentation to 
determine turbine blade temperatures and for control of coolant 
flow. Following the component tests, the turbine was assembled 
as an overhung machine as shown in Fig. 8. 

It had been planned to run the three-stage turbine for a total 
of 500 hours. Mechanical difficulties were such that the program 
was discontinued after 119 hours. At that time the engine needed 
minor repairs to be operable, but it was concluded that the gen- 
eral purposes of the program had been accomplished. 

A review of the design features and operating experience indi- 
cate that the design is fundamentally sound. 

A Rotor Blades. The early tests were plagued with a rash of 
rotor-blade failures, These failures occurred almost wholly in 
the blade tip caps. Alternate methods of brazing were studied 
and proof pressure tests up to three times operating level were 
performed. Blades in the first stage were especially prone to 
these failures. It was found that blade strength was impaired 
by a loss of braze fillet at the time the rotor tip diameter was 
ground. In the second stage, the cap was recessed so that gener- 
ous fillets remained after the tip diameter was ground. This can 
be seen in Fig. 8, in which typical blade tips are shown. The 
major force on the blade cap is from fluid pressure, and the fact 
that blades passed leak and proof pressure tests was taken to in- 
dicate that the shear area was adequate. When operating in a 
turbine, thermal effects caused by temperature gradients could 
easily produce local tension loading on the braze and result in 
failure. 

Toward the end of the program, tip caps were welded in place. 
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Fig. 15 Straddle-mounted Stage | rotor assembly 


Careful technique was required to minimize burn-through and a 
loss of the airfoil profile. In addition, the alloy used in the blades 
was air-hardenable, and welding cracks were frequently en- 
countered. At the conclusion of the program, blade caps were 
retained adequately, but fine leaks detectable only under high 
pressure were found in several blades. The brazed joints in the 
blade manifold plugs and O-ring studs caused no problems. 

B_ Rotor-Blade Attachment. The pin method of rotor-blade at- 
tachment proved very successful. No difficulty was found with 
the feature. Blades were readily interchangeable, and the pin 
joint could withstand a large amount of abuse. Blades could be 
readily changed and at times were changed without rebalancing 
the rotor. Several pins were drilled out with a hand drill and a 
number of holes had drill drifts that caused uneven pin loading, 
but no pin failures occurred. 

C Disk and Shaft Design. The separate solid disks with O-ring 
seals proved to be very practical. Individual rotors could be 
assembled, pressure tested, and then stacked onto the shaft. 
By removing the exhaust collector, the rotor disks could be 
readily removed as shown in Fig. 16. With adequate chamfers on 
the disk bores, the assembly could be readily made without cut- 
ting the O-rings. 

Special techniques were developed for balancing the rotors. 
Blades were weighed individually and selected for assembly. 
For balancing, the entire rotor was filled with water. Filling was 
done by evacuating the system, and then allowing the water to be 
drawn in. Because of the nature of the blade fabrication, some 
rather large out-of-balance conditions were found, and compensa- 
tion would require the removal of a considerable amount of ma- 
terial from the rotor. This was objectionable as it removed the 
nickel-plating and subjected the disk to corrosion. Since, one 
time, holes had been drilled through the disks just below the 
wheel rims for instrumentation, balancing was achieved by in- 
serting tungsten plugs. 

In almost all cases for which satisfactory machine balancing 
was accomplished, no vibration troubles were encountered. As 
was noted, blades were changed on the test stand without re- 
balancing and usually but not always without producing critical 
unbalance. 

D Stators. The stators were damaged by rotor-blade tip caps 
coming loose. Typical damage can be seen in Fig. 17. As re- 
ported in reference [1], air cooling of the stators was generally in- 
adequate. One major failure occurred with a rupture of the weld 
joining the cover plates to the inner shroud ring. This failure was 
attributed to overheating. New plates were welded into place 
after internal baffles were added to provide more effective internal 
cooling in that area. : 
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Fig. 17 Typical damage to first-stage stator 


E Bearings. Conventional antifriction bearings were used with 
good results. Some initial difficulty was encountered with the 
roller bearing as a result of differential thermal expansion. With 
the outer race in a warm housing and the inner race in contact with 
a cold shaft, an increase in radial clearance would take place. 
This resulted in. loss of one bearing, and the assembly was 
changed to provide for a lower initial radial clearance. 

F Coolant Seals. After extensive research and testing, various 
types of coolant seals in the earlier part of the program, standard 
carbon-face Sealol seals were selected for use in the engine. A 
minimum diameter seal was selected to be compatible with the 
shaft design to minimize circumferential rubbing velocity. From 
the manufacturer’s recommendation and the seal test rig ex- 
perience, the flatness and surface quality of the rub ring mating 
face was found to be extremely critical in its effect on seal life. 
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Very close stack-up dimensions were taken on engine as- 

sembly and shims installed to assure proper seal spring compres- 
sion. During actual operation many seal failures occurred. 
These failures were attributed either to foreign matter in the 
coolant, an out-of-flatness on the rub ring caused by distortion on 
stack-up, or excessive engine vibration. Individual seals, which 
were operated under ideal conditions where the mating rub ring 
surface was not distorted and foreign particles did not enter the 
system, displayed excellent wearing capabilities for long life. 
Because of the experimental nature of the program, seal lives of 
50 hours or less were attained. 
r t Coatings. Examination of the turbine 
parts indicated various degrees of effectiveness of the corrosion- 
resistant coatings. Parts in contact with the hot gas had been 
chromized and generally evidences of residual surface chrome 
were apparent. Most of the loss in protection is believed to have 
resulted from excessive vapor blasting. 

Machined parts that had been protected with electroless nickel 
plate showed similar irregular residual protection, again caused 
mostly by handling and cleaning. 


V Performance 

Extensive data were taken to confirm design performance and 
heat-transfer calculations. Aerodynamically, the turbine was 
found to yield satisfactory performance when run with reasonable 
tip clearances. Because of the early failure of blade tip caps, 
honeycomb seals with generous clearances were used except during 
tests to determine performance. 

A Heat Transfer. Heat-transfer data for the liquid-cooled rotor 
indicated that design calculations were more than adequate. At 
the engine design point, a total heat rejection rate to the coolant 
of 187.7 Btu/sec had been estimated. As the measured flow 
rates were higher than design, the theoretical curves of rotor-blade 
heat transfer (Fig. 1) were used with the flow rates from the 
actual test to re-estimate heat-transfer rates. These results are 
presented in Table 4. 


Table 4 Theoretical rotor heat-transfer rates for measured flow 
Coolant flow, Blade Q, Disk Q, 
Btu/sec Btu/sec 
14.9 
10.8 
7.9 
33.6 
Thus the total theoretical rotor heat rejection rate is 261.4 
Btu/sec. The measured rate was 237.0 Btu/sec, 9 per cent below 
the theoretical value. This difference can be explained by the 
fact that the theoretical values were based on the shaft-power 
configuration, whereas the measured values were obtained with 
an air-bleed configuration. The lower pressure for the air-bleed 
configuration would produce a lower heat rejection rate. 
An over-all heat balance is presented in Table 5. 


Table 5 Design point heat balance from test data 


Gas flow rate 6.88 
Turbine inlet temperature 1750 
Turbine outlet temperature 1050 
Turbine enthalpy drop 

Shaft power output 
Compressor enthalpy rise 
Mechanical losses fi per cent) 


Heat rejection to coolant water 
Heat rejection to coolant air 
2.7 Ib/ 


: sec 
249 F AT 
Unaccounted heat losses 


Total enthalpy drop, Btu/see 1368 
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Fig. 18 Rotor-blade trailing edge temperatures 


For the over-all power plant, reduction of the measured per- 
formance data to conditions comparable to design for a shaft 
power machine resulted in an output of 488 hp at a specific fuel 
consumption of 0.987 lb/hp-hr. This does not reflect the power 
required to cool the stators. At the design point, 2.7 Ib/sec of air 
at 75 F and a pressure of 20 psig were supplied. This is equivalent 
to about 175 hp. The higher pressure was required to keep hot 
gas from leaking into the cooling passages. Even with this 
cooling air pressure, evidence of leakage was found in the sooty 
deposits in the cooling air collector. 

D Coolant Contro!. A study was conducted of methods of con- 
trol of coolant flow rates. Although for optimum efficiency, 
constant blade temperatures would be desirable, a system to 
provide the control was considered to be too complicated. It 
was decided to adopt a control system that utilized natural con- 
vection pumping for control. 

As the return fluid is heated, and thus is at a lower density, it 
can be demonstrated analytically that a net pressure rise should 
occur in the rotor. This pressure rise would be proportional to the 
coolant temperature rise and to the square of the speed. During 
single-stage tests, attempts were made to measure this pressure 
rise. At constant speed the gas temperature was raised and the 
coolant flow rate adjusted to a constant value. The resulting 
pressure rise was plotted against temperature rise as shown in 
Fig. 19. Though the convection pumping effect can be seen to 
occur, its magnitude was only a fraction of that anticipated. 

For the three-stage turbine, this convection pumping effect 
was used to provide coolant control. Fixed area orifices were in- 
stalled in the discharge line as shown in Figs. 8 and 9. These 
were sized individually for each stage, but because difficulty was 
experienced with Stage II, it was run with a wide open discharge. 
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Fig. 19 Net coolant pressure rise versus temperature rise, Stage | 


B Blade Temperature Measurements. Data taken during the 
single state tests confirmed calculated rotor-blade temperatures. 
A comparison of blade temperatures with and without copper 
inserts is shown in Fig. 18. Inasmuch as maximum temperatures 
even without the insert were relatively low, it was concluded that 
the inserts were not necessary. 

C Turbine Efficiency. The turbine efficiency calculated from 
the total pressure ratio and the work output was 78 per cent, as 
compared to a design value of 82 per cent. The design value 
corresponded to the higher pressure ratio of the shaft-power 
machine. The measured value was obtained under bleed condi- 
tions. It is felt that the deterioration of the stator, due to 
damage, and the high heat loss were responsible for the lower ef- 
ficiency. 
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Data taken during engine tests demonstrated that, as the gas 
temperature and coolant temperature were increased, the coolant 
flow rate automatically increased. These data are shown in Fig. 
20. 

The increase in coolant flow with coolant temperature rise (and 
gas temperature) indicates that the convection pumping effect is 
occurring. Again, the magnitude is less than expected theo- 
retically. For an increase in the coolant temperature rise from 
70 to 90 degrees, the actual coolant flow rates increased 7.3 per 
cent, as compared to a theoretical value of 16.3 per cent. 


VI Conclusions 


Though the test operation of the liquid-cooled gas turbine was 
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Fig. 20 Cooling system performance, liquid-cooled turbine 


not carried through the scheduled number of hours, most of the 
design details have been adequately proved. The following con- 
clusions are drawn both from the present paper and the earlier 
data of reference [1]. 


1 It was demonstrated that’ a low alloy turbine could be 
operated successfully at temperatures ovér 1700 F. 

2 Adequate corrosion protection can be given to low alloy 
steels for use in liquid-cooled gas turbines. 

3 Available data on heat-transfer coefficients are accurate for 
prediction of cooling losses and blade temperatures. 

4 High temperature brazing can be used effectively for the 
fabrication of components of cooled turbines. 

5 The pin method of blade attachment is effective in produc- 
ing a practical design. 

6 (O-rings were satisfactory for coolant seals between shaft 
and disk and between disk and blades. 

7 Natural convection effects were measured though of 
lesser magnitude than expected theoretically. 

8 The radial coolant inlet to the shaft could not be developed 
satisfactorily. Additional work will be necessary to ascertain if 
it is feasible. The axial coolant inlet was entirely satisfactory. 

9 Face-type coolant seals can be used if cleanliness, installa- 
tion accuracy, and low vibration can be maintained. 

10 The design tested did not wholly solve the stator cooling 
problem. ‘The possible use of ceramic or liquid-cooled stators 
should be considered. 

11 The turbine inlet scroll was never successfully cooled, as 
high alloy liners were used. For high temperatures, turbine inlet 
scrolls are undesirable unless heavy, furnace-type construction 
can be used. An engine arrangement with a straight-through flow 
path is required. 

In summary, the experience with this development program 
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supports the conclusion that liquid-cooled turbines are feasible if 
the added weight and complication can be tolerated. Studies of 
reliability of gas-turbine plants have indicated that the environ- 
mental stress resulting from high temperatures is one of the major 
causes of malfunctions. Liquid cooling, or any cooling, can reduce 
the deleterious effects of high temperatures, but with the penalty 
of added complexity. In the field of small gas turbines, the com- 
plexity does not appear to be warranted for most known applica- 
tions. For larger, stationary machines, the picture may be differ- 
ent. With the use of high alloy steels, turbine inlet temperatures 
over 2000 F appear possible. With high inlet gas temperatures, 
exhaust gas temperatures are also high, and cycles utilizing 
exhaust heat recovery can show attractive theoretical perform- 
ance. 

In conclusion, the authors would like to express agreement with 
the statement made by Dr. von der Nuell in a recent paper [4], 
“Tn general, resort should not be taken to higher working tem- 
peratures until all other methods of improvement have been ex- 
hausted, especially when such methods mean greater safety, 
longer life, greater freedom in materials, and numerous other 
significant advantages.” 
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Effect of Ambient and Fuel Pressure 
on Spray Drop Size 


After a brief review of existing knowledge concerning the effect of ambient pressure on 


spray drop size, results are presented for a swirl nozzle at fuel Ap’s of 25 and 100 pst 
for ambient gas pressures of 0.5, 14.5, and 114.5 psia. The liquid sprayed is diesel 
fuel, the nozzle capacity being 45 gal/hr at 100 psi Ap, and the nominal spray angle, 
80 deg. The photographic method by which drop size was determined is described. 
Curves are presented which show the spatial variation in fuel flow rate, spray-stream 
velocity, and drop size. The notable effects are large drop size and velocity variations 
across the spray stream, and for the total nozzle output, an increase in the drop size as 
the ambient pressure goes from 14.5 to 114.5 psia. Some ramifications of the results 


are discussed. 


1 Introduction—The Need for Drop-Size Information 


| LACK OF KNOWLEDGE concerning the drop size 
distribution of nozzle sprays has long been a stumbling block in 
the path of progress in many fields where nozzles are used and in 
particular in combustion processes. This lack of knowledge 
about fuel spray characteristics makes it difficult to predict ac- 
curately the performance of gas turbine, jet engine, or rocket 
combustion chambers when the operating conditions are varied. 
In these fields much importance is attached to similitude studies 
which would enable one to test at convenient and economical 
engine flow conditions or with scaled-down combustion cham- 
bers; and also to scale up or down in size a successful combustion 
system. Here fuel atomization has an important role and suc- 
cessful scaling of combustion systems depends on knowing the 
behavior of fuel atomizers. 

The effects on drop size of nozzle pressure and fluid prop- 
erties have been studied by numerous investigators, but with few 
exceptions this work has been done at atmospheric ambient 
conditions. A few studies of drop size have been done at elevated 
ambient gas pressures. Lee [1, 2]! reported negligible effect of 
ambient gas density on drop size at a fuel pressure drop (Ap) of 
4120 psi at chamber air density values of 0.31, 0.94, 1.56 lb/ft*. 
He compares his results with those of Sass [3] who reported a de- 
crease in drop size with increasing gas pressure at a fuel Ap of 
4120 psi at chamber air densities of approximately 0.07, 0.38, 
0.77 lb/ft?; and with Woltjen [4] who was unable to detect 
changes in the atomization due to changes in the air density. 
More recently, Giffen and Lamb [5] conducted experiments at a 
fuel Ap of 1800 psi for air densities of 0.076, 0.336, 0.597, 1.116, 
2.155, and 3.20 lb/ft*. They found a decrease in drop size with in- 
creasing gas density. They compare their results with those of Lee 
and Sass and, in addition, Retel [6]. Retel’s data show a de- 
crease in drop size followed by an increase as air density goes from 

0.25 to 2.25 lb/ft’. Thus there appear to be conflicting results 
concerning the effect of ambient gas pressure on drop size; and, 
in addition, for all the work listed, the injection pressures were so 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Gas Turbine Power Division and presented 
at’ the Gas Turbine Power Conference and Exhibit, Cincinnati, 
Ohio, March 8-11, 1959, of Tue AMERICAN SocieTY OF MECHANICAL 
ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, Decem- 
ber 31, 1958. Paper No. 59—GTP-3. 
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high that the resulting Sauter Mean Diameter (S.M.D.) of the 
sprays were all less than 50u. This leaves unexplored a wide 
range of interest in the effects of gas density at lower fuel pres- 
sures, an example being the gas-turbine range at fuel pressures of 


25 to 400 psi. Work has already been done at this laboratory to 
find the effect of ambient and fuel pressure on fuel distribution 
[7]. 

Now the objective is to study the drop size distribution of fuel 
sprays at fuel Ap’s of 25 to 400 psi, while varying ambient gas 
properties. This paper reports the results so far obtained in this 
endeavor, namely, drop-size data for a 45-gph—80-deg nozzle at 
fuel Ap’s of 25 and 100 psi, and ambient gas pressures of 0.5, 
14.5, and 114.5 psia. The ambient gas used was nitrogen at room 
temperature and the fuel sprayed had the following properties: 
Specific gravity at 60 F = 0.84; surface tension at 77 F, 30 
dynes/cm; viscosity at 77 F = 2.4 centistokes. The nozzle is of 
the pressure atomizing centrifugal type which has found wide 
usage. Fig. 1 shows the construction of such a nozzle. The 
evaporation rate from the spray is negligible. 


Nozzle Spray 

Cone Angle 

Fuel Spray 
Sheath 


Swirl Orifice 
Chamber 
Swirl Slot 
Element 


Flow 
Fig. 1 Typical centrifugal nozzle 
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I! Description of Method and Apparatus 


The method utilized to obtain drop size will be briefly described 
here. In reviewing the techniques that have been utilized for 
drop-size studies, it was seen that no one method had been de- 
veloped to the point where it had overpowering advantages over 
other techniques. Consequently the choice of method has de- 
pended on the objective at hand. One might list a set of charac- 
teristics of an ideal method: 


a Should not disturb the fuel spray pattern or atomization. 

b Should provide a rapid means of sampling. 

c Should provide a rapid means of counting the samples ob- 
tained. 

d Should have good size distinction over the whole range of 
drop sizes to be measured and diameter size intervals should be 
small enough to provide meaningful data in the less than 50u 
range. 

e Should permit the variation of (a) liquid and (b) ambient gas 
properties. 

f Should permit one to obtain drop size distribution i in space, 
i.e., in a given volume of the spray. 

g Should permit one to obtain the nozzle drop size distribution 
produced by the nozzle in a given time. 


Since one of the prime objectives was considered to be a study 
of ambient gas effects, this implied spraying within a sealed tank 
which could be pressurized or evacuated. Because of this con- 
sideration, the photographic method was chosen as a means of 
extracting data (samples) without having ready access to the tank 
interior. Thus the primary requisites of (b) and (e) were met. 
This method is also strong in categories (a) and (f). Its major 
weak point is (c), speed in counting, since a manual counting 
procedure is used, Another weak point (g) is that the temporal 
drop size distribution for a nozzle must be obtained indirectly by 
consideration of local spray velocities and mass flows. 

Existing photographie techniques were found to have two 
major weaknesses and development work was concentrated on 
these: 


1 Fouling of optical surfaces by spray 
2 Uncertainty in depth of field determination 


The first problem was satisfactorily solved by developing a dual- 
shutter system which exposes the optical surfaces for very short 
intervals so that for one set of photographs (24) the surfaces are 
exposed a total of approximately 0.12 sec. The second problem 
which has long been the major drawback of photographic tech- 
niques was met by utilizing an optical system which produces dif- 
fraction patterns around the drop images on the photographic 
film. The diffraction theory for Fresnel patterns indicated that 
the width of the first diffraction ring around the drop image 
could be used as a criterion for determining the position of drops 
in the camera field. This procedure was verified experimentally 
so that a depth of field could be specified for any value of fringe 
width. i 

A drawing of all the essential components in the system is shown 
in Fig. 2. The camera, light source, light source tube, and shut- 
ters are located on a common axis. The axis of the nozzle holder 
is located so that the nozzle axis, with the nozzle at the vertical 
position (i.e., zero degrees) will pass through the camera axis at 
the location of the plane of best focus. 

The nozzle holder has a worm and gear arrangement which 
allows the nozzle (A) to be rotated about an axis which passes 
through the nozzle tip and is perpendicular to the nozzle axis and 
the camera axis. With the tank under pressure, by rotating 
nozzle positioner (B) to a particular setting of a graduated dial 
the nozzle may be inclined at any desired angle between zero and 
56 deg. @ is the angle between the nozzle axis and the vertical as 
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shown in Fig. 2. The nozzle could be set at any one of several 
heights above the camera axis by choice of the extension (C). 

The optical path may be traced along the camera axis as fol- 
lows. The light originates at the Xenon flash tube (E) and pro- 
ceeds through an adjacent pinhole and condensing lens (F). 
Thence the path is through the window (G), light tube (H), light 
shutter (I), shield tubes (J), and fuel spray camera shutter (K), 
camera window (L), lens (M), post camera grid (N), ending at 
photographic film (O). 

The camera lens (M) is threaded into the tip of the conical sec- 
tion of the camera and in front of the lens is a plexiglas window 
which is held against a flange thus forming a pressure-tight seal. 
The compound lens is 1.25 in. in diameter and has a focal! length of 
3.5 in. At the back of the camera a wire grid is installed directly 
in front of the film holder. This grid is made up of 0.0005-in- 
diam wire and produces a fine-lined silhouette on each photograph 
which assists greatly in the counting procedure. 

It was found that suitable photographs could not be obtained 
when any light-transmitting surfaces were exposed to impinge- 
ment of spray particles. Allowing fuel to completely wet these 
optical surfaces was not satisfactory because of the variation in 
thickness of the fuel film as particles impinged on it. This dif- 
ficulty led to the use of a shielding system consisting of two 
Alphax shutters, one on the camera side (K) and one on the light 
side (I). These shutters are of standard leaf type with strength- 
ened spring action to overcome the surface tension forces of the 
fuel wetting the leaves. By setting the shutter diaphragms down 
to approximately 0.4 in. in diameter and for a shutter exposure 
time of 4/99 second, it is possible to obtain a set of 24 or more 
photographs without cleaning of the exposed optical surfaces. 
The total exposure time of the optical surfaces for 24 photographs 
is then 0.12 sec. Also this allows unlimited time for steady spray 
conditions to be reached, while the shutters remain closed. The 
shutters are tripped simultaneously from outside the tank by - 
action of a solenoid-operated valve which supplies gas pressure 
to a cylinder-piston arrangement, Fig. 2 (P), within the tank. 
The light source is flashed by the synchronizing contacts of one 
of the shutters, through the synchronizing cable, Fig. 2 (Q). 

Under certain spray conditions the spray is so dense that ad- 
ditional shielding is needed to reduce the length of the illuminated 
path in the spray in order to obtain satisfactory photographs. 
Two tubes approximately 1/2 in. in diameter are used for this 
purpose, being attached to the shutters, see (J) in Fig. 2. Be- 
tween the tubes a gap remains through which the spray passes. 
At the end of each of these shield tubes is a thin projection which 
prevents particles which bounce on the tube from being carried 
down into the gap. In obtaining some of the earlier data this 
shielding was accomplished by two truncated conical sections held 
between the shutters. To reduce recirculation and bouncing of 
the fuel drops from the tank and camera walls, matting was used 
as shown in Fig. 2 (R). 

The spray tank has an inside diameter of 18 inches and the 
height of the portion shown in Fig. 2 from flange to flange is 30 
inches. Below the section shown in Fig. 2 is additional tank vol- 
ume. The optical effects of distortion of the plexiglas camera 
and light windows, Fig. 2 (G), (L), due to varying tank pressure, 
were found to be negligible in the range of 0.5 to 114.5 psia. This 
was checked by examining photographs of a ronchi grating taken 
at various tank pressures. The camera, light source, and nozzle 
holder flanges are all doweled to the tank to insure that proper 
alignment is maintained. 

The fuel system and gas system for pressurizing the tank are 
essentially the same as described in reference [7]. 

In taking photographs of the spray it is desirable to know the 
fuel flow-rate distribution in space as a function of @ (nozzle 
angle) and R (the distance between nozzle tip and camera axis). 
To obtain this distribution an open-end sampling tube is placed 
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Fig. 2. Drawing of drop-size apparatus 


vertically directly below the nozzle with its tip at the camera axis. 
The flow distribution is then obtained by rotating the nozzle 
(varying @) and recording the flow into the sampling tube for 
timed ‘intervals in the manner described in reference [7]. Once 
the flow-rate distribution is known, considerable time is saved by 
photographing only regions where the flow rates are not negligible. 
The flow-rate distribution is also required in order to obtain 
values of spray-stream velocity which are necessary for converting 
the space drop-size distribution into a temporal value, that is, the 
distribution produced by the nozzle in a given interval of time. 

Eastman Kodak Contrast Process Ortho film in the 4 x 5-in- 
sheet size is used to obtain the drop photographs. For this film 
the manufacturer gives a resolving power of 145 lines/mm. On 
the basis of a camera magnification of 9.12 this resolving power in 
relation to the object photographed is 1320 lines/mm or 1.32 
lines/u. The film has a light sensitivity only in a fairly narrow 
band in the blue-green region of the spectrum. 

The illumination for the photographs is produced by the dis- 
charge of capacitors across a Xenon-filled quartz tube. The light 
source was adapted to give either of two duration times, 2 or 5u- 
sec, with the shorter flash having less intensity. 

Drop sizes were tabulated by an operator who examined the 
negatives on an optical comparator at a magnification of 20/1, 
thus providing an over-all magnification of 182/1. A vexing 
problem which arises in the counting procedure is the question 
of how many drops must be counted before one can be assured 
that he has an accurate representation of the drop-size distribu- 
tion for any given set of physical conditions. The procedure used 
was to apply the t-test familiar in statistics considering each 
negative as an individual sample of the spray drop population. 
The application of the t-test to obtain confidence limits is similar 
to the technique using the standard normal distribution but is 
applicable for small samples. The t-test is applied to obtain con- 
fidence limits on the mean of the drop volume found in each 
sample (i.e., each negative). The drop volume is the pertinent 
quantity of interest which is most sensitive to changes in drop 
diameter determination and hence it was chosen a's the quantity 
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to be used in the t-test. It occurs in the determination of the 
Sauter Mean Diameter (S.M.D.) which is of interest in combus- 
tion applications. The 8.M.D. = 2nj,d,°/Zn,d;2 where d, refers 
to the drop diameter and n, is the number of drops found of that 
size. The drop count from each negative is tabulated in histo- 
gram form, the size intervals being 5u. From these histograms, 
then, any required characteristic diameter can be calculated and 
cumulative volume or number curves can be plotted. 

The operator who does the drop counting in any photographic 
method has to make decisions on which drops to count, i.e., decide 
which drops are outside a predetermined depth of field. It is 
known that an operator counting drop images that appear to have 
the same degree of sharpness will count large drops over a greater 
depth of field than will be covered for the count of smaller drops, 
thus giving a count biased toward the large drops. Manson, etal., 
[8], carried out a calibration to correct this effect and found that 
in going from drop diameters of 12.5 to 150u the depth of field 
over which drops were counted varied from 3.7 to 8 mm. York 
and Stubbs [9] compensated for this effect by the use of images 
of varying diameter having standard sharpness which an opera- 
tor could compare with the drop images being counted. In our 
development work it was found that the formation of Fresnel 
rings or fringes around a drop image could be utilized to provide 
a measure of the location of the drop in the field. The width of 
the first fringe was then used as a criterion for deciding whether 
the drop should be counted. The essential elements needed to 
produce Fresnel] diffraction patterns are a point source of mono- 
chromatic light, an obstacle, and a screen; or, if white light is 
used, the screen must be receptive only in a narrow band of the 
spectrum. The latter is the case with the system used here. 
The receiving screen, which is actually the film, receives the image 
formed at the plane of best focus. When an obstacle is located on 
the light-source side of this plane, a real fringe pattern is formed at 
the plane; when the obstacle is located between the plane and 
camera, a virtual pattern is formed at the plane. By use of diffrac- 
tion theory the Fresnel diffraction patterns were calculated for 
several cases of interest. In addition, comprehensive calibration 
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work was done over a drop size range of 35 to 500u using opaque 
circular spots and fuel droplets suspended on fine wire. This 
calibration work resulted in the curve of Fig. 3 which shows the 
fringe width as a function of distance from the plane of best focus. 
If it is desired to count drops over a depth of field of 1280, for 
example, then all drops having a fringe width of 27.5 or less 
would be counted. Fig. 4 is a print of a negative taken of an 
actual spray, with the fringes clearly visible around the drops. 

The Rayleigh limit of resolution for the shutter-camera system 
was calculated to be 3.14. This is a limit imposed by diffraction 
at the lens aperture. Calibration work with fine-lined reticles 
and diffraction gratings indicated an actual resolution limit of 3 to 
4u. One should be mindful of this limit in assessing the accuracy 
of the system. 

On the calibration work which preceded the spray photography, 
it was found that the camera magnification factor for droplet or 
spot photographs was greater than the nominal value of 9.12 by 
amounts which ranged from 10 per cent for 5-10y size to 0.5 per 
cent for the 500u size. While this represented small absolute 
changes in measured drop size, corrections could be easily incor- 
porated into the counting procedure, and this was done. 


Hl Results of Drop-Size Experiments 


The system used to identify a particular spraying condition can 
be illustrated by this example. The notation 45/80°|25|14.5|35° 
denotes the following situation: A nozzle rated at 45 gal/hr flow 
(on diese] fuel at 100 psi Ap) and of nominal spray angle equal to 
80 deg, spraying at a fuel pressure drop of 25 psi into gas which is 
at a pressure of 14.5 psia, with the nozzle inclined at an angle 0 = 
35 deg. Fig. 2 shows the angle 6. 

The quantities of interest (m, d,, u, S.M.D., N) are tabulated 
for various values of @in Table 1. In addition, m, u, and 8.M.D. 
are plotted versus @ in Fig. 5 and Fig. 6. Some explanation of 
these quantities is needed. The quantity m represents the 
volumetric fuel flow rate at some point in the spray stream and is 
obtained by the sampling tube procedure previously described. 
The quantity d, represents the average over the photographs 
counted of the total volume of fuel drops in a photograph of 
standardized size (10 in.*), for some given spraying condition. 
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Since the depth of field is known and the area of the photograph is 
fixed (10 in.*), this represents the fuel volume in a given space 
volume (obtained essentially instantaneously) at some chosen 
point in the spray. The quantity u is the velocity of the droplet 
stream. This can be obtained once the following four quantities 
are determined: m; 4,; Az = cross-sectional area of sampling 
tube = 0.170 cm*; V, = space volume represented by one 
standard photograph. This latter quantity, V, = area of photo- 
graph X depth of field/(camera magn)*. The depth of field is 
1280u, camera magnification is 9.12, and photograph area = 10 
in?. Thus V, is found to be 0.0993 cc. From these: 


m 
u = 0.000547 rari (ft/sec) 


u = 32.0 X 0° 


where m is in cc/min and 4, isin cc. The next quantity in Table 
1, S.M.D., refers to the Sauter Mean Diameter of the spray, a 
characteristic diameter of special value in combustion applica- 
tions. The S.M.D. is the diameter of a drop which has the same 
ratio of volume to surface area as that of the actual spray itself. 
It is determined from the drop count and is equal] to a summation 
over the spray of 2n,d,3/Zn,d,*, where d; refers to the drop 
diameter and n, is the number of drops found of that size. As 
shown in Table 1 and in Figs. 5 and 6, 8.M.D. represents an 
instantaneous or space distribution of fuel drop size without re- 
gard to the mass flow of fuel at any point. The distance between 
the nozzle tip and the camera axis for the data presented is 7.24 
inches. The last quantity shown in Table 1 is NV, the drop number 
density, given in number of drops/in®. 

Several points of interest appear in Figs. 5 and 6: (1) The shift 
of fuel flow distribution with changing ambient pressure is again 
evident as was pointed out in reference [7]; (2) except for the 
45/80°|100|114.5 case, there is a decided decrease in drop size 
at the inner portions of the spray. Thus a segregation of par- 
ticles appears to be taking place. This process may be attributed 
to the same induced ambient gas crossflow which causes the 
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Fig. 4 Print of spray photograph 


Table 1 Spray variables 
45/80° nozzle, 25 psi, fuel Ap 


Tank pressure, psia 0.5 14.5 ~ 114.5- . 
6, deg 44.5 28.5 35 41.5 3.75 11.25 18.75 26.25 
m*, cc/min 1 95 0.334 0.835 0.366 0.47 0.9 1.29 0.5 Pa 
G,, cc X 10~* 0.563 0.0753 0.144 0.855 0.80 1.51 6.26 3.85 
u, ft/sec iA 14.2 18.5 1.37 1.88 1.91 0.66 0.42 
S.M.D., » 206.5 90.8 164.5 196.5 101 170 229 247 
N, no. drops/in.? 25820 31550 10940 4600 2130 3330 1860 1520 
45/80° nozzle, 100 psi, fuel Ap Pi 
Tank pressure, psia 0.5 14.5 - 114.5 ~ 
0, deg 38.5 41.5 44.5 12 24 30.5 35.5 40.5 2.5 7.5 12.5 
m*, cc/min 1.10 3.8 1.5 0.4 0.6 1.05 2.6 0.65 19.0 12.0 2.65 
G,, cc X 10-4 0.0865 0.395 0. 260 0.22 0.46 0.553 0.639 0.218 1.70 1.875 2.50 
u, ft/sec 40.6 30.8 18.45 5.8 4.1 6.1 13.02 0.954 35.8 20.5 3.39 
S.M.D., u 57.0 149.6 104.5 40.7 54.0 80.7 94.9 72.3 103 114 108.4 
N, no. drops/in.* 131800 228000 96200 48900 51000 32270 75520 . 6370 18990 26550 11650 


*m values obtained from m versus @ curve. 
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spray flow pattern to close in as tank pressure is increased. This 
segregation of particles was also noted by Manson, et al. [8]; (3) 
there is considerable variation of stream velocity across the 
spray. In particular note the increase in velocity at the 114.5 
tank pressure condition as @ decreases, i.e., as one moves toward 
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the center of the spray. At the atmospheric condition (14.5 psia) 
the maximum velocity values coincide with the peaks in the flow- 
rate curves. The maximum velocity for the 45/80°|100|0.5 case 
would appear to be located inside the flow-rate peak. Since 
photographs were obtained for only one value of #, no u or 8.M.D. 
curves can be drawn for the 45/80°|25|0.5 condition. 

The S8.M.D. values of Table 1 and Figs. 5 and 6, as stated pre- 
viously, represent measures of the drop size existing at various 
points in the spray. As opposed to this, a second type of 8.M.D. 
value may be obtained which represents the total drop output of 
the nozzle, for any given spraying condition. This is obtained 
from the data on hand by weighting the numerator and de- 
nominator of the spatial S.M.D. values by the stream velocity at 
the sampling point and by a representative stream-flow area for 
that sampling point. If we denote the spatial 8.M.D. value for 
a particular spraying condition and value of 6, as 

a 


8.M.D. = 


then the “nozzle” S.M.D. is 


which is summed up for the number of sampling points (i.e., @ — 
values) at any spray condition. A, is the representative flow area 
at 6;. The u-weighting accounts for the velocity variation in 
the spray stream and the A-weighting the variation of flow area 
as different rings or bands of the spray stream are considered. 
The flow areas, A;, are taken to be ring-shaped portions of a 
sphere of radius R (where R = distance from nozzle tip to camera 
axis). The area which includes the 9 = 0 point would be a 
spherical segment. When these calculations are carried out the 
following values of (S.M.D.),y result: 


Table 2 (S.M.D.)y values 


Tank pressure (psia) 
Fuel Ap, 25 psi 
Fuel Ap, 100 psi 


These values are plotted in Fig. 7.. The striking thing here is 
the increase in size indicated in going from 14.5 to 114.5 psia 
ambient pressure. From the aspects of spray breakup a con- 
tinuing decrease in drop size with increasing ambient pressure 
would be expected, since the drag force on a drop increases with 
increasing density. Thus as ambient pressure rises, a decrease 
would be expected in the critica] drop size, i.e., the maximum size 
that can withstand breakup as is indicated by Hinze [10] and 
Lane [11]. The explanation for the anomalous increase in drop 
size of Fig. 7 may lie in increased coalescence of the spray droplets 
as the ambient pressure is increased. 


250 T 


10 100 
Ambient Pressure, psia —> 


Fig. 7 (S.M.D.)v versus ambient pressure 
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It can be seen from the flow distribution curves of Figs. 5 and 6 
and from the curves of reference [7] that the action of the induced 
gas flow is such as to increase the probability of drop collisions 
especially at the higher ambient pressures where the fuel spray 
has been collapsed into a small volume. It appears reasonable to 
postulate that the curve shapes in Fig. 7 result from the following 
processes operating simultaneously. In the breakup process, in- 
creased ambient pressure tends to produce smaller drops. But the 
coalescence of drops also is increasing and becomes dominant in 
going from 14.5 to 114.5 psia in pressure. 

A considerable effect of drop coalescence is reported by Dunskii 
[12] for the case of an air-atomizing nozzle. He also finds ana- 
lytically the change in drop size distribution for orthokinetic 
motion of a system of drops which move through a gas with dif- 
ferent velocities. Turbulence or crosswise gas flow is not con- 
sidered. 

The accelerating force necessary to make a drop travel along a 
circular are toward the center of the spray is 


T u? 
— 


and the drag force exerted by induced gas flow, in the same direc- 
tion as F,, is 


1 wd? 
F, = Capyu,? (=) 


Here d = drop diam, p = fluid and p, gas density, « = drop 
velocity, r = radius of curvature of path, u, = gas velocity rela- 
tive to drop in direction of F,, C, = drag coefficient based on 
u,. Then 


3p u,\? r 
« 4 p (*) d 


From this it can be seen that the dependency of this ratio on d 
will vary essentially from 1/d* to 1/d depending on whether C, is 
evaluated for laminar or turbulent flow. Thus the smaller the 
drop the more susceptible it is to entrainment toward the nozzle 
axis, which is borne out by the data obtained. 

Various experiments have fitted data to the expression 
(S.M.D.)y = KP* at some constant value of ambient pressure 
(usually atmospheric). Here K = a nozzle constant and P is the 
fuel pressure drop. Although only two points are available at 
each tank pressure yet it is of interest to find the values of z for 
the data of Fig. 7. These are, for tank pressures of 0.5, 14.5, and 
114.5 psia, z = —0.48, —0.50, and —0.45, respectively. 

The data in this paper are based on a count of 31,375 drops. 
The drop size distributions obtained are presented in the form of 
curves of cumulative per cent under size by number plotted on 
arithmetic probability paper. These are shown in Figs. 8 and 9. 
A normal] distribution of drop size would plot as a straight line on 
this paper. Attempts were also made to fit the data to the Rosin- 
Rammler equation and the log-normal] distribution, but curves 
were obtained which varied considerably from straight lines. 
The advantages of working with a straight-line plot of the data 
are great but unfortunately spray drop-size data seldom conform 
to any of the drop-size distribution functions that have been ad- 
vanced. For each spray condition in Figs. 8 and 9, a combined 
curve is shown which represents the drop-size distribution of the 
total nozzle output at that condition. These are obtained by 
weighting the individual distributions at each value of 6 by the 
appropriate stream velocity and flow area corresponding to the 
6-value. This procedure is similar to that used to obtain the 
(S.M.D.)y values, the velocities and flow areas used being identi- 
eal. Of course, for the 45/80°|25|0.5 condition only one curve 

appears because only one sampling point was used. Many of the 
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distributions in Figs. 8 and 9 and, in particular, the combined dis- 
tributions, appear to be of the type that could be satisfactorily 
represented by two straight lines. This suggests that the atomiza- 
tion process is occurring in such a way as to produce the sum of 
two distinct distributions. Fig. 10 is drawn to show how curves 
of the type seen in Figs. 8 and 9 can occur as the sum of two nor- 
mal distributions, possibly added in varying proportions. The 
norma] distributions used in Fig. 10 are one of mean = 30, stand- 
ard deviation = 10, and one of mean = 100, and standard devia- 
tion = 100. These are designated as (30, 10) and (100, 100). 
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IV Discussion of Errors 


A discussion of possible errors in the results presented is in 
order. When sampling from a particle stream is attempted there 
exists the danger that the sample will be biased. Here the 
sampling structure can be considered as a cylinder with a section 
sliced out so that the sample passes down through the cutout 
section, see (J) Fig. 2. The available information dealing with 
flow around solid objects indicates that a decrease in sampling ef- 
ficiency occurs as: (1) Flow velocity decreases; (2) sampling 
structure size increases; (3) as drop size decreases. Of secondary 
importance is the decrease in sampling efficiency as the ambient 
pressure increases. Keeping these effects in mind, ‘no trend in the 
drop size results was discerned which could be attributed to 
sampling error. 

A second source of error arises in the calculation of the fuel 
volume per spatial volume, which in turn affects the stream ve- 
locity value. This is due to an “edge effect’’ arising from the fact 
that drops at the edge of the sample volume (i.e., just barely 
within the depth of field) will have a varying amount of volume 
projecting outside what is considered to be the sample space or 
sample volume. Since all the drops counted are considered to be 
within the sample space, it is seen that an error is present which 
would increase as the drop size encountered increases and as the 
depth of field is decreased. This error has been found to be a 
linear function of depth of field and drop diameter. For the 
depth of field of 1280, the calculated error goes from zero at drop 
diameter = 0 to 7 per cent at 500u. Because application of this 
correction would be rather difficult and because of its limited size, 
the correction was not applied in the calculation of stream ve- 
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locity. However, it is essential to establish its magnitude and, if 
a smaller depth of field is specified, the correction might have to 
be considered. 

For any given spray condition, d, represents the mean of the 
total drop volumes per sample (negative). 4, figures in the calcu- 
lation of stream-velocity and thus in the (S.M.D.)y values. The 
deviation of the sample mean 4, from the actual mean volume 
can only be estimated. The statistical procedure is first to assume 
some distribution function for the quantity to be evaluated, which 
in this case is the drop volume per negative. For the t-test a 
normal distribution is assumed. The expected error at any chosen 
confidence level is then a function of the size of the sample and the 
scatter of the individual a, values. In our case the size of sample 
refers to the number of negatives counted. The 8.M.D. accuracy 
also depends on the magnitude of the drop count but to a lesser 
extent than d, since the S.M.D. value is not as sensitive to varia- 
tions: in the number of drops occurring in individual negatives. 
The percentage errors in the dG, or 4 figures are listed in Table 3, 
calculated at the 95 per cent confidence level, according to the 
t-test. 

A review of the change in per cent error at each condition as the 
sample size was increased indicates that the figures of Table 3 
are on the pessimistic side. The stream velocity figures would be 
most affected by any variations in d,. In the calculation of 
(S.M.D.)y any variation in @, would be partially compensated by 
a corresponding opposite variation of the local 8.M.D. value. 


V Additional Discussion of Results 
Some implications of this work will be presented here. From 
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45/80° 25 0.5 44.5° 21.8 
45/80° 25 14.5 28.5° 39.5 
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45/80° 25 14.5 41.5° 42.0 
45/80° 25 114.5 3.75° 31.8 
45/80° 25 114.5 11.25° 46.5 
45/80° 25 114.5 18.75° 8.7 
45/80° 25 114.5 26.25° 15.4 
45/80° 100 0.5 38.5° 42.3 
45/80° 100 0.5 41.5° 40.0 
45/80° 100 0.5 44.5° 42.2 


the drop-size segregation which occurs it is evident that photo- 
graphs of a spray which show only the outer portions of the spray 
stream will be misleading and will generally give a drop size 
measure larger than the actual. An appreciable coalescence of 
drops has been postulated to occur in these sprays. If this is the 
case then any drop-measuring methods which prevent or hinder 
this effect (such as the molten-wax method) should be used with 
caution. 

At high fuel and ambient gas pressures the coalescence effect 
appears to be fully as important as the spray breakup process 
and deserves increased study. For example, this effect may help 
to explain the well-known increase in drop size as the nominal 
spray angle of a nozzle is decreased, since this decrease provides 
increased opportunity for drop collisions. 

Translated into the area of combustion chamber design, it is 
hoped that this work will help to explain anomalous behavior of 
combustion systems as spray conditions vary. Also with increas- 
ing knowledge of environmental effects on sprays, the designer 
will have at hand additional necessary tools for rational design 
and scaling of liquid-fired combustion systems. 
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Table 3 Per cent error in a, at 95 per cent confidence level 


Per cent error 


condition 


45/80° 100 14.5 12° 18.2 
45/80° 100 14.5 24° 32.2 
45/80° 100 14.5 30.5° ae 
45/80° 100 14.5 35.5° 14.7 
45/80° 100 14.5 40.5° 16.8 
45/80° 100 114.5 2: 5° 15.6 
45/80° 100 114.5 29.4 
45/80° 100 114.5 12.5° 20.3 
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Axial-Flow Turbines, Including Design and 
Test Results of a Typical Axial-Flow Stage 


A method is presented to determine the shape of stationary nozzle blades and rotor blades 
for an axial-flow-type turbine in a generally consistent manner based on the concept 
of aerodynamic blade loading. The mean blade load is a typical design parameter 


which predominantly determines the blade curvature. It depends in particular on the 
rate of change of momentum across the blade row. By applying the design method, 
airfoil shapes are obtained which satisfy the momentum requirements regardless of 
what blade-load distribution is assumed as long as the mean blade load remains con- 
stant. A specific application of the design method is described and test data are pre- 
sented which show that good agreement between design goal and test results was achieved. 


Bis IN THE DEVELOPMENT PROGRAM for small gas 
turbines, undertaken by the authors’ company, the need was 
recognized for a consistent method to determine the blade shapes 
in turbines and compressors of axial and radial-flow type. The 
method had to be complete enough to define the correct blade 
shape uniquely for any given application, yet it had to be clear 
and simple enough to permit the average turbine designer to 
evaluate the basic problems involved without applying compli- 
cated mathematical procedures or incomprehensible experience 
factors. It appeared to be desirable to combine experimentally 
determined blade-design information with a simple mathematical 
formulation of the flow problem to obtain a clear picture of the 
flow pattern. 

Extensive testing by the authors’ company on airfoil cascades 
furnished some basic information on pertinent design parameters 
such as blade loadings, angles of incidence, blade-overturn angles, 
pitch-chord ratios, aspect ratios, and soon. The flow problem was 
formulated according to the classical fluid-dynamics theory. The 
derivation of the basic flow equations is given in Appendix IT. 
These equations, together with empirically established design 
parameters, may be used to calculate the blade profiles for tur- 
bines and compressors of either axial or radial-flow type. De- 
tailed procedures to cover all these specific applications of the 
same basic theory have been worked out and used successfully by 


Contributed by the Gas Turbine Power Division and presented 
at the Gas Turbine Power Conference, Cincinnati, Ohio, March 8-11, 
1959, of Tae American Society OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
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the Society. Manuscript received at ASME Headquarters, December 
31, 1958. Paper No. 59—GTP-4. 


the authors’ company. Applying such a strict mathematical 
treatment to the three-dimensional compressible-flow pattern of 
a radial-type turbine or compressor involves, of course, a more 
complicated procedure than applying it to a two-dimensional 
flow pattern of axial-flow-type turbines and compressors. 

It should be mentioned here that the same analytica] approach 
to blade-shape calculations has been used by other designers. A 
notable example is presented in a paper by Concordia and Dowell,! 
where development work on impellers for centrifugal air compres- 
sors is described. 

In the following discussion we shall limit ourselves to the case 
of the axial-flow-type turbine. 


Description of Method 


Before the blading design for any particular turbine stage can be 
started the basic design parameters must be established. First, 
the flow rate must be determined, based on power output, turbine- 
inlet conditions, turbine-exit conditions, and expected over-all 
stage efficiency. Next, an appropriate flow pattern is selected. 
For theoretical considerations, a free-vortex flow pattern is 
usually selected because it requires the least amount of kinetic 
energy for a given flow rate and therefore represents the most 
stable flow condition. In a free vortex the tangential-velocity 
components of the stream particles are inversely proportional to 
their radial distance from the center line and the axial-velocity 
components are constant across the whole flow area. A free- 
vortex flow pattern is maintained across the rotor if equal 
amounts of work are taken out at any radius station. 

1C. Concordia and M. F. Dowell, ‘‘Analytical Design of Centrifu- 
gal Air Compressors,” Journal of Applied Mechanics, vol. 13, 
Trans. ASME, vol. 68, 1946, p. A-271. 


Nomenclature 
r = radial co-ordinate, ft U = circumferential velocity of rotor Subscripts 
z = axial co-ordinate, ft blades, ft/sec 0 = nozzle inlet 
o= angular co-ordinate, radians V = absolute flow velocity, ft/sec 1 = noskle exit—also rotor inlet 
y = r@, circumferential co-ordinate, ft |W = flow velocity relative to moving 2 . 
3 3 = rotor exit 
w = dd/dt, angular velocity, radians blade, ft/sec ; 
& = static pressure 
per sec p = pressure, lb/ft? ; 
V, = rw, circumferential velocity, ft/sec  y = specific weight, lb/ft* total pressure 
a = dw/dt, angular acceleration, ra- p = Y/Q, specific mass, slugs/ft* rel = relative hes a blades 
dians/sec* g = 32.2, gravity acceleration, ft/sec? r = radial direction 
ra = dVu/dt, circumferential accelera- J = 778, work-heat-equivalent, ft-lb/ u = circumferential direction 
tion, ft/sec? Btu z = axial direction 
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The degree of reaction that can be built into the turbine stage is 
dictated to a large extent by stress considerations. A high degree 
of reaction is desirable if the highest possible efficiencies are to be 
attained. A high degree of reaction, however, requires high 
blade speeds and causes high rotor stresses. In addition to this, 
the blade temperatures in the rotor increase with increasing de- 
grees of reaction for any given turbine-inlet temperature. This 
further aggravates the material requirements. If a maximum 
power output per stage is desired, then it is necessary to design 
the blading for zero reaction at the blade-root station to keep the 
stresses within acceptable limits. Under these conditions about 
20 to 30 per cent reaction is obtained at the tip station of a typical 
small gas-turbine wheel. The best compromise between high 
power output per stage and high efficiency is usually obtained by 
designing for about 15 to 20 per cent reaction at the blade-root 
station, whereby about 50 per cent reaction at the blade 
tip station is obtained. 

A typical velocity diagram, a corresponding entropy diagram, 
and definitions of the most commonly used performance parame- 
ters for axial-flow turbines are shown in Appendix I. 

After the basic design parameters are established, aerodynamic 
design of the blading can be started. It should be noted that 
velocity diagrams usually apply to a mean streamline, located 
about halfway between two adjacent blades. In the immediate 
vicinity of the blades, the flow directions will vary to some extent 
from those of the mean streamline. The deviations of the concave 
side of the flow boundary from the mean streamline at entrance 
and exit are defined, respectively, as angle of incidence and angle 
of overturn. Design information about these angles must be ob- 
tained from cascade tests. Fig. 1 gives typical design values as 
function of pitch-chord ratio. 

For practical purposes, it is convenient to arrange the velocity 
diagram as shown in Fig. 2, where the total velocity change in 
the tangential direction across the rotor blades is represented by a 
single vector B. Vector B together with the flow rate determines 


INCIDENCE AND OVERTURN ANGLE, DEGREES 


5 6 9 10 
PITCH-CHORD RATIO ty, 


Fig. 1 Typical design values for incidence angle and overturn angle 
of nozzles and rotor blades, based on cascade tests by authors’ company 


Fig. 2 Typical velocity diagram for cxial-flow turbine stage 
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the momentum change in the rotor and fixes the torque which 
the gas exerts on the rotor. 

If the velocity diagram, Fig. 2, were used directly to determine 
the shape of the rotor blades, a blade would result that would 
provide neither an incidence angle at the inlet nor an overturn 
angle at the exit. The gas would not be turned enough by the 
blades, the required velocity changes would not be obtained and, 
as a result, the power output of the turbine stage would fall short 
of the requirements. To remedy the situation, a corrected velocity 
diagram is drawn up that takes into account. an overturn angle. 
By doing this, the vector B is increased which means that a larger 
momentum change is used for the blade design than is needed to 
satisfy the required power output. This correction would be suf- 
ficient for obtaining a blade profile that turns the gas an adequate 
amount. The gas now would enter the blades without any inci- 
dence and the blade load would increase gradually as the gas is 
turned by the blade. Cascade tests have shown, however, that 
any given blade chord can be utilized most effectively if the blade 
loadings are increased fairly rapidly near the inlet by means of a 
certain incidence angle. This practice results in more uniform 
blade loadings across the whole blade chord. Fig. 3 shows a de- 
sirable blade-load distribution for an axial-flow-type blade which 
can be obtained by proper selection of basic profile and angle of 
incidence. Thus the ideal velocity diagram is subjected to a 
second modification which takes into account a certain incidence 
angle at the inlet. Figs. 4 and 5 show modified velocity diagrams 
for nozzles and rotor blades, respectively. 

To determine the airfoil shapes of the blades, the general dif- 
ferential equations for the three-dimensional compressible flow 
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Fig. 3 Typical blade-load distribution for axial-flow turbine blades 


Fig. 4 Nozzle-velocity diagram modified for incidence and overturn 


Fig.5 Rotor-blade velocity diagram modified for incidence and overturn 
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are used. The derivation of these equations is shown in Appendix 
II. For convenience these equations are repeated here: 


For axial-flow turbines the radial velocity V, may be neglected 
and the equations can be simplified as follows: 


(1) 


(3) 


Equation (1) expresses radial equilibrium and may be used to 
calculate the fluid density at any radius station. Equation (2) is 
used to calculate the blade shape as explained later, while equa- 
tion (3) in conjunction with (1) and (2) provides a means for 
checking the density changes of the gas as it passes through the 
blades. 

Flow rate, gas densities, and velocity pattern determine the 
required flow areas. Based on the peripheral blade velocities and 
the flow areas, a first approximation of flow-annulus dimensions 
and wheel configuration can be made as shown in Fig. 6. The 
blade-chord length is selected so that satisfactory wheel propor- 
tions are obtained and the blade spacing is selected such that 
pitch-chord ratios below 1.0 are obtained at the outer diameter of 
the flow annulus. The main feature of the present blade-design 
method now consists of selecting a suitable airfoil and shaping it 
into a blade profile such that momentum requirements are satis- 
fied, together with requirements of optimum blade-load distribu- 
tion, flow rate, and stress. At first sight, this looks like an almost 
endless task based on trial and error. However, after the method 
has been applied a few times and the designer has developed a 
certain feeling for proper blade proportions, the method proves to 
be quite straightforward, and trial-and-error solutions can be 
avoided almost completely. 

It has been found most convenient to start by selecting a basic 
airfoil suitable for operation at high subsonic velocities. This air- 
foil is then shaped into a blade profile such that the tangents to 
the concave face of the blade satisfy the modified velocity diagram. 
Fig. 7 shows a typical rotor-blade profile that is obtained with 
this approximation method. Then the preliminary wheel con- 
figuration should be used for a careful stress study to avoid the 
need for major changes later, after the blade profiles have been 
calculated in detail 

The concave side of the blade profile is now checked according 
to equation (2), which can be rewritten as follows: 


_ _ do _ _ de 
prod dt dt dz dt 


dV, 
The pressure gradient —dp/(rd@) between two adjacent blade 
surfaces may be assumed to be constant. Therefore, it can be re- 
placed by a constant g/c, where g represents the blade load or the 


pressure difference between the two blade surfaces, and ¢ is the 
distance between the two blade surfaces. Solving for dV, we find 
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(4) 


This is the differential velocity increase in circumferential direc- 
tion over a blade length increment dz measured in axial direction 


for a given blade load q. 
The blade load g must now be determined from the momentum 


changes as shown in the velocity diagrams Figs. 4 and 5. The 
total momentum force on the blade area bh, Figs. 6 and 7, can be 
expressed as 
_ 

9 


Q 


where f,, is the mean flow rate for one blade passage of height 
h. The velocity change B may be obtained from Fig. 4 
for the nozzles and from Fig. 5 for the rotor. 

The mean flow rate f,, can be calculated from the inlet and exit 
conditions of the blade. Due to centrifugal forces the flow lines 
are crowded toward the outside in regions of high circumferential 
absolute gas velocities. When the gas flows in purely axial direc- 
tion, the flow lines are uniformly spaced. Because of this slight 
radial displacement of the streamlines which is caused by changes 
in flow direction, the flow rate f varies as the gas flows through the 
blades. For instance, in a nozzle diaphragm which imparts to the 
gas a considerable amount of whirl, f increases from inlet to exit 
at the outer shroud, it decreases from inlet to exit at the inner 
shroud, and remains constant at some intermediate radius station. 
In the turbine wheel these conditions are reversed because the 
whirl of the gas stream is reduced to obtain usable work at the 
rotor shaft. Since the change of f is small from blade inlet to 


AXIS OF ROTATION 
Fig. 6 Typical blading configuration of axial turbine stage 


Fig. 7 Typical rotor-blade profiles developed into a flat pattern 
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blade exit, a linear function may be assumed to establish inter- 
mediate values and a mean value of f for any given radius station 
can be calculated as follows: 


fiths 
2 

fi = 

fa = 


The mean blade load for a given radius station can be expressed 
as 


Im 


where 


and 


(5) 


With the mean blade load fixed, the actual design blade-load dis- 
tribution ean now be established as desired, for example, according 
to Fig. 3 


= InP (6) 


such that the total area under the blade load curve is equal to 


md. 
Combining equations (4), (5), and (6), we find 


The tangential velocity V, at any point along the blade chord 
is obtained by integration: 


7 


The initial velocity component Vy is given in the velocity dia- 
grams Figs. 4 and 5 for nozzle diaphragm and rotor as + Ay and 
— Ap, respectively. 

The displacement of the flow particles in the tangential direc- 
tion within a time element dt is 


V, = f/(¥ch) 


Equations (7) and (8) can be integrated with sufficient accuracy 
in finite steps if the blade chord b is divided into at least ten equal 
increments. Values for the variable c are obtained from Fig. 7; 
values for Y are established by assuming a linear drop in density 
from blade inlet to blade exit if the corresponding pressure 
changes are small. If the pressure drop from blade inlet to blade 
exit is large, then it is safer to assume a linear temperature drop 
to calculate the corresponding pressure and density decrease. 
Integration of equation (8) produces the concave side of the 
blade in accordance with the corresponding velocity diagram at 
any particular radius station. The final blade profile can now be 
completed by adding the blade thicknesses d as well as leading- 
edge radius and trailing-edge radius as assumed in Fig. 7. If the 
integration of equations (7) and (8) is carried out with sufficient 
accuracy, then the tangents of the concave side of the blade at 
blade inlet and blade exit will agree with the corresponding angles 
of the modified velocity diagrams Fig. 4 or Fig. 5 within a frac- 
tion of 1 deg. Also, the value of V, integrated over the whole 
blade chord, must check with the vector By, Fig. 4, or Bp, Fig. 5. 
To obtain a well-defined blade at least three sections per inch 
of radial height should be calculated. The station closest to the 


where 


and 
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hub should be located above the blade-root fillet so that the 
blade section can be checked dimensionally. The blade sections, 
calculated for the various radius stations, may be stacked up 
along any convenient radial line as long as a continuous blade 
shape can be obtained by smooth fairing between the calculated 
sections. For nozzle blades, it is usually advantageous to select 
as stacking points the centers of all trailing-edge radii. This way, 
a nozzle blade is obtained with a straight radial trailing edge. 
For rotor blades, it is most common to select as stacking points 
the centers of gravity of all sections. With this type of blade 
construction, blade stresses and danger of vibration are minimized 
and bending stresses due to centrifugal forces are completely 
eliminated. 

After the general blade shape has been essentially completed, a 
final stress check is made. If some minor changes in blade thick- 
ness are required to satisfy stress requirements, the changes 
should be made on the convex side of the blade. As a final step 
it may then be necessary to make some adjustments in over-all 
radial blade height in order to satisfy flow-area requirements. 

It is obvious from the foregoing discussion that a close initial 
approximation of the final blade shape is essential for successful 
application of the design method. It was found that even a 
relatively inexperienced designer was able to obtain a satisfactory 
blade according to this method after not more than three at- 
tempts. An experienced designer usually obtains the best possible 
blading arrangement in the first attempt with only minor adjust- 
ments being required between the initially assumed and the 


final blade shape. 


Application of Method to a Specific Case 


The blade-design method just described has been used at the 
authors’ company for many years to design the blading of axial- 
flow-type turbines. The most recent application was the power 
output stage of a 400-hp, two-shaft, gas-turbine engine. To 
present a clearer concept of what can be accomplished with the 
blade-design method, a more detailed account is given of the de- 
sign features of this turbine stage. 

An ideal free-vortex-flow pattern was selected for the blading 
The nozzles were designed for constant pitch-chord ratio at all 
radius stations. These assumptions resulted in a blade with very 
little twist and with very similar blade sections from root to tip. 
Hence it was only necessary to calculate the exact blade contour 
for the mean radius station. The blade profiles at root and tip 
stations were obtained by changing all principal dimensions of the 
calculated blade contour, by the ratio of the respective station 
radii. The three principal sections were stacked on a radial line 
connecting the centers of the trailing-edge radii in order to obtain 
a straight radial trailing edge. The positions of the sections with 
respect to the turbine axis were determined from velocity dia- 
grams suitably corrected to obtain the desired overturn angles. 

Fig. 8 shows a scale view of the three principal sections stacked 
up as described. The discharge flow area between the nozzles 
was made about 5 per cent larger than required ideally to account, 
for boundary-layer effects and to assure that the basic design as- 
sumption of an unchoked nozzle diaphragm would be achieved. 

The rotor blades needed considerably more cross-sectional area 
change from root to tip to satisfy aerodynamic as well as struc- 
tural requirements. Safe stresses demanded a considerable area 
taper from root to tip which resulted in higher pitch-chord ratios 
and higher overturn angles for the tip station than for the root 
station. Aerodynamic requirements called for a blade with a 
considerable amount of twist. With such drastic profile changes 
it became necessary to calculate the exact blade contour for three 
radius stations. Blade contours for two additional radius stations 
were obtained by interpolation. All five blade sections were 
stacked on a radial line connecting the centers of gravity of the 
sections. Again the sections were positioned with respect to the 
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Fig. 8 Principal sections of nozzle for a 400-hp turbine stage 
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Fig. 9 Principal sections of a rotor blade for a 400-hp turbine stage 


TURBINE AXIS 


axis of rotation according to velocity diagrams suitably corrected 
for angles of incidence and overturn. 

Fig. 9 shows a scale view of all five principal sections stacked 
in their relative positions. The effective discharge flow area in 
this case was estimated to be 96 per cent of the physical area. 
It should be emphasized here that an accurate prediction of the 
effective flow areas of nozzle diaphragm and rotor is indispensable 
if the desired degree of reaction in the turbine stage is to be ob- 
tained. A high degree of reaction at the tip was desired in this 
case requiring near sonic gas velocities relative to the blade, yet 
no choking effect at the blade exit could be tolerated. The physi- 
cal flow area at the blade inlet proved to be about 10 per cent 
larger than the required effective area. This may appear to be 
somewhat too large, but considering that the inflow conditions to 
the rotor blade are somewhat undetermined because of the wakes 
from the nozzles, an ample inlet flow area was considered de- 
sirable. Moreover, slightly too much area at the blade inlet is 
not expected to impair the turbine stage efficiency. 

The stage efficiency was predicted to be at least 80 per cent on a 
total-to-static basis and 90 per cent on a total-to-total basis. 
Enough energy was left in the exhaust gases leaving the stage to 
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raise the total-to-static efficiency to 84 per cent by means of a 
properly designed exhaust diffuser. 

Blade sections and reference co-ordinate systems, established 
by the aerodynamic designer, were transferred directly to produc- 
tion drawings. The engine designer only had to concern himself 
with actual manufacturing problems, such as shrink allowances 
for cast blades, tooling techniques, and methods of attaching 
the blades to rotor hub or nozzle shrouds. All fairing between 
the principal stations shown on the production drawings was left 
up to the pattern maker without further help from Engineering. 
Very satisfactory blade shapes were obtained with this procedure. 


Test Results 


A turbine stage according to the design just described was 
fabricated for performance evaluation in a test rig. Fig. 10 shows 
the assembled turbine wheel. 

The turbine was tested over as wide an operating range as 
available test facilities permitted. Expansion ratios and velocity 
ratios ranged from 1.6 to 2.7 and from 0.40 to 0.80, respectively. 
The design point with an expansion ratio of 2.10 and a velocity 
ratio of 0.57 lay just about in the center of the whole test range. 
No exhaust diffuser was used for the test. Figs. 11 and 12 are 
performance maps calculated from the test data. U,, is the blade 
speed at the mean radius station between hub and tip. 

Fig. 11 shows that the predicted total-to-static stage efficiency 
of 80 per cent at the design point was exceeded by about 3 per 
cent. This fact indicates that an aerodynamically very clean 
blade configuration had been obtained which provided 
higher flow coefficients in the rotor passages than expected. The 
larger than predicted effective flow area resulted in lower dis- 
charge velocities and correspondingly higher total-to-static stage 
efficiencies. It still would be possible to raise the total-to-static 
efficiency an additional 2 or 3 per cent by using a properly designed 
exhaust diffuser. The possible gain that could be expected must 
be weighed carefully, however, against the additional size and 
weight of the required duct work. 

The predicted total-to-total efficiency of 90 per cent was 
achieved or even slightly exceeded as shown in Fig. 12. Again 
there is a clear indication that the effective flow area of the rotor 
was larger than expected. This resulted in more pressure drop 
across the nozzles and less pressure drop across the rotor com- 
pared to design conditions. The degree of reaction thus was lower 
than intended, as borne out by a slight shift of the optimum ve- 
locity ratios toward lower than design values. To establish a 
reference for typical values of optimum velocity ratios, it is helpful 
to remember that for a pure impulse stage, having a total-to- 
static efficiency of 81 per cent, the optimum velocity ratio is 0.45 
by definition. The corresponding value for a 50 per cent reaction 
stage is 0.63. 

Theoretically, then, the throat area of the rotor should have 
been reduced a certain amount to obtain the intended amount of 
reaction and therefore highest possible efficiency. In practice, 


Fig. 10 Assembled rotor for a 400-hp turbine stage 
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Fig. 11 Test performance of 400-hp turbine stage. Total-to-static 
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Fig. 12 Test performance of 400-hp turbine stage. Total-to-total 


however, the obvious remedy was to run the wheel slightly slower 
than planned because only a minor improvement could have been 
expected from a modified wheel. A further justification for ac- 
cepting the wheel without modifications was found in the in- 
tended application of the turbine as a power-output stage of a 
free-turbine engine. In order to obtain the desired gas-producer 
turbine-inlet temperature in such an engine, it is often necessary 
to make certain final flow-area adjustments in the nozzles because 
of unavoidable manufacturing discrepancies. A slight increase in 
the nozzle flow area of the power turbine will reduce the inlet 
temperature of the gas-producer turbine. At the same time, such 
a flow-area change would have the tendency to restore the cor- 
rect degree of reaction in the power turbine. 

Figs. 11 and 12 indicate that some very high efficiencies were 
owtained at lower than design expansion ratios. It must be men- 
tioned, however, that these favorable results are not beyond 
question, because at low expansion ratios the test readings become 
relatively less accurate. Of more significance is the fact that the 
efficiencies do not decrease at expansion ratios considerably higher 
than design. This means that the oblique nozzles of the dia- 
phragm assembly can produce higher than sonic velocities ef- 
ficiently in the vaneless space between the stator and rotor 
blades. Also, the rotor, which in effect is a rotating nozzle, seems 
to be capable of producing transonic relative exit velocities ef- 
ficiently. 

Transonic velocities are obtained behind any set of oblique 
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RAD U Vi | V2 | Wy | We | ot: | B, B, 
769 | 1480 | 505 | 791 | 912 | 18-1°| 90.7°] 36.1°| 33.6° 
875 | 1318 | 467 | 580 | 1007 | 20.3°| 88.6°| 52.5°| 27.6° 
‘3 | 982 | 1150 | 448 | 452 [1038 | 22.7°| 96° | 81.4°| 25.2° 


Fig. 13 Velocity diagrams for radius stations r), r2, and r;, calculated 
from test data. 


RAD U vi V2 | Wi | | | Bi | B2 
1344 | 517] 655 | 929: 18.1°| 90° | 39.5°/ 341 


Fig. 14 Design-velocity diagrams for radius stations 1, r2, and r3 


annular nozzles, designed for subsonic flow if a high enough ex- 
pansion head is available. The velocity at the nozzle throat then 
becomes sonic and further expansion takes place downstream be- 
cause more flow area becomes available as the gas leaves the 
throat. Overexpansion can be utilized to good advantage in the 
design of axial-turbine stages with very high expansion ratios. 

The performance maps further show that the stage efficiencies 
are more sensitive to underspeeding than to overspeeding. This 
effect also can be explained quite logically by considering the 
varying flow conditions, set up at off-design speeds by the fixed 
flow-passage configuration. However, a thorough discussion of 
off-design operating conditions is beyond the scope of this paper. 
Instead, it is considered more appropriate here to discuss briefly 
the flow pattern that was obtained at the test point closest to the 
design point. 

To provide a clear idea of the flow pattern, the following infor- 
mation was recorded: Static pressure, total pressure, and total 
temperature ahead of the nozzles; static pressures at inner and 
outer shroud at nozzle exit; total pressure and flow direction at 
three different radius stations about 2 in. downstream from the 
rotor blades; static pressure and total temperature in the dis- 
charge chamber; air flow rates by means of a standard flow nozzle 
in the discharge duct. No attempt was made to measure total 
pressure and flow direction at the nozzle exit because previous 
experience had shown that such readings, to be reliable, must be 
taken at least 2 in. downstream from the blades. Even the static 
pressure readings at the nozzle exit are sometimes difficult to in- 
terpret because of prevailing transonic-flow conditions. With no 
flow-direction readings at the nozzle exit available, the assump- 
tion was made that the gas was leaving the nozzles exactly ac- 
cording to design calculations. This assumption was justified 
because the nozzle overturn angles had been established very 
carefully in earlier cascade tests. 
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The assumed nozzle-exit flow directions, together with the 
measured test data, permitted construction of velocity diagrams 
for the three radial stations where turbine-exit flow-directions 
readings were available. Fig. 13 shows these calculated velocity 
diagrams. For comparison, the design velocity diagrams are 
shown in Fig. 14 for the same three radius stations. Evidently 
a very close approximation of the design flow pattern was ob- 
tained in the test. The absolute discharge velocities were nearly 
equal in magnitude and parallel to the axis of rotation as planned. 
The incidence angles at the rotor inlet were somewhat larger than 
intended; about 3 deg larger near the root and 6 deg larger near 
the tip. No loss in efficiency was experienced, however, because 
the relative rotor-inlet velocities were low, especially near the tip 
where the incidence angle was a maximum. Well-rounded rotor 
leading edges also helped, preventing a serious loss. Significantly, 
the rotor turned the air more than expected; about !/, deg more 
near the root and 2!/, deg more near the tip. These test results 
furnish evidence that an increase in rotor efficiency not only pro- 
vides an increase in effective rotor flow area but also reduces the 
required amount of blade overturn. 

Close examination of Fig. 13 reveals that the work output at 
the mean station is about 5 per cent and 10 per cent higher than 
at the hub station and tip station, respectively. Theoretically, 
the work output should have been equal at all radius stations be- 
cause the blading design was based on an ideal free-vortex flow 
pattern, It cannot be expected, of course, that in any actual 
test perfect free-vortex velocity diagrams will be obtained be- 
cause the accuracy of any test data is never perfect, and the flow 
pattern of a compressible gas always shows some unpredictable 
irregularities. 


Review and Conclusions 


A blade-design method has been developed that takes into 
consideration the basic theoretical concepts of fluid dynamics as 
well as experimentally established design parameters. The 
method can be applied by the average turbine designer after a 
short period of special training. No special computing equip- 
ment is required to establish the blade profiles unless a compre- 
hensive study is undertaken to show how the blade profiles change 
for given variations of aerodynamic blade loadings. The integra- 
tions necessary to calculate the blade profiles can be carried out 
with sufficient accuracy in finite steps according to conventional 
methods if the blade chord is divided into at least ten equal in- 
crements. By applying the method to a specific case, it is ex- 
plained how well-proportioned blades can be obtained that 
satisfy aerodynamic as well as structural requirements. The 
basic blade profile is fully determined from aerodynamic con- 
siderations, while experience factors are utilized only for refine- 
ments. The blade shapes are defined in a co-ordinate system that 
lends itself readily to manufacturing purposes. Test data, ob- 
tained from a turbine, designed according to this method, show 
that a very close approximation of the design flow pattern can be 
obtained. Specifically, it was demonstrated that nearly ideal free- 
vortex flow can be maintained across the blading with suitably 
designed blade profiles. Aerodynamic performance of the turbine 
stage exceeded expectations. 

Success of any specific blade design, of course, still depends to 
a very large extent on the good judgment of the designer, but the 
described method provides him with a powerful design tool which 
virtually eliminates all guesswork. 


APPENDIX 


Definition of Performance Parameters 
for Axial-Flow-Type Turbine Stage 


Generally used performance parameters for an axial-flow-tur- 
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Fig. 15 Typical enthalpy-entropy diagram for axial-flow turbine stage 


bine stage are most conveniently defined by referring to a typical 
enthalpy-entropy diagram and a corresponding velocity diagram 
as shown in Fig. 15. Each velocity, indicated in the velocity 
diagram, has its corresponding enthalpy head in the entropy dia~- 
gram. The dimensionless performance parameters—nozzle ef- 
ficiency, rotor efficiency, over-all stage efficiency, degree of reac- 
tion, and velocity ratio—are defined as follows: 


hy/ha nozzle efficiency, total-to-static 
Au/hs rotor efficiency, total-to-static 
rotor efficiency, total-to-total 
he/ho = stage efficiency, total-to-static 
he/hio = stage efficiency, total-to-total 
hs/(he + he) = degree of reaction 
U/Vo = velocity ratio 
Vo? 


——,, total available expansion head 


where ho 


2 
hi nozzle-exit, absolute-velocity head 


297 


V2? 
—, rotor-exit, absolute-velocity head 


Wi? 


——,, rotor-inlet, relative-velocity head 


2gJ 


2 
= rotor-exit, relative-velocity head 


2_ Ww 
stage-reaction head 


exp 
Symbols for velocities and station designations, indicated by 
subscripts, are defined in Fig. 15. 


APPENDIX Il 


Derivation of Differential Equations for 
Three-Dimensional Compressible Flow 
For the aerodynamic design of turbomachinery, it is most con- 
yenient to set up the basic flow equations in terms of cylindrical 
co-ordinates. To simplify the problem, the case of a mass element 
moving at a constant radial velocity V, along a radial beam, which 
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hy 


1 
de Yu=rw 
| Vu=(r+ dr) w 
4 =Vu+wdr 
p 


Fig. 16 Velocities and displacements for uniform differential motion of 
a mass element from point 1 to point 2 


rotates at constant angular velocity w around a fixed center, is 
considered first. Acceleration forces created by these uniform 
motions are calculated. In the second step, forces caused by ac- 
celerations in radial, tangential, and axial direction are superim- 
posed to make the equations applicable to any type of three-di- 
mensional flow pattern. 

Fig. 16 shows the velocities and displacements to be considered 
if a mass element moves from point 1 to point 2 at constant V, 
and constant w. 

The velocity increments in the tangential direction which are 
caused by the velocity V, and by the radial displacement dr, 
dV,, and dVy,, respectively, can be calculated as follows: 


dV,, = V, sin do 
= — V, cos dd 
= V, + wdr — V, cos dd 
= wdr + V,(1 — cos dd) 
The velocity increments in the radial direction caused by V, 
and dr, dV,,, and dV,,, are 


V,’ — V, = V,{1 — cos dg) 
= —V,sindd 


For an infinitely small d@¢ these expressions can be simplified as 
follows: 
wdr 
dV,, 0 
—V,do 


The acceleration in the tangential direction can now be written 


dt 


a, = Vw+ovV, = 


This is the Cariolis acceleration. 
The acceleration in radial direction is 
dVu, 
dt 


(10) 


a, 


This is the centripetal acceleration. 

If, in addition to the accelerations just described, the mass ele- 
ment has linear accelerations ag in the direction of the beam and 
a, in the direction of the axis while the beam itself has an angulag 
acceleration @ in the tangential direction, then the resulting ac- 
celerations in the three principal directions become 
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Fig. 17 Mass element and differential forces in the three principal direc- 
tions 


dV, 


d 
a, = ar + WV, = + Ww, 


a, = | 


dt 


With the accelerations known, the forces dP acting on a mass 
element dm can be expressed according to Newton’s first law: 


dP, = dma, ) 
dP,, = dma, 
dP, = dma, 


(12) 


where 


dm = rd@¢ dr dz p (13) 


Fig. 17 shows a mass element in cylindrical co-ordinates with 
differential forces acting in the three principal directions. The 
forces dP also may be expressed in terms of pressure gradients: 


dP, = dz 


2p 
dP, = ae rd@ dr dz 


op 
dP, = de dr 


By substituting expressions (11) and (13) in (12) and equating to 
(14), we find: 


dV, 
rat) = dr dz 


d 
r d@ dr dzp (2 r+ rd@ dr da 


dV. 
dr dzp = dz dr 


Cancelling common factors, we find the general differential equa- 
tions for three-dimensional compressible flow: 


(15) 
(16) 
(17) 
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The Basic Heat Transfer and Flow 
Friction Characteristics of Six 
Compact High-Performance 

Heat Transfer Surfaces 


In this paper the basic characteristics of six new compact high-performance surfaces 
are presented. Four are of the strip-fin type, including three stacked in a unique mul- 
tiple sandwich arrangement. One of the surfaces has perforated fins formed by punch- 
ing circular holes through the fins, and one has wavy fins. The basic characteristics 
are presented in terms of NgNp,/* and f versus Reynolds number. A suffi- 
cient geometrical description is provided so that the data are readily usable in heat ex- 
changer design analysis. The triple-sandwich surface introduces some new complica- 
tions in the analysis of the effectiveness of the fins, and this problem is considered in an 


Appendix. 


Introduction and Objectives 


I. 1947 A RESEARCH PROGRAM was established at 
Stanford University, under the sponsorship of the Office of Naval 
Research, the Bureau of Ships, and the Bureau of Aeronautics, 
for the purpose of investigating the basis heat transfer and flow 
friction characteristics of compact high-performance heat transfer 
surfaces. The original impetus for this work came from the ad- 
vent of the gas turbine as a promising transportation-type prime 
mover, and the further promise offered by inclusion of heat ex- 
changers in the cycle. The shell-and-tube heat exchanger, while 
perhaps satisfactory for a stationary power plant, looms very 


Contributed by the Gas Turbine Power Division and presented 
at the Gas Turbine Power Conference and Exhibit, Cincinnati, Ohio, 
March 8-11, 1959, of THe American SocteTy oF MECHANICAL EN- 
GINEERS. 


large and bulky for the transportation application, and it was felt 
worth while to investigate some of the much more compact heat 
exchanger surfaces, especially of the extended-surface variety, 
that have been used extensively in aircraft intercoolers, air- 
conditioning systems, and many other primarily low temperature 
applications. 

The objectives of the program were to first determine by direct 
test the basic characteristics of a large variety of different sur- 
faces and, secondly, to investigate the flow and heat transfer 
mechanisms so as to form a rational basis for the design of new 
surfaces of superior characteristics. Most of the surfaces tested 
were of low temperature construction, usually aluminum, but the 


Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, December 
31, 1958. Paper No. 59—GTP-2. 


Nomenclature 
English Letter Symbols 


A = total heat transfer surface area, ft? 
= free-flow area on one side of a heat 
exchanger core, ft* 
Ay, = surface area of the fins, ft? 
= a fin element length, ft, see Fig. 8 U 
= plate spacing, ft; also a fin element 
length, ft, see Fig. 8 V 
specific heat at constant pressure, w 
Btu/(Iby deg F) 
total drag force, lbr 
core mass velocity based on A,, 
/(hr ft?) 
proportionality factor in Newton’s 
Second Law, 32.2 (Ibuft)/(Ibr- 
sec?) 
convection heat transfer conduc- 
tance, Btu/(hr ft? deg F) 
thermal conductivity, Btu/(hr ft? 8 
deg F/ft) 
total flow length through heat ex- 


A.L/A 


Greek Letter Symbols 


Fig. 8 


Fig. 8 
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changer core, ft; also a fin ele- 6 
ment length, ft, see Fig. 8 
M = afin element length, ft, see Fig. 8 7 
hydraulic radius, ft, defined as 


= over-all unit conductance, based on 
A, Btu/(hr ft? deg F) 

= velocity, ft/sec 

= flow rate, lb/hr 


a = thickness of a fin element, ft, see 


8 = surface area density, surface area 
per unit of volume between the 
plates of one side of the heat 
exchanger, ft?/ft?; also thick- 
ness of a fin element, ft, see Fig. 


y = thickness of a fin element, ft, see 


thickness of a fin element, ft, see 
Fig. 8 

viscosity coefficient, lby/(hr ft) 

Ny fin effectiveness, dimensionless 

n. = over-all surface effectiveness, di- 
mensionless, see Equation (2) 

p = fluid density, Iby/ft® 

T. = equivalent shear stress at the 

surface, F/A, lbr/ft? 


Nondimensional Groupings 
NgNpr/* = convective heat transfer 
parameter, (h/Gc,) 
(uc,/k)*/*, Stanton num- 
ber times the two-thirds 
power of the Prandtl 
number 
Nr = Reynolds number, 47,G/u 
f = friction factor, To/(pV2/2g,) 
NTU = number of transfer units, 
(AU)/(we,) 
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surface characteristics have nothing to do with the material of 
construction, and it was felt that if such surfaces showed suf- 
ficient promise for the gas turbine application there would be in- 
centive to develop methods of construction using high-tempera- 
ture materials. 

As the program progressed, and the performance data became 
available in a form useful to heat exchanger designers, and to 
those making feasibility studies for various applications, it be- 
came apparent that such data would be useful in a far wider 
variety of applications than merely in the gas turbine. Thus the 
program was broadened to include surfaces that might not be of 
direct interest to the gas turbine application. Subsequently the 
AEC joined in sponsorship of the program. 

The test results and the various conclusions that have been 
reached relative to the basic mechanisms and the ultimate possi- 
bilities of compact high-performance surfaces have been published 
in a number of papers in the ASME Transactions. In 1954 the 
characteristics of 88 surface configurations were summarized in a 
report that has subsequently been made available in book form 
{1].!. This summary includes not only all of the data obtained as 
a part of this project, but also a substantial quantity of similar 
data obtained elsewhere. 

Since 1954 a number of new and interesting surfaces have been 
introduced by various manufacturers, and it was decided to re- 
institute the test program on a modest scale. To date six addi- 
tional surfaces have been tested, and it is the objective of this 
paper to present these new data so as to make them generally 
available for heat exchanger design and feasibility studies, and 
to supplement and augment conclusions previously reached. 


Experimental Apparatus and Reduction of Data 


Each of the test surfaces are built into steam-to-air test cores, 
8*/; X 9%/, in. frontal dimensions, with the surface under con- 
sideration on the air side. The test apparatus consists of an in- 
duced-draft air duct in which air flow rate can be accurately con- 
trolled and metered, and in which air temperature can be ac- 
curately determined. Slightly superheated steam is supplied to 
the steam-side of the test cores and condenses on the steam-side 
surface. A considerable excess of steam is passed through the 
core to prevent the build-up of a thick condensate layer. Instru- 
mentation is provided to measure air, condensate, and blow-steam 
rates, air entering and leaving state, and steam entering and lJeav- 
ing state. A complete and detailed description of this apparatus 
is given in reference [2] and thus will not be elaborated upon here. 
It will suffice to say that provision is made to establish an energy 
balance between the steam and air sides of the test cores, and the 
energy unbalance for the test cores considered here was consist- 
ently less than 3 per cent. An estimate of the experimental un- 
certainty is also given in reference [2], with uncertainty intervals 
established at +5 per cent for the friction factor f and the heat 
grouping NaNp,’*. The uncertainty interval in the Reynolds 
number is +2 per cent. 

The details of the method of reduction of data are presented in 
reference [2]. It should be noted here that the steam-side heat 
transfer resistance and the ineffectiveness of the extended surfaces 
(fins) are taken into account so that the reported results pertain 
to the surface geometry alone and are independent of the steam- 
side and the material of the fins. 


Method of Presentation of the Results 


The heat transfer and flow friction characteristics, respectively, 
are presented using the dimensionless correlations, Ng:Np,'/* 
versus Np and f (Fanning friction factor) versus Nx. These re- 
sults are given in tabular form in the Tables and also in graphical 


' Numbers in brackets designate References at end of paper. 
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Tables Summary of basic heat transfer and friction data (h/Gc,)Npr’/* 
f versus Nr from smoothed curves 


Surface 1/4-15.4(D) 


(n/Ge, 


00510 
200537 
00570 
-00617 
-00650 


00692 
00756 
00809 
00864 
00952 


-01107 
01227 
-01407 
*,0169 


Surface 1/8-16.12(D) 


00721 
-00764 
.00622 
00859 
-00908 


-00987 
-01060 
201123 
201205 
201352 
20150 
-0176 
-0226 


Surface 1/8-16.12(T) 


(n/ae, 


.00705 
-00773 
-00820 
00880 


-00967 
-01045 
201122 
-01225 
01395 
20155 


0219 


Surface 13,95(P) 


00503 
00547 
00600 
00631 
-00660 
00695 
.00708 
00700 
200736 
00817 
00893 
01010 
-01180 
20132 


form, Figs. 1 to 6. Two sets of tests are run for each test core. 
One is called a hot-core test, in which heat transfer takes place 
and both heat transfer and friction measurements are made; 
the other is an isothermal test in which only friction measurements 
are made. On Figs. 1 to 6 the curves representing the best in- 
terpretation of the friction factor data are drawn through the 
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ak 
tA ‘ 
5000 .0167 
5000 
20202 
.0213 
1500 .0228 
1200 -0255 
1000 3 
000 
600 .0362 
400 "0507 
700 0601 
000 
i000 .0310 
3000 0 315 
2500 2033 
2000 10357 
1500 0379 
1200 -0400 
"B00 
800 
600 
00 
300 0781 
-0937 
5000 
20358 
3000 
1200 
.0496 
600 -0552 
.0721 
400 
.0830 
.1003 
= 
.01168 
8000 
6000 01285 
.0146 
3000 .0155 
.0169 
1500 .0178 
“B00 .0232 
800 
600 


Tables (continued) 
Surface 1/8-13.95 
2/3 
Ne (h/Ge, ) f 

8000 

6000 .0110 

5000 -01170 0664 

4000 .01250 0684 

3000 .0712 

2500 20144 00733 

2000 20155 20765 

1500 0817 

1200 .0181 

1000 .0192 092 
800 1020 
600 .0223 -1170 
500 -0233 
400 0247 2154 

Surface 11.5=-3/8(Ww) 

10000 00686 0331 
8000 00746 
6000 .00831 
5000 -00890 0427 
4000 00970 
3000 ~01077 20525 
2500 201155 -0567 
2000 
1500 
1200 0150 0779 
1000 0158 0845 

800 
600 .0178 21035 
500 .0185 e111 
400 


test points for the isothermal flow tests, since for these there is no 
question of the effects, if any, of the temperature varying fluid 
properties, The curves representing the best interpretation of 
the heat transfer data lie somewhat above the actual test points 
due to an allowance for an estimate of the steam-side heat transfer 
resistance. This method of presentation has been arbitrarily 
adopted simply because the uncertainty in the steam-side estimate 
is much greater than for any of the other quantities, and it is 
possible to see at a glance the effect of this uncertainty. 

The Reynolds number Ng = (4r,G/u) is evaluated on the 
basis of the flow passage hydraulic radius and a viscosity evaluated 
at the logarithmic mean temperature, which is essentially the 
mean of the bulk air temperature with respect to flow length. 
This temperature is also used to evaluate Prandtl number in the 
heat transfer correlation. The hydraulic radius is defined as 
follows: 


= A,/A 
r, = A,L/A 


For cylindrical tubes of any cross-section geometry, A,L/A is 
identical with the flow cross section divided by the wetted 
perimeter, which is the conventional definition of hydraulic radius. 

The friction factor is evaluated on the basis of an equivalent 
shear stress on the heat transfer surface area A which is in turn 
simply defined as the total drag forces opposing the flow divided 
by the heat transfer area. Thus, 


t = F/A = foV*/2g, 


This pseudo shear stress can be made up of both true shear and 
pressure forces acting on elements of surface not parallel to the 
flow direction. 

The tabular ‘Summary of basic heat transfer and friction 
data’’ is taken directly from the curves representing the best in- 
terpretation of the data. 
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Description of the Test Surfaces 


The six surfaces to be considered here are all of the plate-fin 
variety, i.e., fins are bonded between parallel plates to form con- 
tinuous flow passages. The plates can then be stacked sandwich- 
type to form crossflow or counterflow heat exchangers with each 
of the two fluids flowing between alternate pairs of plates. A 
virtue of the plate-fin arrangement lies in the great flexibility af- 
forded the designer who can employ entirely different fin arrange- 
ments and configurations, and different plate spacing, on the two 
fluid sizes of the heat exchanger. But more important, this ar- 
rangement lends itself to very great area compactness (heat trans- 
fer area per unit of volume), and it becomes a simple matter to 
substantially increase the heat transfer conductance by the use of 
various methods of boundary layer interruption and turbulence 
promotion. When high heat transfer conductances are obtained 
by such methods the surface is often termed a high-performance 
surface. High performance is generally bought at the expense 
of increased friction factors, but it can be readily demonstrated 
that in terms of over-all heat exchanger weight and volume, a 
very substantial increase in friction factor can be tolerated as the 
price for a modest increase in the heat transfer parameter NgNV. pr? 

Diagrams of the surfaces considered here are shown in the 
upper right-hand corner of Figs. 1 to 6. All six can be classed 
as high-performance surfaces in that they all incorporate methods 
of boundary layer interruption to increase the convective heat 
transfer conductance. The surfaces shown on Figs. 1, 2, 3, and 4 
will be referred to as strip-fin surfaces. The fins are formed from 
a continuous sheet of material that is folded to form flow passages, 
cut and offset at intervals in the flow direction to provide bound- 
ary layer interruption, and then is brazed between parallel plates. 
The surface shown on Fig. 4 is the more conventional type of 
strip-fin surface with a single set of fins between the base or prime 
surface. The characteristics of a number of surfaces of this type 
are presented in reference [1]. 

The surfaces shown on Figs. 1, 2, and 3 are unique in that they 
are multiple sandwich arrangements. Two and, in the case of 
Fig. 3, three sets of fins are brazed between a single pair of base 
plates. The sets of fins are separated by splitter plates to which 
they are also brazed. The splitter plates provide additional ex- 
tended or fin surface, and in the case of the triple sandwich, Fig. 3, 
the heat to the center passage fins must first be conducted through 
the splitter plates. The idea behind the multiple sandwich ar- 
rangement is the greater flexibility allowed the designer and 
manufacturer. In a gas-to-gas heat exchanger where there is a 
large disparity in gas density on the two sides of the heat ex- 
changer, a balanced design will usually call for more surface area 
and more flow area on the low density gas side. With the multiple 
sandwich arrangement, a triple sandwich might be used on the 
low density side while a single sandwich could be used on the high 
density side. The actual fin material on each side could be 
identical, which has obvious advantages to the manufacturer in 
terms of the number of dies which must be stocked and main- 
tained. The gas turbine regenerator is a classic example of a gas- 
to-gas heat exchanger with a large density difference on each 
side, and this scheme, which has so far been employed only for 
low temperature applications, is offered for consideration in the 
gas turbine application. 

The individual flow passages in a multiple sandwich surface 
are identical, so one would expect the same basic characteristics, 
Ng. p,'/* and f versus Np, regardless of whether the arrange- 
ment is of the single, double, or triple-sandwich type, provided the 
same fin material is used. However, the fin effectiveness analysis 
is more complex for the double and triple sandwich than for the 
single sandwich (see Appendix), and the triple sandwich appears to 
introduce other difficulties to be discussed later. The double and 
triple sandwich surfaces shown on Figs. 2 and 3 have identical 
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INTERPRETATION 


O03 04 0506 08 10 152.0 30 40 5060 80 


Fig. 1 Heat transfer and friction data, surface 1/4-15.4(D) 
Surface Data, Air Side: 


Fin pitch, 15.4 fins per inch 

Plate spacing, b = 0.206 in. 

Splitter symmetrically located 

Fin length (flow direction), 0.250 in. 

Flow passage hydraulic diameter, 4r, = 0.00527 ft 

Fin metal thickness, 0.006 in. 

Total heat transfer area/volume between plates, 8 = 642 ft?/ft® 
Fin area (including splitter)/total area, 0.816 


©, Hot core test points. Heat-transfer data evaluated on the basis of zero 


steam-side resistance. 
x, Cold core friction factors. 
Best interpretation of friction data based on cold core data. 


Best interpretation of heat-transfer data based on an allowance for 


steam-side resistance. 


| | 


0.010 2 
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Fig. 3. Heat transfer and friction data, surface 1/8-16.12(T) 

Surface Data, Air Side: 

Fin pitch, 16.12 fins per inch 

Plate spacing, b = 0.314 in. 

Splitters symmetrically located for three even passages 

Fin length (flow direction), 0.125 in. 

Flow passage hydraulic diameter, 4r;, = 0.00514 ft 

Fin metal thickness, 0.006 in. 

Total heat transfer area/volume between plates, 8 = 650 ft?/ft® 

Fin area (including splitters)/total area, 0.882 


©, Hot core test points. Heat-transfer data evaluated on the basis of zero 


steam-side resistance. 
x, Cold core friction factors. 
Best interpretation of friction data based on cold core data. 


Best interpretation of heat-transfer data based on an allowance for steam- 


side resistance. 
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Fig.2 Heat transfer and friction data, surface 1/8-16.12(D) 
Surface Data, Air Side: 
Fin pitch, 16.12 fins per inch 
Plate spacing, b = 0.206 in. 
Splitter symmetrically located 
Fin length (flow direction), 0.125 in. 
Flow passage hydraulic diameter, 4rn = 0.00509 ft 
Fin metal thickness, 0.006 in. 
Total heat transfer area/volume between plates, 8 = 660 ft*/ft® 
Fin area (including splitter)/total area, 0.823 
0, Hot core test points. Heat-transfer data evaluated on the basis of zero 
steam-side resistance. 
x, Cold core friction factors. 
Best interpretation of friction data based on cold core data. 
Best interpretation of heat-transfer data based on an allowance for steam- 
side resistance. 
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Fig. 4 Heat transfer and friction data, surface 1/8-13.95 
Surface Data, Air Side: 
Fin pitch, 13.95 fins per inch : 
Plate spacing, b = 0.375 in. 
Fin length (flow direction), 0.125 in. 
Flow passage hydraulic diameter, 4r, = 0.00879 ft 
Fin metal thickness, 0.010 in. 
Total heat transfer area/volume between plates, 8 = 381 
Fin area/total area, 0.840 
Note: The fin surface area on the leading and trailing edges of the fins 
has not been included in computations. 
0, Hot core test points. Heat-transfer data evaluated on the basis of zero 
steam-side resistance. 
x, Cold core friction factors. 
Best interpretation of friction data based on cold core data. 
Best interpretation of heat-transfer data based on an allowance for steam- 
side resistance. 
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Fig. 5 Heat transfer and friction daia, surface 13.95(P) 


Surface Data, Air Side: 

Fin pitch, 13.95 fins per inch 

Plate spacing, b = 0.200 in. 

Fin center material perforated with 0.079-in-diam holes spaced 32 per 
square inch in square pattern (16 per cent open). Position of holes in 
walls of various passages differs. 

Flow passage hydraulic diameter, 4r, = 0.00822 ft 

Fin metal thickness, 0.012 in 

Total heat transfer area/volume between plates, 8 = 381 ft*/ft® 

Fin area/total area, 0.705 

Note: Because of the very thick fins the surface area in the rim of the 
holes has been included in computations. 

©, Hot core test points. Heat-transfer data evaluated on the basis of zero 
steam-side resistance. 

x, Cold core friction factors. 

Best interpretation of friction data based on cold core data. 

Best interpretation of heat-transfer data based on an allowance for 
steam-side resistance. : 
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Fig.6 Heat transfer and friction data, surface 11.5-3/8(W) 


Surface Data, Air Side: 
Fin pitch, 11.5 fins per inch 
Plate spacing, b = 0.375 in. 
Flow passage hydraulic diameter, 4r, = 0.00993 ft 
Fin metal thickness, 0.010 in. 
Total heat transfer area/volume between plates, 8 = 347 ft?/ft? 
Fin area/total area, 0.822 
o, Hot core test points. Heat-transfer data evaluated on the basis of 
zero steam-side resistance. 
x, Cold core friction factors. 
Best interpretation of friction data based on cold core data. 
Best interpretation of heat-transfer data based on an allowance for 
steam-side resistance. 
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flow passages in all respects, the difference being only in the 
sandwich arrangement. The purpose in testing these two was to 
determine what, if any, new problems would be encountered in 
the triple sandwich arrangement. The double sandwich surfaces 
shown on Figs. 1 and 2 differ materially in that the uninterrupted 
fin-flow lengths are respectively 1/, and-1/; inch. 

The surface shown on Fig. 5 employs a novel system of boundary 
layer interruption. The fin material is perforated with holes be- 
fore being folded and brazed between the plates. The surface 
shown on Fig. 6 has fins bent in the flow direction to form wavy 
passages. Boundary layer separation apparently occurs at the 
turns to provide a high conductance, so this is also a high-per- 
formance surface. Other examples of this type of surface are pre- 
sented in reference [1]; in fact, this particular surface differs only . 
in a minor way from the surface designated 11.44—/,(W) in that | 
reference. 

A semidescriptive method of surface designation is employed, 
consistent with the designations used in reference [1]. For the 
strip-fin surfaces the first number refers to the fin-flow-length in 
the flow direction, e.g., 1/5 inch or 1/, inch. The second number 
refeis to the number of fins per inch transverse to the flow direc- 
tion. The letter in parenthesis refers to the double or triple 
sandwich construction. Thus the strip-fin surfaces shown on 
Figs. 1, 2, 3, and 4 are designated, respectively, 


1/,-15.4(D) 
1/s-16.12(D) 
1/_-16.12(T) 
1/_-13.95 


The perforated fin surface shown on Fig. 5 is designated by the 
number of fins per inch with the letter P referring to the per- 
forated fins. Thus this surface is designated, 


13.95(P) 


The wavy fin surface shown on Fig. 6 is designated by the number 
of fins per inch, the wave length, and the letter W to designate the 
wavy fin construction. Thus this surface is designated, 


11.5-*/,(W) 


The captions to Figs. 1 to 6 contain the complete dimensional 
data for each of the surfaces, including the hydraulic diameter and 


_ the surface area density. 


Discussion of the Test Results 


The Strip-Fin Surfaces. All of the strip-fin surfaces, Figs. 1 to 4, 
exhibit the general characteristics of such surfaces noted pre- 
viously. The heat transfer mechanism is primarily conduction 
through a laminar boundary layer that originates and grows on 
each individual fin. Turbulence appears to play a relatively minor 
part, although there is some evidence that the very high heat trans- 
fer performance of surface '/s-13.95 is in part attributable to a 
disturbance of the otherwise laminar boundary layer by turbulent 
eddies shed from the blunt leading and trailing edges of the rela- 


_ tively thick (0.010 inch) fins. All of these surfaces have rather 


high friction factors attributable to the form or pressure drag 
offered by the blunt leading and trailing edges of the fins. 

On Fig. 7 the test results for the three multiple sandwich sur- 
faces are compared. Comparing the two double sandwich sur- 
faces, it will be noted that cutting and separating the fins every 
1/, inch rather than every '/, inch (this is the only major dif- 
ference between the surfaces) results in a substantial increase in 
NsNp,’/* and a very large increase in f. The increase in heat 
transfer performance is a result of the thinner boundary layer. 
The increase in friction factor is attributable to both the thinner 
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Fig. 7 Comparison of the test results for the multiple sandwich surfaces 


boundary layer with its larger average shear stress, and to a 
greater degree the increased form drag resulting from there being 
twice as many fin edges. 

Surfaces 1/s-16.12(D) and !/s-16.12(T) are identical save that 
one is in a double sandwich and the other is in a triple sandwich. 
Close to identical performance should be expected. At Reynolds 
numbers less than 1000 this is certainly the case, and the discrep- 
ancies noted can be easily attributed to manufacturing tolerances. 
These are compact surfaces and a 0.001-inch variation in fin spac- 
ing or plate spacing can have a rather large effect, especially on 
the friction factor. However, the 15 per cent difference in 
friction factor at the higher Reynolds numbers is more than one 
would like to charge to dimensional tolerances, and the triple 
sandwich surface tends to yield lower heat transfer than its 
double-sandwich counterpart. The fin conduction problem in the 
case of the triple sandwich is considerably more complex than for 
the double sandwich, but this has been taken into consideration. 
The results of an analysis of the fin effectiveness in a triple sand- 
wich surface are presented in the Appendix. Besides which, the 
conduction problem can have nothing to do with the friction 
factors. One possible explanation could be a flow mal-distribution 
resulting in a disproportionate share of the flow going through 
the center of the three passages. Such a flow distribution could 
be caused by the abrupt contraction entrance to the test core 
where a two-dimensional vena-contracta undoubtedly occurs, re- 
sulting in a higher total pressure at the entrance to the center 
section than at the entrance to the two outside sections. Calcu- 
lations based on an assumed contraction ratio indicate that a very 
large increase in apparent friction factor could occur, but that 
the apparent Ns.N’p,’/* remains virtually unchanged. However, 
the higher flow velocity in the center section would result in a 
higher conductance, and this in turn would decrease the fin effec- 
tiveness. This latter possibility has not been investigated, but 
still remains the most plausible explanation. If poor flow dis- 
tribution is truly the explanation of the differences noted, then it 
should be further noted that the air flow length of the test cores 
was only 23/, inches. In applications where flow length is sub- 
stantially greater than 23/, inches the flow distribution for the 
triple sandwich would be expected to be better, and performance 
approaching that of the comparable double sandwich should be 
anticipated. Also some form of flow distribution control could be 
provided, such as a nozzle entrance. At any rate, this is a prob- 
lem that should be considered in application of the tripie sand- 
wich arrangement to heat exchanger design. 

The results for surface !/,-13.95, Fig. 4, are of interest primarily 
in that they show the effect of a substantially thicker fin than any 
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of the comparable surfaces previously considered in reference [1]. 
The fin thickness is 0.010 inch as compared to the 0.006 inch of 
the otherwise very comparable surface !/,-15.2 previously re- 
ported. In high temperature applications, such as in the gas 
turbine regenerator, the alloy steels generally employed necessi- 
tate a thick fin because the low thermal conductivity of the steel, 
coupled with the characteristic high convection conductance of 
the high-performance surface, results in low fin effectiveness. The 
thicker fin in this case results in a 25 per cent higher friction fac- 
tor, undoubtedly attributable to the increased form drag. Sur- 
prisingly, the heat transfer performance is similarly improved, 
evidently meaning that the boundary layer is strongly influenced 
by vortexes shed from the blunt fin edges. 

The Perforated Fin Surface. Surface 13.95(P), Fig. 5, exhibits 
performance quite unlike the other high performance surfaces. 
A dip in the N sVp,*’* curve in the 1500 to 3000 Reynolds-number 
range suggests the laminar flow-turbulent flow characteristic of 
flow through smooth tubes, although the entire curve is shifted 
upward at least 50 per cent above a smooth tube. The most in- 
teresting feature of the performance of this surface is the friction 
factor curve. The friction factors in the “turbulent” region are 
only a little more than twice the Nsg.Np,'/*, indicating little if any 
form drag. Such a characteristic is very desirable in heat ex- 
changer design, for it results in a minimum heat exchanger core 
frontal area, a particularly bothersome problem in many applica- 
tions of very compact high-performance surfaces, and most 
especially in the gas turbine regenerator application. The holes 
apparently accomplish an effective interruption of the turbulent 
sublayers without the high velocity part of the stream impinging 
upon blunt edges normal to the flow direction. Part of the sur- 
face area is of course lost by this expedient, nullifying part of 
the gain as far as heat exchanger volume is concerned, although 
the weight advantage remains. From the results of these tests 
alone it is not possible to tell whether an optimum has been 
achieved with respect to size of hole and frequency of hole. 

The Wavy-Fin Surface. Surface 11.5-*/;(W), Fig. 6, differs geo- 
metrically only very slightly from a previously reported surface, 
11.44~-%/,(W), in reference [1]. The test results differ no more 
than 5 per cent over the entire test range. While adding nothing 
new, these results do demonstrate remarkable repeatability, con- 
sidering the facts that the test core dimensions and steam-side 
arrangement differed considerably and methods of holding manu- 
facturing tolerances have undoubtedly improved in the eight 
years separating the two sets of tests. 


Summary and Conclusions 

In this paper the basic heat transfer and flow friction charac- 
teristics have been presented for six very compact high-per- 
formance heat transfer surfaces. Four are of the strip-fin variety 
and develop their high performance by the use of short discon- 
tinuous fins. One has fins with holes punched through to provide 
boundary layer interruption, and one has wavy fins which cause 
a boundary layer breakdown at each turn. 

Three of the strip-fin surfaces are constructed in multiple sand- 
wiches, and the tests indicate that flow distribution difficulties 
may be encountered in a triple sandwich arrangement, especially 
if an abrupt contraction entrance is used and the total flow length 
is small. In the Appendix the special problem of the fin effective- 
ness of a triple sandwich is considered. 

One of the strip-fin surfaces demonstrates the results of rela- 
tively thick fins, which tend to substantially increase friction 
factor, but which also materially increase heat transfer per- 
formance. 

The perforated fin surface results indicate a distinct ad- 
vantage in this type of boundary layer interruption, for a sub- 
stantial increase in heat transfer performance is attained without 
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introducing an additional form-drag component to the friction 
factor. 

The results for the wavy fin surface serve primarily to demon- 
strate the repeatability of the data contained in reference [1]. 
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APPENDIX 


The Fin Effectiveness. The effectiveness of a fin is defined as the 
ratio of the heat actually transferred through the base of the fin 
to what that transfer rate would be if the entire fin surface were 
at the base temperature. In any actual fin the effectiveness is 
some number less than unity because of the temperature drop 
along the fin due to conduction resistance. Thus only with a fin 
material thermal conductivity of infinity could the fin effective- 
ness be unity. 

For simple fins the effectiveness is a function of the fin length 
and thickness, the thermal conductivity of the fin material, and 
the convection conductance along the fin surface. These rela- 
tionships are well known and will not be repeated here. How- 
ever, the double and triple sandwich surfaces present a complex 
fin geometry which requires its own special solution. In this 
Appendix a genera! solution applicable to any double or triple 
sandwich arrangement will be presented, and then further con- 
sideration will be given to a fin effectiveness problem peculiar to 
the triple sandwich alone. 

Any double or triple sandwich finned surface, symmetrical 
about the center line of the passage, can be represented for pur- 
poses of analysis by the diagram shown in Fig. 8. To maintain 
generality each of the fin elements are considered to be of different 
length and thickness. The analysis of this system for the fin 
effectiveness follows precisely the same procedures as are em- 
ployed for analysis of simple fins; the problem is simply more 
complex, and the resulting expression for effectiveness is more com- 
plex. If the fluid surrounding the fin is considered to be 


CENTERLINE 
OF MIDOLE 
PASSAGE 


b 
L 


SYMMETRY PLATE 


idealization of the fins in a multiple sandwich surface 


Fig. 8 
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everywhere at the same temperature, if the surface convection con- 
ductance is considered to everywhere be the same, and if fin tem- 
perature is a function only of fin length (thin fins), the following 
expression for the effectiveness results: 


1 
m(L + M +a +b) 
+(B+ D+ di it (1) 
lam, + (B + D + E) tanh (m,L)]} cosh (mL) 


B = Bm: tanh (ma) 
D= tanh (m;b) 
dm, tanh (mM) 


V 2h/ka 

V 2h/kB 

V 2h/ky 

= 2h /kd 
This solution can be reduced to the simple fin solution if a, b, and 
M are set equal to zero, or to the case of the double sandwich sur- 

face if M is set equal to zero. 

Since a part of the total heat transfer surface area is a base or 


prime area (not extended area), an over-all surface effectiveness, 
is defined as, 


where 


= 1 —(1 — 9,(A,/A) (2) 


7, then expresses the ratio of actual heat transferred by the entire 
surface, A, to what this rate would be if 7, were 1.00. Note that 
1. < 

Equation (1), together with (2), is probably an adequate solu- 
tion for most double sandwich surface applications, but for a 
triple sandwich arrangement a new problem arises, An idealiza- 
tion basic to the analysis was that the fluid temperature is the 
same along the entire fin surface. If in a triple sandwich ar- 
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Fig. 9 Effect of the heat exchanger NTU on a triple sandwich over-all 
surface effectiveness 
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rangement the fin effectiveness tends to be low the fluid tem- 
peratures in the center and outer passages will only be close to 
equal near the entrance to the heat exchanger. At points re- 
moved from the entrance there will be less heat transfer to the 
fluid in the center passage, and a temperature difference between 
the fluid in the center passage and in the outer passages will be 
established. . This solution no longer applies, and it is easy to 
see that the apparent fin effectiveness will be lower than given by 
this solution. 

An approximate analysis for this effect has been carried out for 
the case of a heat exchanger with the base surface at a uniform 
temperature. A triple sandwich surface was chosen and was 
simplified such that (referring to Fig. 8) M = L/2,a = 6 = L/4, 
and a = B = 6 = y. This corresponds approximately to the 
proportions of the triple sandwich surface shown on Fig. 3. The 
details of this analysis are too lengthy to present here, but the re- 
sults are of interest. The over-all performance of the heat ex- 
changer was evaluated, taking into consideration the different 
behavior of the center passage as compared to the two outer 
passages. From this, working backward, an over-all surface 
effectiveness 7, was evaluated such that the simple heat ex- 
changer solution would be applicable. Labeling this new effective- 
DESS Nornd approx, ANd the over-all effectiveness from Equations (1) 
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and (2) as Moist approx, the ratio of the second approximation to the 
first approximation is plotted on Fig. 9 as a function of the first 
approximation, with the heat exchanger NTU as a parameter. 
As expected, the triple sandwich surface will have a lower surface 
effectiveness than the more simple solution, Equations (1) and (2) 
would indicate, and this effect will be more pronounced for low 
surface effectiveness and for high NTU. The worst condition 
met in the tests of surface !/;-16.12(T) is a surface effectiveness of 
about 0.90 at an NTU of about 1.0. From Fig. 9 the correction 
factor is 0.996, which is hardly worth worrying about. However, 
in an application of a triple sandwich surface where low conduc- 
tivity high temperature materials must be used, as say in a gas 
turbine regenerator, Moist approx Might be in the range of 0.70 and 
NTU around 3.0. In this case the surface effectiveness as evalu- 
ated from simple theory would be about 7 per cent too high, and 
the true surface effectiveness would be nearer to 0.93 X 0.70 or 
0.65. Note, on the other hand, that this analysis has been based 
on a heat exchanger with one fluid side at a constant temperature 
(or at least the base surface at a constant temperature) as in a 
condenser or evaporator. Somewhat different behavior would be 
expected of a counter or crossflow heat exchanger with equal 
capacity rates, but it is expected that Fig. 9 can still be used as a 
good approximation to the effect. 
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ASME Researcu Commitree on High Tem- 
perature Steam Generation has a program in progress to deter- 
mine the performance characteristics of selected ferritic and 
austenitic tube materials when exposed to the action of steam at 
high temperatures and under pressure. The test units for this 
program are installed at the Philip Sporn Station of American 
Electric Power Service Corporation near Pomeroy, Ohio. In this 
research, selected materials in tubular form are subjected to high- 
pressure steam (2000-psi) at temperatures of 1100 to 1500 F for 
time periods up to 36 months to determine the resistance to cor- 
rosion and creep, the effects of oxide films on thermal conduc- 
tivity, the metallurgical stability, and the resistance to thermal 
shock. The test tubes are exposed to air on the outside and there 
is no temperature gradient across the tube wall. The present 
paper is concerned with the metallurgical aspects of this program. 

Two papers have already been presented before The American 
Society of Mechanical Engineers on certain phases of the ASME 
Research Committee on High Temperature Steam Generation 
program to determine the performance characteristics of selected 
ferritic and austenitic tube materials when exposed to the action 
of steam at high temperatures and under pressure. At the 1951 
Annual Meeting,’ the apparatus constructed for this program 
was described *n detail,® while at the 1955 Annual Meeting, re- 
sults of tests on the first prototype unit were given. A schematic 
arrangement of the test installation and a typical section of a test 
rack are reproduced in Figs. land 2. Five units are involved with 
the two containing the ferritic tube materials operating at 1100 
and 1200 F, and the three with the austenitic tubes at 1200, 1350, 
and 1500 F. The steam pressure in each unit is approximately 
2000 psi and the tubes are 2-in. OD X '/:-in. wall X 42-in. long. 
The internal and external surfaces were as would be normally 
supplied to ASTM Specification A213. The type of steels tested 
and their arrangement within the unit are shown in Fig. 3. 

The tube dimensions were arbitrarily selected. The size was 
considered representative of that generally used in superheaters; 
and all specimens were made identical to facilitate the final scale 
analyses. 
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Metallurgical Evaluation of Superheater- 
Tube Alloys After Six-Months’ Exposure at 
Temperatures of 1100 to 1500 F 


As a result of this arbitrary selection, it is evident that all 
tubes are subject to the same stress and, thus depending upon 
the materials selected and the operating temperatures, the ex- 
pected life of the specimens will vary considerably. 

The pressure stress was established on the basis of the Barlow 
formula 

PD 


where 
s = stress, psi 
P = pressure, psi 
D = inside diameter, in. 
t = tube-wall thickness, in. 


Thus for the steam pressure of 2000 psi established for this re- 
search, the tubular specimens are subjected to a stress of 2000 psi. 
The long-time strength properties of the materials were based 
on data published by The Babcock & Wilcox Company, The 
Timken Roller Bearing Company, and The United States Steel 
Corporation. Two criteria were used in selecting the materials 
for the test rack. First, the expected stress of 2000 psi was com- 
pared with the stress required to produce a creep rate of 0.02 per 
cent in 1000 hr and, secondly, the expected stress was compared 
with 80 per cent of the stress necessary to produce rupture in 
25,000 hr. However, in the final selection, some of the materials 
selected did not meet these criteria, and a calculated risk was 
taken in the decision to include such materials in the test racks. 
After the first six-months’ operation, the first set of tubes, desig- 
nated as Selection A of Fig. 3, was removed from each unit and 
the present paper compares the metallurgical characteristics of 
these tubes with those of the original materials and also presents 
findings with respect to the degree of tube swelling and the mag- 
nitude and type of attack on both the OD and ID surfaces. 


Operating and Test Procedure 


The operating conditions during the first six months of opera- 
tion at the Philip Sporn Power Plant were as follows: 


Test units 


No.2 No.2A No.3 No.4 


1201 1206 1353 1497 


1350 
1341 


1494 
1485 


1194 

Outlet-probe read- 

ing, deg 1202 
Average inlet steam 


830 


*In soot-blower section, average metal-temperature drop 45 F 
during 10-sec blow. 

> The pressure in the 1500-F unit is between 1500 and 2000 psi, 
because of the pressure drop through the test assembly. 
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steam pressure, 
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900 TO 1000°F STEAM IN 


2000 p.s.1.-1000 /bs./ hr. 


>— 


900 TO 1000°F STEAM 
BACK TO SYSTEM. 


LEGEND 


HE HEAT EXCHANGER 
SH* SUPERHEATER. 
DESUPERHEATER. 


PROTOTYPE INSTALLATION 


RACK 
NO. 3 


1500° F. 


Schematic arrangement of test installation 


ROTATING AIR 


BLOWER 


STEAM OUTLET). 


SPECIMEN 


| -MEATING ELEMENT 
MUFFLE SHEET 


oF 


34° Dia 


LENGTH 


24° 


TEST RACK No | @2 


3-6" 


THERMOCOUP! 


4-6" HEADER TO HEACER & 


FERRITIC ALLOYS 
To 1200°F 


| yy 


= = - STEAM INLET 
2°00 1172" THICK TUBE 


TEST RACK No 26,384 
AUSTENITIG ALLOYS 
1200°F, I350°F, 


Fig. 2. Typical section of test rack 


As would be expected for a program of this type, the operation 
for the first six months was not continuous. The boiler furnish- 
ing steam was out of service part of the time, for various reasons. 
There also were failures in the 1500-F test unit. These failures 
will be discussed later. 

The metallurgical examination of the tubes removed was under 
the supervision of the metallurgical subgroup of the committee 
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and was done on a co-operative basis: The testing program 
undertaken on the tubes before and after exposure is given in 
Table 1. This table also lists the co-operating laboratories and 
shows the type of steel assigned to each. All testing details were 
in accordance with ASTM Specification A370. 

Fig. 4 shows a photograph of one of the tube elements with the 
end connections and the section exposed to the soot-blower 
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A) Type Steels Tested In Experimental Units 


Unit No. 1 Unit No. 2 Unit No. 2A Unit No. 3 Unit No. 4 
° 
1100 F 1200°F 1200°F 1350°F 1500°F 


1. Cr-Mo-V TP 304-(18 Cr - 8 Ni) 

2-1/4 Cr - 1 Mo TP 347 (304+Cb 

5 Cr Mo + Ti TP 321 (304+T1 

2-1/4 Cr - 1 Mo* TP 316 (16 Cr - 13 Ni+Mo) 
Mo FP 310 (25 Cr - 20 Ni) 
Mo 16 Cr - 25 Ni - 6 Mo 
Mo* 17 Cr - 14 Ni + Cu + Mo 
- 1 Mo 15 Cr - 15 Ni + Nitrogen 
Mo + Si Inconel {35 Cr - 78 Ni 
Mo + Ti* Incoloy (20 Cr - 32 Ni 


OW ON OVW Sw 


*Welded at center of length with manual deposit. 


Note: 1) Complete chemical analyses and details with respect to heat 
treatment and mechanical properties given in Figures ( ) through ( ) 
of the Appendix. 


2) Tubes are 2.000" OD by 0.500" wall and were produced in accordance 
with ASTM Specification A213 (SA 213) 


B) Arrangement Of Tubes In Unit 


Fig. 3 Type steels tested and arrangement in the test units 


Journal of Engineering for Power JANUARY 1960 / 3] 


3 
A 
© 
~~ 
@ 
vo) 
9 Vv 
8 
SECTION 


action; and a schematic diagram for the location of the specimens 
required for the metallurgical examination. The particular tube 
element shown contains a center weld, and from Fig. 3, it will be 
noted that three tubes with center welds were installed in the 
ferritic-tube units. 


Test Results 


The results obtained from the metallurgical examination 
before and after the six-months’ exposure at the designated tem- 
peratures are summarized for each of the individual] steels in 
Tables 4 through 20 of the Appendix. Similar data with ref- 
erence to each of the individual failures to date are given in Figs. 
12 to 16, inclusive of the Appendix 

Visual Examination. Representative outside surface appearances 


Table 1 Outline of tests conducted on tubes before and after exposure 


I Tests on original tubes 
A Tensile tests (duplicate 0.357-in. diam X 1'/:-in. gage 
from tube wall) 
1 Room temperature 
2 temperatures (1100, 1200, 
1300, 1500 F) 
Charpy impact tests (0.394-in-sq V notch) 
m temperature 
2 Proposed operating temperature 
Hardness (Brinell on cross section) 
Microexamination 
1 Structure 
2 Decarburization 
3 Grain size (structural and McQuaid-Ehn) 
) Flattening test (to complete flatness) (1-in. ring) 
F Size measurements 


ya 


II Tests on tubes after exposure 
Based upon a review of the results obtained after the six- 
month period, a decision will be made on the completeness 
of the testing required after the 18, 24, and 36-month 


periods. 
A All tests as for original tubes 
B Examination of racks after removal from unit 
C Visual examination of general surface conditions—4- 
in. length 
1 Incipient cracking 
2 Scaling at soot-blower area 
D Scale thickness and nature of attack (micro) 
1 Intergranular versus general attack 
E Heat-transfer tests—Professor Solberg at Purdue 


(a) Superheater tube showing end 


of two ferritic tubes after the exposure in air at 1100 and 1200 F 
are shown in Fig. 5 and of two austenitic steels after exposure at 
1200, 1350, and 1500 F in Fig. 6. Unfortunately, surface details 
can not be shown too clearly by photographs and accordingly the 
following comments as to the conditions are included. 

Visual inspection at the outer surface cf the 1100-F test unit 
indicated the titanium and silicon modifications of the 5-Cr Mo 
steel scaled the least amount while the 9-Cr 1-Mo steel was second 
best. The most heavily scaled appeared to be the 2'/,-Cr 1-Mo 
and the Cr-Mo-V tubes. However, in no case was the scaling 
very heavy and subsequent microscopic examination revealed 
it to be of the order of 0.005 in. 

Inspection of the 1200-F ferritic test unit again indicated the 
titanium modification of 5-Cr Mo steel scaled the least amount 
but the 9-Cr 1-Mo steel appeared to be slightly superior to 5-Cr 


III Assignment of work to co-operating laboratories 
International Nickel 


Company, Inc. Inconel —3 tubes 
Incoloy —3 tubes 
Type 304 —3 tubes 
Type 310 —3 tubes 

Babcock & Wilcox Com- 
pany 21/,-Cr 1-Mo —2 tubes 
9-Cr 1-Mo —2 tubes 
15-15N —3 tubes 
321 —3 tubes 
National Tube Division 5-Cr '/.-Mo —2 tubes 
United States Steel 3-Cr 1-Mo —2 tubes 
Corporation Type 347 —3 tubes 
Type 316 —3 tubes 

Cr-Mo-V —1 tube 

Timken Roller Bearing 
Company 5 Cr-Mo-Ti —2 tubes 
Steel and Tube Division 5 Cr-Mo-Si —2 tubes 


16-25-6 —3 tubes 
Armco 17-i4 Cu-Mo—3 tubes 

Naval Experiment Sta- 
tion All welds (seven tubes plus 


and header 
welds) 
Purdue University All tubes to be tested for heat- 


transfer properties. 24-in. 
piece required. These will 
then be available to the 
above laboratories for ten- 
sile, impact, and other teste. 
General Electric Com- 
pany Analysis of scale 


oti oti. 


A is soot blower section 


(b) Schematic diagram for locati 


of test specimens after exposure 
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Fig. 4 Photograph of superheater element and sketch showing test specimen location 
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A Height of scale deposit in tubes, in. 


Table 2 Scaling characteristics of superheater-tube materials based on scale accumulation in tubes 


Height of scale deposit in tubes ~ 

1350 F oh 1500 ~ 
Type alloy Section B Section C Section B Section C 
om: Tubes removed 
17-14 Cu-Mo 31/5 
15-15N... 9'/, 7/3 


* No measurable scale was detected by radiographic examination of tubes where no value is shown. 
In the radiographic examination, exposure was determined for 1-in. total thickness, 18-in. focal 


distance, 200 mg radium, aiming for 
B_ Analyses of the scale samples 


Temp, 
Type alloy deg F_ Al Mo Cu 
15-15 1200 0.01 0.04 0.05 
1350 0.01 0.38 0.06 
1500 0.02 0.60 0.15 
304 1500 0.01 0.09 0.26 
316 1500 0.02 0.85 0.18 
17-14 Cu-Mo 1350 0.01 0.09 0.20 


C Magnetic and microscopic examination 


m density of about 1.5. 


Si Cr Cb Ti Ni Mn 
0.03 0.22 .*. %O.01 0.90 0.60 
0.01 .*. 0.90 
0.15 12.50 0.35 0.03 8.50 0.94 
0.11 14.50 .*. 0.01 6.20 1.05 
0.13 19.00 .*. 0.01 8.70 1.30 
0.01 1.00 .*. 0.03 2.30 0.52 


In all ten samples examined, the scale appeared homogeneous, black or gray-black in color, 


Mo + Si. The most heavily scaled appeared te be the 5-Cr Mo, 
2'/.-Cr 1-Mo, and Cr-Mo-V tubes. Evidently the nature of the 
scale formed, that is, whether it is dense and adherent or tends 
to be flaky, greatly affects the visual appearance. Microscopic 
examination showed the thickness of scale to range from 0.010 
to 0.017 in. with the 9-Cr 1-Mo showing the least, and the 21/,-Cr 
1-Mo tube the greatest thickness. 

The only distinguishing feature of the tubes in the austenitic 
units was the fact that the 16-25-6 tube in the 1500-F unit was 
heavily scaled and pitted around the weld at the inlet end and also 
exhibited a number of slight blisters over its length. The pitting 
might have resulted from weld slag, incompletely removed. 
Microscopic examination showed the scale thickness on the out- 
side of these tubes to be the order of 0.001 in, or less. 

The inside surface of two representative ferritic tubes is shown 
in Fig. 7 and of two austenitic tubes in Fig. 8. Insofar as the 
ferritic tubes are concerned, the scaling on the inside surface re- 
sulting from steam attack was of the same order of magnitude as 
the air attack on the outside. Likewise, the scale formed must 
have been adherent for no appreciable amount of loose scale was 
found in any of these tubes. 

While microscopic examination indicated the scale formation 
on the inside surface of the austenitic tubes to be of a low order 
of magnitude, it was found that certain of these tubes contained 
appreciable amounts of loose scale when the units were dis- 
mantled. Consequently an x-ray examination was made of the 
tubes in Section B and C to determine the height of the scale de- 
posits. These results are summarized in Table 2. The tube 
showing the maximum amount of scale was 15-15N at both 1350 
and 1500 F. The only other tube showing measurable scale at 
1350 F was 17-14 Cu-Mo while at 1500 F the 18-8 tube in Section 
C scaled appreciably. The analyseg of the scales, also given in 
Table 2, show the scale formed at 1500 F came from the tubes in- 
dicated and was not carried over from one tube to another. No 
ready association can be made with the scale formed at the other 
temperatures. Before the units were restarted for the second 


six-months’ exposure, the scale was removed from Sections B and 
© through a drilled hole at the base of the tube containing scale. 
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flaky, and hard. The samples were magnetic, indicating the main constituent present to be iron. 


Scale-thickness measurements made on both the air and 
steam side of two of the ferritic and two of the austenitic tubes 
are given in Fig 9. This shows that there is no significant dif- 
ference in the scaling resistance in steam and air at 1100 F be- 
tween 2!/,-Cr 1-Mo and 9-Cr 1-Mo, whereas at 1200 F there is 
a considerable difference between these two materials. This 
figure also shows type 321 is more resistant to scaling than 15-15N 
in both air and steam at 1200 to 1500 F. 

Tensile Properties. Room-temperature and elevated-temperature 
tensile tests at the exposure temperatures were conducted on the 
original tubes and also after the six-months’ exposure time. The 
results, for each of the individual steels, are contained in Tables 4 
through 20 of the Appendix. 

The exposure conditions had little influence on the tensile or 
yield-strength values at either room or elevated temperatures but 
did have an influence on the ductility values of certain of the 
steels, especially at room temperature after exposure at the ele- 
vated temperatures. This condition is illustrated by Fig. 10, 
which is based on the elongation values. The two steels showing 
the greatest loss in room-temperature ductility after exposure at 
elevated temperatures are type 310 and 16-25-6. The high- 
temperature exposure likewise reduced the elongation values at 
elevated temperatures of certain of the alloys, and especially 
Incoloy and 17-14 Cu-Mo at 1350 F and 17-14 Cu-Mo and type 
310 at 1500 F. 

Impact Resistance. The six-morths’ exposure time had a marked 
effect on the shock resistance of certain of the alloys, both at room 
temperature after exposure at elevated temperatures and at the 
elevated temperatures after exposure. This condition is shown 
in Fig. 11. The alloys showing the greatest decrease in room- 
temperature impact resistance after exposure at elevated tem- 
peratures were 16-25-6, type 310, 15-15N, and 17-14 Cu-Mo. 
These same alloys likewise showed the greatest loss in high- 
temperature impact resistance after the six-months’ exposure. 
For both of these conditions, Incoloy showed its minimum re- 
sistance after the 1200-F exposure. 

Flattening Characteristics. All of the tubes considered met the 
flattening requirements of ASTM Specification A213 in the origi- 
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Fig. 5 Photograph showing outside surface of representative ferritic tubes after six-months’ 
exposure at 1100 and 1200 F 4 

(B11, B21, 5-Cr Mo Ti, after 1100 and 1200 F, respectively. F11, F21, 5-Cr Mo Si, after 1100 
and 1200 F, respectively.) 
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Fig. 6 Photograph showing outside surface of representative austenitic tubes after six-months’ exposure at 1200, 1350, and 1500 F 
(M2, M3, M4, 16-25-6, after 1200, 1350, and 1500 F, respectively. N2, N3, N4, 17-14 Cu-Mo after 1200, 1350, and 1500 F, respectively.) 
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(a) 2'/,-Cr 1-Mo, 1100 F 


(b) 2'/,-Cr 1-Mo, 1200 F 


(d) 9-Cr 1-Mo, 1200 F 
Fig. 7 Photograph showing the inside surface of representative ferritic alloy tubes after six-months’ exe 


posure at 1100 and 1200 F 


nal condition and, with but few exceptions, the same was also 
true after the six-months’ exposure at the various temperatures. 
In fact, most underwent complete flattening without showing 
any signs of cracking. 

The exceptions were 16-25-6, which could not be flattened to 
the calculated H value after any of the exposure temperatures, 
and types 310 and 347 which would not meet this requirement 
after the 1500-F exposure. The silicon modification of 5-Cr Mo 
satisfactorily met the H value but could not be completely flat- 
tened without cracking. 

Bulging During Exposure Tests. Measurements of the tubes after 
the exposure tests, as well as the calculated per cent deformation 
of the OD’s are given in the figures of the Appendix and scme are 
included in Table 3. Table 3 also includes the available creep 
(0.01 per cent per 1000 hr) and 100,000-hr rupture-strength 
values, as taken from the literature, for the ferritic steels at 1100 
and 1200 F and the austenitic steels at 1350 and 1500 F. 
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The values of Table 3 would indicate that there is no definite 
relationship between the amount of bulging and the reported 
high-temperature properties. With the ferritic steels at 1100 and 
1200 F a certain correlation exists in that the alloy with the maxi- 
mum rupture strength underwent the least amount of bulging 
and at 1200 F the steel with the lowest creep and rupture strength 
underwent the most. However, all of the remaining steels did 
not align themselves in proper order. With the austenitic steels, 
the correlation was less definite. : 

The lack of a more definite correlation may be a result of several 
factors such as: (a) The magnitude of the values involved, (b) 
uncertainties with respect to the actual diameter of the original 
tubes, (c) possible scale effects influencing the readings, and (¢) 
the high-temperature data used for comparison not being repre- 
sentative of the test specimens. A more definite correlation may 
exist after the second six-months’ exposure. 

Failures During Test. To date several tubes have failed by ruptur- 
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(a) Type TP321, 1200 F 


(b) Type TP321, 1350 F 


(c) Type TP321, 1500 F 


(f) Type 15-15N, 1500 F 


Fig. 8 Photograph showing the inside surface of representative and 
austenitic alloy tubes after six-months’ exposure at 1200, 1350, and 1500F 


ing, or excessive bulging, but all, with one exception, have been 
confined to the 1500-F test unit. Likewise there was one failure 
of a section of the type 316 connecting pipe Details with respect 
to these failures are given in Figs. 12 to 16 of the Appendix. 

The one failure not in the 1500-F test unit was one of the 5-Cr 
Mo + Ti tubes tested at 1200 F. As shown in Fig. 12 this failure 
resulted from a marked defect on the ID of the tube and was not 
caused by this alloy having insufficient strength under the given 
test conditions. The room-temperature mechanical properties 
of the failed tube were identical to those of the original material. 
The remaining two tubes of this analysis are still under test. 

Failure by cracking of one of the type 347 tubes occurred at 
2178 hr and the other two tubes of this analysis were bulged to 
such an extent (about 7 per cent) that they were removed from 
the unit. 

One of the type 310 tubes failed by cracking after about 3000 
hr. From Fig. 13 it will be seen that a series of cracks occurred, 


Journal of Engineering for Power 


that sigma was formed, and that there was a marked pick-up of 
nitrogen on the outer surface. Likewise the room-temperature 
physical properties were changed, with both the ductility and im- 
pact resistance being greatly reduced. The remaining two tubes 
of this analysis were also removed because of excessive bulging. 

A type 321 tube failed, by cracking, after 6812 hr, Fig, 14. 
In sections removed from the actual failure, there was considera- 
ble cracking on both the OD and ID and, as in the case of the 
type 310 tube, marked nitrogen pick-up on the outer surface. 
The deformation to failure was 5.0 per cent and the room tem- 
perature properties were not appreciably influenced. The re- 
maining two tubes are still under test. 

A 17-14 Cu-Mo tube likewise failed by cracking after 6812 hr, 
Fig. 15. The deformation to failure was 3.8 per cent and room- 
temperature impact resistance was decreased. The nitrogen 
pick-up at the surface was not as pronounced as with the type 
310 and 321 tubes. The remaining two tubes are still under test. 

The Incoloy and Inconel alloy tubes were removed from the 
1500-F unit because of excessive enlargement. One tube of each 
was reinstalled to verify the findings on these compositions. 

A section of the type 316 connecting pipe failed, Fig. 16, after 
6%12 hr in a very ductile manner, its elongation to fracture being 
398 per cent. The grain size was extremely fine and there is a 
question as to whether or not this pipe was heat-treated so as to 
develop maximum strength at elevated temperatures. Large 
amounts of the sigma phase were formed and there was a marked 
pick-up of nitrogen at the outer surface. 

In addition to the preceding failures resulting from insufficient 
strength, or in one case to a defective tube, a rather unusual 
failure occurred in 16-25-6 as a result of the thermocouple’s pro- 
tective coating (Sauereisen cement) reacting with and eating 
through the tube wall. This hole was repaired by welding and 
this tube is still under test. 

There appears to be no correlation between the time of these 
failures and the reported high-temperature strength character- 
istics given in Table 3. This probably would not be true if high- 
temperature data were available on the actual tubes used in these 
tests. The type 347 and 310 tubes in the 1500-F unit failed in 
approximately the same time. Yet type 347 has one of the lowest 
rupture strengths at 1500 F and type 310 about the same strength 
as type 316, yet none of the type 316 elements have failed. Like- 
wise 17-14 Cu-Mo has the highest rupture strength and one of 


' these tubes failed in 6812 hr, while none of the type 304 tubes 


have, even though the reported rupture strengh for type 304 is 
much lower, 

The lack of agreement between the failures and the high-tem- 
perature data may result in part from differences in heat-treat- 
ment. Unfortunately, exact data are not available with respect 
to the heat-treatment given the tubes and it is now known that 
heat-treatment may have a marked effect on the high-tempera- 
ture-strength characteristics. However, the microstructures of 
the tubes that have failed, with the exception of the type 316 
connecting pipe, would indicate that a satisfactory temperature 
was used for solution treating. All the tubes were heat-treated 
in accordance with standard commercial practice for ASTM 
Specification A213. 


Conclusions 

The results from the visual and metallurgical examination 
after the first six-months’ exposure permit the following conclu- 
sions: 

1 Visual examination indicated differences in the appearance 
of the scale on the outside surface of the ferritic tubes after ex- 
posure at 1100 and 1200 F but microexamination indicated rela- 
tively slight differences in the amount of actual attack, especially 
at 1100 F. At 1200 F the measured depth of scale for these 
steels varied from 0.010 to 0.017 in. with 9-Cr 1-Mo showing the 
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Fig. 9 Steam oxidation of some superheater materials after si tion (ASME steam research) 


Table 3 Comparison of maximum tube bulging during exposure tests with high-temperature characteristics 
A Ferritic steels 
————-1100-F exposure tests -———--1200-F exposure tests——— 
Creep Creep 
strength 100,000-hr strength 100,000-hr 
Max 0.01%/1000 rupture Max 0.01%/1000 rupture 
exp. % hr strength exp. % } strength 
2'/4,-Cr 1-Mo 0.15 3600 200 2 4500 
0.80 3600 5400 1.30 2 wink 
§ 0.90 3200 6700 1.50 3000 
5-Cr Mo+ Ti 0.60 3500 4000 3.80 1925 
5-Cr Mo+Si 2600 5400 1.70 j 3100 
9-Cr 1-Mo... j 4100 7800 0.70 4600 


B Austenitic ¢ 
exposure tests——— —-—1500-F exposure tests——— 
2150 3150 1.26 1450 1700 
3500 3600 2.76 1250 1600 
3500 5500 Stes 1800 1600 
2500 2500 te 625 1250 
3000 5000 7.00 625 1250 
sents 9500 0.50 5500 
16-25-6...... 5 te 6250 0.70 3000 3300 
Inconel 
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Fig. 10 Short time elongation values of austenitic steels before and after exposure at 1200, 1350, and 1500 F 


| 
ROOM TEMPERATURE IMPACT 
AFTER 6 MONTHS EXPOSURE 


a 


8 


& 


8 


> 
g 


& 


Fig. 11 


least and 2!/,-Cr 1-Mo the greatest amount of scaling. At both 
temperatures the scale thickness on the ID was comparable to 
that on the OD. This is in agreement with the known effect of 
increasing chromium content in improving the resistance to 
scaling in both air and steam, especially at temperatures in excess 
of 1100 F, 

2 With the austenitic steels the scale on the OD was slight 
at each of the three temperatures. While visual appearance of the 
ID indicated little attack, certain of the tubes, especially 15-15N 
at 1350 and 1500 F, contained several inches of flaky scale, thus 
indicating appreciable attack by the steam. 

3 The mechanical properties of the ferritic steels were not 
appreciably changed during exposure at 1100 and 1200 F. The 
same was also true of the austenitic steels with respect to their 
tensile and yield-strength values. However, with certain of the 
austenitic steels, a marked decrease occurred in the room-tem- 
perature-ductility and impact properties as a result of exposure 
at the elevated temperatures. The steels undergoing the great- 
est changes were 16-25-6, type 310, 15-15N, and 17-14 Cu-Mo. 
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Impact characteristics of austenitic steels before and after exposure at 1200, 1350, and 1500 F 


4 A definite correlation was not found to exist between the 
high-temperature rupture properties and the degree of bulging 
or, in the case of the tubes that have rupture, in the rupture time. 

5 To date failures have occurred in the 1500-F tests of type 
347 and 310 and the remaining tubes of these analyses were re- 
moved because of excessive bulging. Incoloy and Inconel alloys 
were also removed because of excessive bulging. One tube of type 
321 and 17-14 Cu-Mo have likewise failed but the two remaining 
tubes of each of these analyses are still under test. 

6 All of the austenitic tubes that failed in the 1500-F test 
unit showed a marked nitrogen pick-up on the surface and this 
probably accounted for some, at least, of the shallow-OD crack- 
ing visible on the macro sections. There is a question if this same 
condition would occur in combustion gases. 

These conclusions apply only to the given test conditions. 
With higher rates of steam flow through the tubes, or with com- 
bustion gases on the outside of the tubes, different results may 
have been obtained 
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APPENDIX 


of failure 


(c) microstructure nital etch 


(b) Tube sections showing type failure 


(d) Room temperature properties before and after failure 


Tensile longation Reduction Maximum OD 
strength, (0.2% in, 15", of area, Brinell Charpy V expansion, 
psi psi % % hardness notch, ft-lb 
Original tube 60,000 25,750 46 76 115 238, 236, 244 
Failed tube 59,000 25,000 43 85 115 224, 237, 254 


Fig. 12 Data on 5-CrMo Ti tube which failed in the 1200 F superheater unit 
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(a) Photograph showing appearance 
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(a) Photograph showing appearance of the failure 


(c) Microstructure showing sigma formation, (d) Microstructure showing nitrogen pick-up at 

100 surface; nitrogen content 0.394 per cent, tube 
proper, 0.058 per cent, 500, chromic acid, 
electrolytic etch 


(b) Macrostructure showing extent of cracking 


(e) Room temiperature properties before and after failure 
Yield Elongation Reduction 
(0.2%), in 1.5", of area, Brinell 
psi % % hardness 


Original tube 32,750 50 74 128 
Failed tube 82,000 35,000 33 31 143 


Fig. 13. Data on 310 tube which failed after about 3000 hr at 1500 F 
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4 
Maximum OD 
Charpy V expansion, 
notch, ft-lb % 
206, 204, 206 Pye 
20, 19, 18 9.75 | 


(b) Macrostructure showing extent 
of cracking 


(c) Microstructure of failed tube GS 5/8 (4), ng cold work in grains, 
xX 100 : X 1000, 75 per cent nitric electrolytic etch 


(e) Nitrogen pick-up at surface 


Original —0.017% 
After exposure —0.210% 


(f) Properties before and after failure 


Tensile Yield Elongation in Reduction Maximum OD 
Test strength, (0.2%), of area, Brinell Charpy V expansion, 
temp. psi psi % % hardness notch, ft-lb % 
Room 


Room 76 , 600 31,000 52.0 73.0 124 165, 216 5.0 
22,000 14,000 41.0 39.0 


Fig. 14 Data on type 321 tube which failed after 6812 hr at 1500 F 
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(2) Photograph showing appearance of the failure 
Failed 
Failed 
ith 
> 
s 


(b) Macrostructure showing extent of 
cracking, 0.10 per cent N: at surface 


(d) Properties before and after failure 


Test 
temp. 
Original Room 
Failed Room 
Original 1500 
Failed 1500 


Fig. 15 Data on 17-14 Cu Mo tube which failed after 6812 hr at 1500 F 
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Maximum OD 
arpy expansion 
notch, ft-lbs % y 
132 
40, 36, 40 3.8 
81, 80, 72 
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(a) Photograph showing appearance of the failure 
X 1000 
(c) Microstructure of failed tube, glyceregia etch 
Tensile Yield Elongation Reduction 
strength, (0.2%), in 1.5", of area, Brinell 
psi psi % % hardness 
89, 250 45,500 46.8 67.7 143 
84,000 38, 500 38.0 51.0 143 
27 ,600 22,300 53.0 58.0 
22,000 17,000 12.0 15.0 
| 


(a) Photograph showing appearance of the failure 


A 


(b) Microstructure showing type of failure 
(c) Microstructure of failed tube; note appearance of sigma phase, 75 per cent nitric electrolytic etch 


(d) Nitrogen pick-up at surface 


Original —0.035% 
After exposure —0.300% 


(e) Properties before and after failure 


Tensile Yield Elongation Reduction Maximum OD 
strength, (0.2%), in 1.5”, of area, Brinell Charpy V expansion, 
psi psi % % hardness notch, ft-lb % 
Original 
Failed 40,000 39.0 48.0 143 51, 61 39.8 


Original 
Failed 16,500 67.0 66.0 


Fig. 16 Data on type 316 connecting pipe which failed after 6812 hr at 1500 F 
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Table 4 Characteristics of Cr, Mo, V before and after six-months' exposure at 1100 F 


Chemical composition 


C Mn 
0.11 0.44 


Heat-treatment 
Annealed 


0.012 


Mo 
0.86 


S Si Cr 
0.023 0.24 1.01 


———Grain size————~ 
Before After exposure 
exposure 1100 F 
7 


7 


Carburized 
grain size 
9/8 


Short-time tensile and impact properties 


Tensile 


Test strength, 
temperature i 


Room (1) 
2 


(2) 

1100 (1) 
(2) 

(1) Unexposed, 

Tube dimensions 

Exposure 

temperature 

Unexposed 2.002/2 .007 

1100 F 2.016/2.016 


* OD at soot blower 
Flattening characteristics 


Exposure 
temperature 


Unexposed 
1100 F 


Scaling characteristics 


Exposure temperature 
1100 F 


Microstructures (etchant: Nital) 


1000 D 
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-———Inlet end, in—— 


0.2% yield, Reduction Brinell 


Charpy 
",% ofarea,% hardness V notch, ft-lb 
81 


62 


183 
163 


30 
31 
39 


(2) 6 months at 1100 F. 


Max 


——Outlet end, in-—— 
Wall expansion, % 


Wall OD OD* 


0.45 


Calculated 


Height at which cracking occurred 
H value 


Inlet end Outlet end 
No cracking No cracking 
No cracking No cracking 


——_Sceale thickness———. 
OD ID Type attack 
General 


0.007 0.007 


4.4.5 


6 eaethes 1100 F 6 — 1100 F 
D D 
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0.22 
95 ,500 62,700 3 
56,700 
Unexposed 
D 


Table 5 Characteristics of 2 1/4-Cr 1-Mo before and after six-months’ exposure at 1100 and 1200 F 


Chemical composition 

Mn 8 P Si 
0.11 0.40 0.021 0.008 0.28 
0.10 0.39 0.018 0.009 0.29 


——— Grain size after exposure ~ 
Original 1100 F 1200 F 

Structural if 6 6-7 

McQuaid Ehn 67 67 6-7 


Short-time tensile, hardness, and impact properties 
Tensile 


Before exposure 
After exposure 


Heat- 
treatment 
Annealed 


ID 


strength, 0.2% yield, Elongation Reduction 


Test temperature psi in 1.4", % area, % 
Room 


temperature 


psi 


42,000 
33,300 
30,200 
22,000 
28 ,500 
18,000 
23 , 000 


(1) 
(2)¢ 


1100 F 
1200 F 


(1) Unexposed, (2) six months at 1100 F, (3) six months at 1200 F. 
Tensile specimens for the other 


Standard 0.357-in. diam X 1.4-in. gage tensile specimens. 
tests were as shown in Fig. 1. 


Tube dimensions 
Exposure 
temperature 
Unexposed 
1100 F 
1200 F 
Flattening characteristics 
Exposure temperature 


Unexposed 
1100 F 
1200 F 


Scaling characteristics 


OD Wall 
2.009/2 .007 0.553/0.547 
2.012/2.008 0.572/0.551 
2.023/2.028 0.570/0. 558 


Height at which cracking occurred 
No cracking 
No cracking 
No cracking 


Scale thickness, in. 
Exposure temperature ID ID 
1100 F 
1200 F 


0.005 0.005 


0.017 


Microstructures (etchant: Picral) 


6 months, 1100 F 
Inside surface 
x 100 
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OD. 


Outlet end, in —-————— 


Cr Mo 


2.11 1.00 
2.12 1.05 


Decarburization——~ 
Original 1100 F 1200F 


0.001 0 
0.001 0.002 


Brinell Charpy 


hardness V notch, ft-lb 


Max % 
expansion 


0.15 
0.95 


Calculated H value 


Type of attack 


General 
General 


6 months, 1200 F 
Inside surface 
«100 


Transactions of the ASME 


. 

0 

66 ,900 3 74.4 143 224, 195 

ang 65 , 200 3 82.7 145 181, 194 

63,000 3 82.6 137 194, 193 

(1) 41,500 4) 81 114, 98 

Ta (2) 41,500 4 83.9 Rac 122, 130 

(1) 27,500 6 85.6 131, 130 

| (3) 32,700 60 89.1: see 125, 112 

Lite 

1.68 

KS 1.69 

1.70 

Unexposed 

100 

“ite 


Table 6 Characteristics of 3-Cr 1-Mo before and after six-months’ exposure at 1100 and 1200 
Chemical composition 

Cc Mn 8 Si Cr 

0.10 0.40 0.013 0.022 0.30 2.93 


Heat-treatment After exposure, 
Annealed Before exposure 1100 F 1200 F 

8 9 9 


Short-time tensile and impact properties 
Tensile Reduc- 

Test strength, 0.2% yield, tion of Brinell Py 
temperature psi psi 1.5”, % area,% hardness V Bere t-lb 
Room (1) 66,000 44,000 120 + 

(2) 65,000 33 ,000 120 + 

64,000 33,300 120 + 
1100 F (1) a 120 + 
2 


(2) 
1200 F (1) i. 120 + 


(1) Unexposed, (2) 6 months at 1100 F, (3) 6 months at 1200 F. 


Tube dimensions 

Exposure Inlet end, in——— -———Outlet end, in ——~ Max % 
temperature OD Wall OD Wall OD* expansion 
Unexposed 2.001/2.006 

1100 F 2.016/2.020 0.7 
_ 1200 F 2.017/2.030 1.2 

* OD at soot blower. 


Flattening characteristics 
Exposure Height at which cracking occurred Calculated 
temperature Inlet end Outlet end H value 
Unexposed No cracking No cracking 
1100 F No cracking No cracking 
1200 F No cracking No cracking 


Scaling characteristics 
Seale thickness, in. 
Exposure temperature OD ID Type attack 
1100 F 0.008 0.008 General 
-1200 F 0.016 0.016 General 


Microstructures 


Unexposed 6 months, 1100 F 6 months, 1200 F 
1000 D 1000 D D 
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Carburized 
grain size 
9/8 
| 
3 17 , 400 78 

Unexposed 
x10 D 


Table 7 Characteristics of 5-Cr 1/2-Mo before and after six-months’ exposure at 1100 and 1200 F 


Chemical composition 
Cc Mn Si Cr Mo 
0.09 0.42 0.025 0.024 0.30 4.93 0.58 


—Grain 


Heat-treatment Before exposure After exposure Carburized 
Annealed 1100 F 1200 F grain size 
5/4 7/5 6/4 9/8 


Short-time tensile and impact properties 
Test Tensile strength, 0.2% yield, Elongation Reduction Brinell Charpy V 
: 


temperature psi psi in 1.5", % ofarea,% hardness notch, ft-lb 
~ Room (1) 68, 500 35,500 37 140 120 + 
’ (2) 64,500 31,000 36 134 120 + + 
(3) 61,000 29, 500 38 126 120 + 
1100 F (1) 23, 200 62 120 + 
(2 23,900 65 
1200 F (1) 15,706 75 120 + 
(3) 15,800 76 


a ) Unexposed, (2) 6 months at 1100 F, (3) 6 months at 1200 F. 
a 


Tube dimensions 


Exposure Inlet end, in.-—— ——Outlet end, in-——— Max % 
temperature OD Wall OD Wall OD* expansion 
Unexposed 2.000/2 .004 
1100 F 2.017/2.020 0.8 
1200 F 2.028/2 .032 1.4 
* OD at soot blower. 
Flattening characteristics 
Exposure Height at which cracking occurred Calculated 
temperature Inlet end Outlet end H Value 
‘ Unexposed No cracking No cracking 
; 1100 F No cracking No cracking 
wa) 1200 F No cracking No cracking 
Scaling characteristics 
hy Exposure temperature Scale thickness, in. Type attack 
OD ID 


1100 F 0.010 0.006 General 
1200 F 0.015 0.017 General 


Microstructures (etchant: Nital) 


Unexposed 6 months, 1100 F 6 months, 1200 F 7 
1000 D 1000 D 1000 D 
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Unexposed 

x 100 D 
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Table 8 Characteristics of 5-Cr Mo -+- Ti before and after six-months’ exposure at 1100 and 1200 F 
Chemical composition 
C Mn P 8 Si Cr Ni Mo Ti 
0.06 0.46 0.019 0.015 0.26 5.02 0.20 0.61 0.57 
Heat-treatment | _——Structural grain size after exposure——~ Carburized 
1100 F 1200 F ‘grain si 


Annealed Original grain size 
6/7 (8) 6/8 6/7 (8) 6/7 


Short-time tensile, hardness, and impact properties 
Test Tensile strength, 0.2% yield, sinention Reduction Brinell Charpy V 


temperature psi psi 1.5”, % of area, % hardness aotch, ft-lb 
Room (1) 59,750 25,750 45.7 75.2 115 188 

(2) 58,300 26, 000 42.5 80.9 116 136, 135, 121 

57,900 114 118, 101, 111 


25,500 99, 99, 95 

(2) 25, 650 143, 139, 143 

1200 F (1) 17, 950 11,450 74.0 93.0 89, 91, 93 
(3) 145, 143, 214 


18,300 
~() ‘Unexposed, (2) after 6 months at 1100 F, (3) 6 months at 1200 F. 


Tube dimensions 
Exposure ————Inlet end, in. — ——Outlet end, in. 
temperature OD Wall OD Wall 
Unexposed 1.998/2.002 0.497/0.503 
1100 F 2.008/2.014 0.549/0.560 2.004/2.008 0.550/0.564 2.005 0.6 
1200 2.007/2.079  0.524/0.536 2.015/2.045 0.540/0.551 2.040 3.8 


* OD at soot blower. 
Flattening characteristics 


Max % 
OD: expansion 


Exposure Height at which cracking occurred, in. Calculated 
temperature Inlet end Outlet end H value, in. 
Unexposed 1.14 1.64 
1100 F 1.24 1.22 1.69 
1206 F 1.14 1.73 


Scaling characteristics 


Scale thickness, in. 
Exposure temperature OD ID Type attack 
1100 F 0.006 0.005 Intergranular, 1 grain 
1200 F 0.015 0.018 Intergranular, 1'/: grains 


Nital) 


Microstructures (etchant: 


6 
x 1000 


x 100 D 


6 months, 1100 F 
D 
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Toble 9 Characteristics of 5-Cr Mo + Si before and after six-months’ exposure at 1100 and 1200 F 


Chemical composition 


Cc Mn P 8 Si Cr Ni Mo 
0.061 0.43 0.019 0.012 1.21 4.80 0.29 0.52 


Heat-treatment Structural grain size after exposure Carburized 
| Annealed Original 1100 F 1200 F grain size 
i 6/7 (5) 8/7 8/7 6/7 (5) * 
- Short-time tensile, hardness, and impact properties 
Tensile 
4 Test strength, 0.2% yield, Elongation Reduction  Brinell Charpy V 
temperature psi psi in 1.5",% of area, % hardness notch, ft-lb 
: Room (1) 71,750 36,500 40.5 78.0 143 156 
(2) 71,000 44,600 36.8 73.1 149 17, 12, 32 
3 69,000 43, 400 36.7 74.5 137 20, 17, 33 « 
1100 F (1) 30,300 18,850 68.5 92.0 122, 124, 122 
: 29,700 18, 600 70.3 88.7 174, 183, 179 
1200 F (1) 19,900 14,000 86.0 95.0 117, 116, 114 
19,700 13,750 73.3 92.7 158.158, 160 


( 1) Unexposed, (2) after 6 months at 1100 F, (3) 6 months at 1200 F. 


Tube dimensions 


bs Exposure Max % 
4 tempera- Inlet end, in. — Outlet end, in.— expan- 
ture OD Wall OD Wall sion 


Unexposed 2.001/2.007 0.501/0.510 
1100 F =. 2.010/2.020 0.555/0.567 2.008/2.013 0.554/0.568 2.010 0.6 
1200 2.037/2.042 0.558/0. 567 2.028/2.030 0.563/0.567 2.028 


a OD at soot blower. 


Flattening characteristics 


Exposure Height at which cracking occurred, in. Calculated 

temperature Inlet end Outlet end H value, in. 
Unexposed No cracking 1.64 
1100 F 1.28 1.62 1.67 


1200 F 


Scaling characteristics 


1.30 1.20 | 


Scale thickness, in. 
I] 


Exposure temperature OD ) Type attack 
1100 F 0.005 0.007 Intergranular, 1 grain 
1200 F 0.014 0.015 Intergranular, 2 grains 


Unexposed Unexposed 6 months, 1100 F 6 months, 1200 F . 
100 D 1000 D 1000 D x 1000 D 
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Microstructures (etchant: Nital) 
Wane 


Table 10 Characteristics of 9-Cr 1-Mo before and after six-months’ exposure at 1100 and 1200 F 


Chemical composition 
0.11 


Before exposure 
0.10 


After exposure 


Mn 8 P 


0.42 
0.42 


0.009 
0.009 


0.021 
0.19 


———Grain size after exposure. 


Heat-treatment 
Annealed Structure 


McQuaid Ehn 


Original 


1100F 1200F 
6-7 6-7 


Not run OD 


Short-time tensile, hardness, and impact properties 


Tensile 
strength, 


Test temperature 


Room 
temperature (1) 
(2)¢ 


1100 F 
1200 F 


Elon- 


0.2% yield, gation 


0.35 
0.34 


Cr Mo 


8.63 0.90 
0.58 0.91 


Decarburization 


Original 1100F 1200F 


ID 0.002 0 0 


0 0 0 


Reduction Brinell Charpy V 


psi in 1.4", % of area, % hardness notch, ft-lb 


130, 165 
104, 116 
123, 127 
166, 170 
149, 146 
140, 142 
139, 129 


(1) Unexposed, (2) six months at 1100 F, (3) six months at 1200 F 


* Standard 0.357 X 14-in. gage tensile specimens. The others were as shown in Fig. 1. 


Tube dimensions 


Exposure temperature 


Unexposed 
1100 F 
1200 F 


Flattening characteristics 


Unexposed 
1100 
1200 F 
Sealing characteristics 


Exposure temperature 


1100 F 
1200 F 


Microstructures 


Unexposed 
Etchant: HCl and 
Picral 
x 100 
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————Outlet end, in — 
OD 


2.002/2.006 
2.011/2.018 
2.016/2.020 


Wall 


0.553/0.559 
0.557/0.561 
0.556/0.564 


Height at which 
cracking occurred 


No cracking 
No cracking 
No cracking 


Scale thickness, in. 
ID OD 


0.008 
0.004 


0.004 
0.010 


6 months, 1100 F 
Inside surface 
Etchant: Picral 

x1 


Max % expansion 


0.60 
0.70 


Calculated 
H value, in. 


1.68 
1.68 
1.69 


Type of attack 
General 


6 months, 1200 F 
Inside surface 
Etchant: Picral 
100 
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Si 
6-7 
EB 71,500 38, 200 35.7 75 148 
72,100 37, 600 31.4 82.4 155 
84,100 59,900 21.4 81.4 148 
(1) 41,500 22,000 45.7 81 shite 
(2) 44,000 22,500 37.1 85 rath 
2 (1) 27,500 18,000 65.7 85.6 es 
(3) 39, 000 30,000 32.9 89.7 ots 


Table 11 Characteristics of TP 304 before and after six-months’ exposure at 1200, 1350, and 1500 F 
Chemical composition 


Cc Mn 8 Si Cr Ni 
of 0.072 1.72 0.012 0.026 0.38 18.48 10.23 
Te Heat-treatment —Structural grain size after exposure—— 
? Solution quenched, 2050 F Original 1200 F 1350 F 1500 F 
1/2 (3, 5) 2/4(1, 5) 1/3 (4) 3/5 (2) 
f Short-time tensile, hardness, and impact properties 
Tensile 
Test strength, 0.2% yield Elongation Reduction Brinell Charpy V 
temperature _ psi psi in 1.5", % of area, % hardness notch, ft-lb 
Reom (1) 83, 500 34,750 65.0 76.0 127 213, 216, 207 
(2) 95,000 40,700 31.0 58.0 163 107, 119, 115 
(3) 91,200 36,300 33.0 62.0 156 120, 120, 122 € 
(4) 88,300 34,800 31.0 60.0 97, 94, 100 
1200 F (1) 46,000 17,000 40.0 58.0 162, 160, 169 
(2) 49,000 21,100 31.0 61.0 181, 181, 177 
1350 F (1)* 43,000 16,250 37.0 50.0 167, 164, 168 
(3) 37,900 16,000 42.0 62.0 178, 169, 181 
1500 F (1) 27 ,000 14,350 41.0 44.0 159, 150, 158 
(4) 41,000 15,800 39.0 65.0 162,147,147 


(1) Unexposed, (2) after 6 months at 1200 F, (3) at 1350 F, (4) at 1500 F. 
21300 F rather than 1350 F. 


Tube dimensions 


Exposure -Inlet end, in. ——Outlet end, in—— Max % 
temperature OD Wall OD Wall OD expansion 
Unexposed 1.977/1.978 0.526/0.560 

1200 F 1.996/1.998 0.553/0.558 1.00 
1350 F 1.996/1.997 0.548/0.557 1.00 
1500 F 1.997/2.003 0.533/0.549 1.26 
Flattening results 
Exposure Height at which cracking occurred Calculated 
temperature Inlet end Outlet end H value, in. 
Unexposed Satisfactory 1.62 
1200 F Satisfactory 1.62 
1350 F Satisfactory 1.62 
1500 F Satisfactory 1.62 


Scaling characteristics 


——Scale thickness———. 


Exposure temperature OD ID Type attack 
1200 F Not reported Assumed very slight 
1350 F Not reported Assumed very slight 
1500 F Not reported Assumed very slight 


Microstructures (x 1000 D) 


\\ 
\ 
\\ 


? 


6 months, 1350 F 6 


Unexposed months, 1500 F 
Etchant: Ferric Etchant: Electrolytic 10 per cent chromic acid 
chloride 
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Table 12 Characteristics of TP-316 before and after six-months’ exposure at 1200, 1350, and 1500 F 


Chemical composition 
Cc Mn Si Ni Cr Mo 
0.06 1.84 0.026 0.014 0.35 13.20 16.77 1.96 
—Grain size 
Heat-treatment Before exposure —After exposure——————. 
Solution annealed 1200 F 1350 F 1500 F 


4/3 1/5 1/5 1/5 


Short-time tensile and impact properties 
Test Tensile strength, 0.2% yield, ey 3 Reduction — Brinell Charpy V 
temperature psi ",% ofarea,% hardness notch, ft-lb 

Room 34,700 

37,700 

38,000 
( 33,300 

1200 F 


1350 F 
1500 F 


(1) Unexposed, (2) 6 months at 1200 F, (3) 6 months at 1350 F, (4) 6 months at 1500 F. 


Tube dimensions 
Exposure ————Inlet end, in-—— ——Outlet end, in—— Max % 
temperature OD Wall OD Wall OD* expansion 
Unexposed 2.003/2.007 
1200 F 2.005/2 .007 
1350 F 2.005/2 .009 
1500 F 2.003/2 .007 


*OD at soot blower. 


Flattening characteristics 


Exposure sa ht oe which cracking occurred 
temperature nlet end Outlet end 
cracking No cracking 
1200 F No cracking No cracking 
1350 F No cracking No cracking 
1500 F No cracking No cracking 


Scaling characteristics 
——Scale Thickness——. 
Exposure temperature OD ID 
1200 F Very slight, not reported 
1350 F Very slight, not reported 
1500 F Very slight, not reported 


Microstructures (etchant: Ovxalic, electrolytically, 1000 D) 


6 months, 1200 F 6 months, 1350F  —«@ months, 1500 F 
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(2) 46,300 50 120 + 
(1) 30,300 72 120 + 
(3) 33,800 57 120 + 
(1) 17,500 57 120 + 
(4) 21,000 57 120 + 
Calculated 
H value 
Unexposed 


Microstructures (etchant: 


JANUARY 1960 


Table 13 Characteristics of TP-321 before and after six-months’ exposure at 1200, 1350, and 1500 F 
Chemical composition 


Mn Si Cr Ni Ti 
Before exposure 0.045 Lt 0.014 0.012 0.47 17.79 12.23 0.45 
After exposure 0.044 1.76 0.011 0.010 0.46 17.80 12.12 0.46 
Grain size after exposure 
Heat-treatment Original 1200 F 1350 F 1500 F 
Solution annealed = Structural 5 5 5-6 5-6 
Short-time tensile, hardness, and impact properties 
Tensile 
Test strength, 0.2% yield, age pe Reduction Brinell Charpy V 
temperature psi psi in 1.4", % of area, % hardness notch, ft-lb 
Room 
temperature (1) 75,100 32,400 29.6 70.4 126 230, 230 
(2)2 81,500 39,900 44.3 80.4 145 130, 132 
79, 500 33 , 600 45.7 79.4 143 116, 115 
(4)2 75,900 29,100 52.9 80.8 137 150, 151 
1200 F (1) 62,000 26, 200 37.1 49 177, 179 
(2) 57,000 34,000 28.6 63.4 155, 151 
1350 F (1) 37,500 18,000 31.4 15.4 153, 145 
(3) 45,300 24,100 32.9 42.4 135, 137 
1500 F (1) 28 ,000 17,000 44.3 32.5 147, 152 
(4) 33,500 22,500 41.4 42 149, 150 


1 Unexposed, (2) six months at 1200 F, (3) six months at 1350 F, (4) six months at 1500 F. 
* Standard 0.357-in. diam X 1.4-in. gage tensile specimens. The others were as shown in Fig. 1. 


Tube dimensions 


Exposure ————Outlet end, in. Max % expansion 
temperature OD 
Unexposed 2.009/2.012 0.517/0.533 uP 
1200 F 2.008/2.013 0.527/0.541 0.04 
1350 F 2.008/2.014 0.517/0.533 0.1 
1500 F 2.007/2.012 0.511/0.537 0 
Flattening characteristics 
Exposure Height at which cracking Calculated H value, 
temperature occurred in. 
Unexposed No cracking 1.63 
1200 F No cracking 1.63 
1350 F No cracking 1.63 
1500 F No cracking 1.63 
Scaling characteristics 
Scale thickness, in. 
Exposure temperature OD ID Type of attack 
1200 F Under 0.002 General 
1350 F Under 0.002 General 
Under 0.002 General 


Aqua regia, X 100) 


6 months, 1350 F 


6 months, 1200 F 
Inside surface 


Inside surface 


6 months, 1500 F 
Inside surface 
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Table 14 Characteristics of TP-347 before and after six-months’ exposure at 1200, 1350, and after 2200 hr 
at 1500 F 


Chemical composition 
C Mn Pp 8 Si Cu Ni Cr Cb Ta 


0.06 1.74 0.022 0.010 0.65 0.17 10.90 17.93 0.56 0.05 


— Grain size ~ 
Before 
Heat-treatment exposure After exposure 
Solution annealed 1200 F 1350 F 1500 F 
9/2 9/2 9/2 8/3 


Short-time tensile and impact properties 
Tensile 
Test strength, 0.2% yield, Elongation Reduction Brinnell Charpy V 

temperature i psi in 1.5", % of area, % hardness notch, ft-lb 
Room 

41,000 

36, 100 

34,300 
1200 F 


1350 F 
1500 F 


(1) Unexposed, (2) 6 months at 1200 F, (3) 6 months at 1350 F, (4) 2200 hr at 1500 F. 


Tube dimensions 
Exposure 
temperature D Yall 
Unexposed 2.004/2.007 0.545/0.571 
1200 F 2.005/2.007 
1350 F 2.005/2 .007 
1500 F 2.045/2.147 2.023/2.027 


* OD at soot blower. 


-————Inlet end, in. ——— Outlet end, in. Max % 
O W OD Wall OD* expansion 


Flattening characteristics 
Height at which cracking occurred, in. Calculated 


Exposure temperature Inlet end Outlet end H value, in. 
Unexposed No cracking No cracking 1.65 
1200 F No cracking No cracking 
1350 F No cracking No cracking 
1500 F 1.75 1.75 


Scaling characteristics 
— Scale thickness—— 
Exposure temperature OD ID Type attack 
1200 F Very slight, not reported 
1350 F Very slight, not reported 
1500 F Very slight, not reported 


Microstructures (etchant: Nitric, electrolytically, x 1000 D) 
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id (2) 44,500 35 120 + 
(1) 35,200 32 37.0 120 + 
(3) 28, 500 66 120 + 
(1) 21,940 39 50.6 120 + 
(4) 16,000 60 120 + 
0.0 
0.0 
7.0 
Unexposed 6 months, 1200 F 6 months, 1350 F 2200 hr, 1500 F 


Table 15 Characteristics of TP 310 before and after six-months’ exposure at 1200, 1350, and 1500 F 


Chemical composition 
C Mn Si Cr Ni 
0.072 1.55 0.007 0.018 0.25 24.56 21.42 


Heat-treatment Structural grain size after exposure 
1200 F 1350 F 


Solution quenched, 2075 F Original 1500 F 
3 (2, 4) 3/6 4/6 4/6 (3) 
Short-time tensile, hardness, and impact properties 
Tensile 
Test strength, 0.2% yield, Elongation Reduction Brinell Charpy V 
temperature psi psi in 1.5”, % of area, % hardness notch, ft-lb 
Room (1) 84,500 32,750 50.0 74.0 128 206, 204, 206 
2) 94,000 45,800 35.0 51.0 163 45, 43, 58 
(3) 92,000 42,400 30.0 36.0 170 18, 17, 17 
(4) 70,300 37,500 13.0 15.0 156 24, 22, 21 ¢ 
1200 F (1) 53,500 15,250 41.0 50.0 167, 176, 182 
(2) 56,000 26,150 30.0 52.0 112, 110, 127 
1350 F (1) 49 ,000 15,800 35.0 44.0 180, 175, 178 
(3) 45,500 , 100 29.0 31.0 45, 38, 43 
1500 F (1) 34,500 14,000 41.0 41.0 162, 170, 165 3 
(4) 45,250 19,800 26.0 28.0 , 26, 


(1) Unexposed, (2) after 6 months at 1200 F, (3) at 1350 F, (4) after 3000 hr at 1500 F. 
2 1300 F rather than 1350 F. 


Tube dimensions 


Exposure Inlet end, in. —-Outlet end, in— Max % 
temperature Wall OD Wall OD* expansion 
Unexposed 1.993/1.998 0.527/0.547 
1200 F 1.992/1.990 0.539/0.549 ie 
1350 F 1.998/2.000 0.522/0.540 0.20 
1500 F 2.046/2.054 0.519/0.538 2.76 
* OD at soot blower. 
Flattening results 
Exposure Height at which cracking occurred, in. Calculated 
temperature Inlet end Outlet end H value, in. 
Unexposed Satisfactory 1.62 
1200 F Satisfactory 1.62 
1350 F Fractured at 1.5 1.63 
1500 F Fractured at 1.75 1.65 
Scaling characteristics 
Scale thickness——— 
Exposure temperature OD Type attack 


1200 F Not reported Assumed very slight 
1350 F Not reported Assumed very slight 
1500 F Not reported Assumed very slight 


Microstructures (1000 D) 


6 months, 1500 F 


Unexposed 6 months, 1200 F 6 months, 1350 F 
Etchant: Ferritic Etchant: Electrolytic 10 per cent chromic acid 
chloride 
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Table 16 Characteristics of 16-25-6 before and after six-months’ exposure at 1200, 1350, and 1500 F 


Chemical composition 
Cc Mn 8 Si Cr Ni Mo 
0.053 1.73 0.020 0.010 0.93 15.70 25.09 5.34 0.18 


Heat-treatment ——— Structural grain size after exposure-————— 
Solution quenched Original 200 F 1350 1500 F 
5/4 (6,3) 5/6 (4) 4/3 (2) 5/6 (4) 


Short-time tensile, hardness, and impact properties 
Tensile 
Test strength, 0.2% yield, Elongation eduction  Brinell Charpy V 
temperature i psi in 1.5",% ofarea,% hardness notch, ft-lb 
Room (1) 67 166 
24 229 
28 229 
( ‘ 197 19, 17, 19 
1200 F 131, 140, 117 


28, 27, 26 
1350 F ‘ 104, 114, 125 
34, 37, 35 
1500 F 103, 80, 88 
51, 50, 52 


(1) Unexposed, (2) 6 months at 1200 F, (3) 6 months at 1350 F, (4) 6 months at 1500 F. 


Tube dimensions 
Exposure ——Inlet end, in——~ ———Outlet end, in. 
temperature OD Wall OD Wall 
Unexposed 1.987/8 
1.995/6 0.503/8 
1.995/8 0.505/6 
1.997/8 0.503/5 


* OD at soot blower. 


Flattening characteristics 


Exposure Height at which cracking occurred, in. Calculated 
temperature Inlet end Outlet end H value, in. 


1.88 
1.81 
3. 


83 
Scaling characteristics 
Scale thickness, in ——. 
Exposure temperature OD 
1200 F 0.0004 0.0015 


1350 F 0.0005 0.0010 
~ 1500 F 0.0010 


Microstructures (etchant: Glyceregia, X 1000 D) 


& & 


6 months, 1200 F 6 months, 1350 F 
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| 
— Max % 
OD? _ expansion 
1.995 0.4 
| 1.995 0.5 
1.998 0.7 
= 
Unexposed 1.34 1.60 
1200 F , 1.84 1.60 
1350 F 1.82 eo 
Type attack 
General 
General 
Unexposed eC 6 months 1500 F 


Table 17 Characteristics of 17-14 Cu-Mo before and after six-months’ exposure at 1200, 1350, and 1500 F 
Chemical composition 
C Mn g 8 Si Cr Ni Mo Cu Cb Ti 
0.128 0.76 0.028 0.014 0.25 15.93 14.54 2.60 2.81 0.53 0.19 
Heat-treatment Structural grain size after exposure-——~—. 
Solution quenched 1200 F 1350 F 1500 F 
8(7) 8(6) 


Original 
7/8 8 


Short-time tensile, hardness, and impact properties 
Test Tensile strength, 0.2% yield, Elongation Reduction  Brinell 

temperature psi psi in 1.5", % of area, % hardness 

Room (1) 67.7 143 
(2) 43.0 197 

51.0 179 

48.5 156 

51.0 

30.3 

53.5 

39.5 


Charpy V 
notch, ft-lb 


120, 112, 116 
53, 58, 63 
95, 102, 89 
64, 62, 64 
58.0 81, 80, 72 
20.0 54, 52, 55 


(1) Unexposed, (2) 6 months at 1200 F, (3) 6 months at 1350 F, (4) 6 months at 1500 F. 


1200 F 
1350 F 
1500 F 


Tube dimensions 
exposure ——-——Inlet end, in. 
temperature 


Unexposed 
1200 F 


Max % 


Outlet end, in.——--—. 
expansion 


OD Wall OD* 


2.003/2.005 
2.003/2.005 
1350 F 2.008/2.011 
1500 F 2.013/2.015 


+ OD at soot blower, 


0.502/0.507 
0.522/0.552 
0.529/0.542 
0.529/0.536 


2.002/2.004 0.532/0.542 2.004 
2.005/2.009 0.520/0.552 2.007 
2.003/2.012 0.525/0.540 2.012 


Flattening characteristics 


Exposure 
temperature 


Calculated 


H value, in. 


1.64 


Height at which cracking occurred, in. 
Inlet end Outlet end 


Unexposed 
1200 F 
1350 F 
1500 F 


No cracking 
1.26 

No cracking 

No cracking 


1.61 
1.64 
1.66 


1.30 
No cracking 
1.43 


Scaling characteristics 


Exposure temperature Type attack 


Scale thickness, in. 
OD ID 
General 

General, pits 0.002 in. deep 
General, pits 0.002 in. deep 


0.0002 
0.0003 
0.0003 


0.0020 
0.0020 
0.0020 


1200 F 
1350 F 
1500 F 


Microstructures (etchant: Glyceregia, 1000 D) 


6 months, 1 
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46, 45, 43 
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Table 18 Characteristics of Croloy 15-15N before and after six-months’ exposure af 1200, 1350, and 1500 F 
Chemical composition 
C Mn 8 Po a Cr Ni Mo Cb-Ta W 


Before exposure 0.12 1.27 0.013 0.010 0.89 14.98 15.42 1.45 1.22 1.24 
After exposure 0.12 1.26 0.016 0.010 0.87 15.14 15.40 1.45 1.12 1.27 


-Grain size after exposure— ~ 
Heat-treatment Original 1200 F 1300 F 1500 F 
Solution annealed Structural 6 6-7 6-7 5-6 


Short-time tensile, hardness, and impact properties 
Tensile 0.2% a Reduc- Charpy V 
Test strength, yield, ‘ion i tion of __ Brinell notch, 
temperature i i ; area, % hardness ft-lb 
Room 
temperature (1) 148 137, 134 
(2) 161 65, 63 
175 
155 


1200 F 
1350 F 
1500 F 


(1) Unexposed, (2) 6 months at 1200 F, (3) 6 months at 1350 F, (4) 6 months at 1500 F. 
* Standard 0.357 in. X 1.4-in. gage tensile specimens. All others as shown in Fig. 1. 


Tube dimensions 
Exposure ———— Outlet end, in. 
temperature OD ID Max % expansion 
Unexposed 2.000/2.007 0.525/0.548 
1200 F 2.000/2.003 0.531/0.553 None (0.19% reduction in size) 
1350 F 2.000/2.002 0.521/0.556 None (0.25% reduction in size) 
1500 F 1.983/1.985 0.513/0.552 None (1.1% reduction in size) 


Flattening characteristics 
Exposure Height at which cracking 
temperature occurred, in. Calculated H value, in. 
Unexposed No cracking 1.63 
1200 F No cracking 1 
1350 F 1 
1 


64 
53 
1500 F No cracking 


.52 


Scaling characteristics 
Exposure ——Sceale thickness, 
temperature OD ID Type attack 
1200 F Under 0.002 
1350 F Under 0.002 0.005 
1500 F Under 0.002 0.006 General with pronounced exfoliation 


Microstructures (etchant: Aqua regia, X 100) 


6 months, 1200 F 6 months, 1350 F 6 months, 1500 F 
Inside surface Inside surface Tnside surface 
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(1) 50, 500 18,500 34 49.5 122, 127 
(2) 64,000 36,500 34 27.9 66, 61 
’ (1) 47 ,000 25,000 35 35.1 113, 108 
(3) 50 , 000 30,000 41 65.1 46, 47 
(1) 36, 500 23 , 500 61 55.1 120, 113 
(4) 37, 200 26, 200 71 60.7 70, 74 
15650 / 65 


Table 19 Characteristics of Inconel alloy before and after six-months’ exposure at 1200, 1350, and 1550 F 
Chemical composition 


Cc Mn P s Si Cr Ni 
0.04 0.24 ... 0.007 0.26 15.70 75.57 
4 Heat-treatment 


—————Structural grain size after 
1350 


Solution quenched, 2050 F Original 1200 F 1500 F 
1/4 (6) 1/3 (4, 5) 1/3 (4, 5) 1/4 (5) 
Short-time tensile, hardness, and impact properties 
Tensile 

Test strength, 0.2% yield, Elongation Reduction Brinell Charpy V 

temperature psi psi in 1.5", % ofarea,% hardness notch, ft-lb 
Room (1) 84,500 24,500 56.0 66.0 128 207, 202, 218 

(2) 89,400 28 ,600 48.0 65.0 132 144, 146, 141 
(3) 88,700 27,100 49.0 64.0 126 184, 217, 200 
(4) 86,000 26,200 50.0 68.0 121 184, 192, 187 
1200 F (1) 64,900 15,800 56.0 63.0 174, 171, 170 
(2) 59, 500 17,350 44.0 56.0 181, 165, 162 
a 1350 F (1)? 61,000 15,500 54.0 58.0 168, 184, 172 
: (3) 61,350 15,400 42.0 47.0 202, 175, 219 
‘4 1500 F (1) 42,750 14,500 67.0 74.0 167, 174, 170 
4 (4) 51,500 14,500 45.0 42.0 167, 170, 162 


(1) Unexposed, (2) after 6 months at 1200 F, (3) at 1350 F, (4) at 1500 F. 
* 1300 F rather than 1350 F. 


Tube dimensions 


Exposure ———Inlet end, in-—-——. —Outlet end, in.—~ Max % 
‘* temperature OD Wall OD Wall OD* expansion 
Unexposed 2.004/2.006 0.498/0.507 
1200 F 2.002/2.004 0.488/0.517 
1350 F 2.000/2 .002 0.489/0.523 bic 
1500 F 2.033 /2.033 0.486/0.503 1.40 


* OD at soot blower. 


Flattening results 


Exposure Height at which cracking occurred, in. Calculated 
temperature Inlet end Outlet end H value, in. 
Unexposed Satisfactory 1.60 
1200 F Satisfactory 1.60 
1350 F Satisfactory 1.60 
1500 F Satisfactory 1.62 
Scaling characteristics 
Exposure temperature Type attack 


Scale thickness, in.— 
OD ID 


1200 F Not reported Assumed very slight 
1350 F Not reported Assumed very slight 
1500 F Not reported Assumed very slight 


Microstructures (1000 D) 


6 months, 1200 F 6 months, 1350 F 
Etchant: Glyceregia 


Unexposed 6 months, 1500 F 
Etchant: Ferric 


chloride 
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Table 20 Characteristics of Incoloy alloy before and after six-months’ exposure at 1200, 1350, and 1500 F 
Chemical composition 


Cc Mn P Ss Si Cr Ni 
0.05 0.79 0.007 0.55 20.69 34.66 
Heat-treatment —————Structural grain size after exposure——————~ 
Solution quenched, 2050 F Original 1200 F 1350 F 1500 F 
Duplexed 3/8 2/5, 6/8 2/5 (6, 7) 2/6 
Short-time tensile, hardness, and impact properties 
Tensile 0.2% Elonga- Reduc- 

Test strength, yield, tion in tionof __Brinell Charpy V 
temperature psi psi 1.5”,% area, % hardness notch, ft-lb 
Room (1) 82,500 28,500 51.0 68.0 112 206, 202, 214 

(2) 111,350 57,000 28.0 36.0 190 38, 32, 32 
(3) 96, 400 40,000 35.0 45.0 154 49, 50, 59 
(4) 84,900 34,300 45.0 66.0 130 129, 112, 116 
1200 F (1) 60, 500 16,500 49.0 69.0 183, 187, 174 
(2) 83,750 50,800 20.0 34.0 57, 56, 61 
1350 F (1) 54,900 16,500 48.0 66.0 169, 181, 178 
(3) 68,750 28,000 30.0 48.0 91, 179, 98 
1500 F (1) 36, 500 15,500 57.0 68.0 180, 173, 174 
(4) 49, 200 18,900 44.0 48.0 110, 176, 158 
(1) Unexposed, (2) after 6 months at 1200 F, (3) at 1350 F, (4) at 1500 F. 
@ 1300 F rather than 1350 F. 
Tube dimensions 
Exposure Inlet end, in————~ ——Outlet end, in— Max % 
temperature OD Wall OD Wall OD* expansion 
Unexposed 0.001/2.012 0.493/0 .503 ais 
1200 F 2.002/2 .003 0.499/0.507 
1350 F 2.003/2 .004 0.495/0.510 
1500 F 2.056/2.061 0.465/0.496 3.00 
* OD at soot blower. 
Flattening results 
Exposure Height at aes cracking occurred, in. Calculated 
temperature nlet end Outlet end H value, in. 
Unexposed Satisfactory 1.60 
1200 F ID cracking at 1.81 1.60 
1350 F Satisfactory 1.60 
1500 F Satisfactory 1.62 
Scaling characteristics 
—Scale thickness, in.— 
Exposure temperature OD ID Type attack 
1200 F Not reported Assumed very slight 
1350 F Not reported Assumed very slight 
1500 F Not reported Assumed very slight 

Microstructures (1000 D) 

Unexposed 6 metie, 1200 F 6 went 1350 F 6 months, 1500 F 
Etchant: Ferric Etchant: Glyceregia 
chloride 
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Blast-Furnace-Gas Cleaning for 
lron and Ferromanganese Production 


A process for removing the superfine smokes from ferromanganese blast-furnace gases 
having low investment and low operating cost is described. This method solves both the 
air and stream-pollution problems. The ferro blast-furnace gas is cleaned from 10 to 15 
grains per cu ft to less than 0.01 grains per cu ft by a wet-washing procedure in a Pease 
Anthony venturi scrubber. The alkali cyanide stream-pollution problem is solved by 
the control of pH and complete recycle of the liquors so that the cyanide is stripped by 
the incoming gas and is subsequently burnt in the stoves or boilerhouse. The liquor ts 
separated from solids in a Dorr thickener. The solids are removed in the underflow 
as a thick siurry which is either returned to the furnace or dumped in a pit in the ore 


J. S. MACKAY 
J. F. KEARNS 


Pittsburgh Coke & Chemical Company, 
Pittsburgh, Pa. 


yard. The overflow 1s pumped back to the venturi. 


The chemistry and physical equip- 


ment of the process are described. 


is REMOVAL of the dense-fine fume which occurs in 
the production of ferromanganese in the blast furnace is much dif- 
ferent from that in an iron blast furnace. C. A. Bishop, et al., 
have very adequately described the problem and the make-up of 
the fume. As they have stated, the solids present in the gas are 
80 per cent between 0.1 and 1.0 micron. At one time, this was 
discharged to the atmosphere as a dense fume. This was not only 
obnoxious but was economically wasteful since the fuel value of 
the gas was not used. This fume is largely condensed vapors 
which were volatilized in the high-temperature part of the furnace 
and condensed in the cooler zone as a smoke. One of the groups 
of materials not made in quantity in the blast-furnace operation 
but which is produced in the ferromanganese furnace due to a 
combination of higher temperature and different chemistry is the 
alkali-metal cyanides. This poses the second problem. Wet- 
washing procedures are more efficient in removing fine fumes but 
cyanide wastes cannot be disposed of in the Ohio River. The 
solution worked out by Bishop, et al.,! was electrostatic precipita- 
tion with subsequent burning of the dust and briqueting. This 
process involves large investment and high maintenance and is 
considered less efficient in removing fine fumes. For these 
reasons and others, Pittsburgh Coke & Chemical decided to solve 
‘the problem, if possible, by a wet-wash system. This has been 
successfully done and will be described. If wet washing is used, 
the problem becomes one of river pollution more than air pollu- 
tion. The river-pollution problem was solved in order to solve 
the air-pollution problem. 

In addition to the method used at the Duquesne Works of 
United States Steel,? it had been shown that a wet-washing sys- 
tem combining a gas washer and a Thiesen disintegrator gave 
satisfactory gas under the air-pollution standards. The cyanide 
problem was not solved and the solid slurry was difficult to 
handle due to its tendency to cement. This latter quality is due 
to the fineness of the material as well as the chemical composition, 


1 Bishop, Brisse, Thompson, Liebel, and Swaney, ‘‘Cleaning Ferro- 
manganese Blast Furnace Gas,’’ Yearbook of the American Iron and 
Steel Institute, 1951, pp. 459-475. 

2 C. H. Good, ‘‘Ferromanganese Gas Cleaning at Duquesne Works,” 
Blast Furnace, Coke Oven, and Raw Material 1954 Proceedings, 
AIME, pp. 169-178. 

Presented at the Pittsburgh Bicentennial Conference of the ASME 
Pittsburgh Section on Pollution Problems and Their Solution, April 
20-21, 1959. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, May 8, 
1959. 
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Not only is the plugging problem existent in the aqueous solutions 
but also in the lines containing the wet gas. The Thiesen disinte- 
grator rotating bars would be coated with a very hard cement- 
like material. For this reason two disintegrators were deemed 
necessary. 

The process involves very low investment, has very low main- 
tenance charges, and has been eminently successful in solving 
both the air and stream-pollution problems. In brief, the process 
involves passing the gases from the normal dust collector through 
a Pease Anthony venturi scrubber to a gas main feeding to either 
the stoves or boiler house. The wash water of the venturi serub- 
ber goes to the Dorr thickener where the solids are discharged as a 
sludge either to the skips or to the ore yard and the overflow ef- 
fluent is recycled to the venturi. This will be described in more 
detail later. 


100 
90 
80 
70) 


Fig. 1 HCN has a vapor pressure and is removable; M* CN~ = HCN + 
MOH 


A study of the two methods used for cleaning the gas showed 
that the advantages were certainly on the side of wet-scrubbing if 
the plugging and cyanide-pollution problems could be solved. 
The alkali cyanides are formed in the furnace by fixation of nitro- 
gen from the air in the hot zone. These are condensed in the 
cooler part of the furnace to a smoke and essentially are nonvola- 
tile from there on. When they are wet by water, they dissolve 
and partially decompose to HCN if the pH is at the proper level. 
The more acidic the solution, the more the decomposition. Fig. 1 
shows the degree of decomposition with pH. Pure HCN boils 
at 27 C, but is readily soluble in water. The vapor pressure is suf- 
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HCN CONCENTRATION - LIQUID (PPM) 
Fig. 2 Equilibrium data for HCN-water-air system at 20, 30, 40, and 50 C 


ficiently high to be removed by aeration, preferably with acidic 
gases so that pH is kept at the lower levels. The solids in the gas 
are largely alkaline and one would presume the liquor to be highly 
alkaline. However, there is considerably more carbon dioxide 
than total alkali so if proper absorption is allowed the pH could 
be acidic enough for HCN stripping. This could not be done 
simultaneously with solution of the cyanide but if the liquor 
could be cycled back to the scrubbers a large part of the HCN 
would be removed in the gases. The vapor pressures of various 
concentrations of HCN in water are shown in Fig. 2. It was 
known that recycle of these liquors was nearly impossible and 
that even on iron-blast-furnace-gas wash water, attempts at re- 
cycle had been unsuccessful. On examining the analysis of river 
water and the time available for treatment of the water in the 
Dorr thickener, it was felt that the Dorr thickener effluent could 
be prepared to be at least as good as that obtained from the river. 
On investigation of the reasons for failure of the recycle on blast- 
furnace water, it was found that at least in one case the liquor 
was pumped to a spray-pond cooling system and then back to the 
gas washer. The plugging occurred at the elbows on the return 
line. This plugging could well have been due to supersaturation 
in the air-cooling pond. As a matter of fact, it was found possible 
to recirculate when the sprays were shut off. At the pH desired 
for cyanide removal from the liquor, calcium-carbonate solubility 
is roughly at its minimum. One of the very obvious ways in 
which cement formation could start in a wet system was by in- 
creasing the calcium solubility either with more acid or more 
alkaline conditions and then reducing the solubility by changing 
the pH. In many cases, this would happen at the points where 
plugging could be the worst problem. If the water was com- 
pletely recycled some cooling advantage obtained in the use of 
fresh water would be lost. As a matter of fact, it is obvious that 
almost all cooling would have to be by evaporation. Calculation 
shows that this was the major method of cooling anyway and an 
equilibrium liquor temperature should be arrived at around 150 F. 

The first proposed system consisted of feeding the gas into a 
stripped-gas washing tower and spraying with approximately 
2500 gpm of water. Sprays were also installed at every piping 
bend. After passing throug) the washing tower, gas would then be 
fed into a rotary Thiesen disintegrator where it would be sprayed 
again with about 250 gpm of water. The gas leaving the disinte- 
grators was expected to be exceptionally clean. The disintegrator 
would clean the gas no matter what the dust loading, whereas an 
electrostatic precipitator efficiency is based upon the dust loading 
fed into it. The water discharging from the cooling tower and 
the disintegrator goes to the Dorr thickener. It was planned 
to remove the cyanide before discharging the liquor to the river 
by chlorination. The Wallace & Tiernan method of continuous 
chlorination in which the cyanide is oxidized to cyanate was to be 
used. This chlorination might be before the Dorr thickener so 
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that it could be used for the retention tank, or on the effluent 
after the thickener, in which case a retention tank had to be 
added. The State Health authorities of Pennsylvania will not 
allow chlorination only to cyanate which is the normal procedure 
elsewhere. They require full chlorination to Nz and CO,. Fur- 
thermore, they would not accept the continuous-chlorination 
process of Wallace & Tiernan but insisted on batch chlorination in 
which all the liquor of one lake is tested free before discharge to 
the river. They would consider proof of a continuous process but 
required at least a long run on a pilot scale before acceptance. 
This would automatically require a batch system anyway. 
Roughly, the cost rose from $6 to $15 per ton ferromanganese and 
it was decided that either the water would have to be partially 
recycled so that only part would be chlorinated, or it would be 
necessary to operate on complete recycle. 
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DUST LOADING—GRAINS/C.F (N.T.P) 
Fig. Summary of blast-furnace venturi data 


It was decided to operate on complete recycle. At the same 
time, another change was made. It appeared for several reasons 
that the use of a Pease Anthony venturi scrubber in place of 
water sprays in the gas washer would be desirable. There was 
plenty of pressure for the specified pressure drop to give a 
guaranteed cleanliness of 0.03 grains per cu ft on basic iron and 
0.08 grains per cu ft on ferromanganese. This required only 600 
gpm of water as compared to 2500 in the old system. Fig. 3 
shows the guarantee of solid content of the gas effluent from the 
Pease Anthony scrubber. At the same time, it was decided to 
combine the two furnace-gas and liquor systems from A and B 
furnaces so that a venturi was to be put on both furnace gases 
with both liquors going to the same Dorr thickener, and the 
effluent recycled back to each venturi. The gas was to go to a 
common main and through the Thiesen disintegrators to the 
stoves or directly to the boilerhouse. In the previous system, 
a new Dorr thickener had to be supplied to handle the amount of 
liquor for the two furnaces. There are two advantages to the re- 
vised system. One, the gas from the ferromanganese unit would 
be diluted before going through the disintegrator. And second, 
the cyanide would be stripped with two furnace gases rather 
than one. The one disadvantage was that the cyanide content 
would be considerably higher due to the lower amount of water. 
Again, of course, dilution from the other furnace would be helpful. 
Of course, a shut down on both furnaces was being risked if the 
system did not work. The amount of water required for the 
sprays in the disintegrator was decreased from 250 to 150 gpm. 
The old gas washers or coolers were retained as knockout cham- 
bers for the gas after the venturi. The venturi was modified as 
shown in Fig. 4 so as to better handle any tendency for solid 
plugging. 

The whole system as finally designed is shown in Fig. 5. Water 
sprays to both of the gas washers had been included so that if 
necessary solution could be pumped to there as well to the venturi 
to get CO, absorption for pH control. Obviously, while solid 
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pickup is excellent in a venturi scrubber, gas absorption would 
be very small because of contact time. Gas stripping from a 
solution did not require this contact time, although again pump- 
ing Dorr effluent to the gas washers would be helpful in lowering 
the cyanide content somewhat more. 

It is this design, shown in Fig. 5, that has been used success- 
fully. The gas at the top of the blast furnace is cooled from 900+ 
to 600 to 700 F by water sprays. This gas passes from the furnace 
to the dust catcher and from the dust catcher into the venturi 
scrubber where it is sprayed with 600 gpm of recycle water. The 
pressure drop was from 30 to 70 in. of water and the velocity 
through the throat of the venturi ranged from 250 to 400 fps. This 
is the reason for the exceptionally good cleaning obtained. 


Conventional design 


,A BLAST FURNACE 

,B' DUST COLLECTOR 

,C’ VENTURI 

,D' COOLING TOWER 
MOISTURE ELIM. 
DISINTEGRATOR 
THICKENER 


—=GAS 
—-—WATER 
SLUDGE 


This gas then goes to the old gas washer or cooling tower with the 
hurdles and packing removed. As mentioned, water may or may 
not be added in sprays in this gas washer. Both furnaces have 
identical arrangement. The gas goes into a common main and a 
portion of it is channeled through the disintegrators for final 
cleaning before going into the blast-furnace stoves. The gas dis- 
charging from the disintegrators ran from 0.0005 to 0.0001 grains 
of solid per cu ft. In fact it was hardly measurable. The scrub- 
bing water from both furnaces feeds to the Dorr thickener. The 
overflow water from the Dorr thickener flows to a sump which is 
maintained at constant level with make-up water and from here 
is pumped back to each venturi or can also be pumped to the gas 
washers and to the disintegrators. Make-up water is required 


P. C. & C. design 


Fig. 4 Venturi redesign 


10 goiLer® 


Fig. 5 Gas cleaning system 
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to replace the evaporation losses and that removed in the Dorr 
underfiow concentrate. The pump capacity at the bottom of the 
thickener was higher than desired to maintain a reasonably thick 
underflow if discharged directly so part of this was recycled to 
the top feed of the Dorr. A bypass from this fed to either the 
skips of both furnaces or into a pit in the ore yard. The only 
change contemplated in this system is a filter for the Dorr thick- 
ener sludge. 

This system was installed in the summer of 1957 and went into 
operation in the fall of 1957 with both furnaces on basic iron. No 
problems were experienced except that there was some difficulty 
in initial periods in getting adequate settling in the Dorr thick- 
ener. The pH steadied out to 9 to 9.5 and the gas cleaning was 
excellent. This was an excellent time to run a contro] experiment 
to determine the efficiency of the cyanide stripping in order to 
know what the hazards could be around the Dorr thickener and 
at the various sumps and how efficient the cyanide stripping 
could be. One ton of sodium cyanide was pumped into the A 
flume as an aqueous solution such that the HCN equivalent was 
about 3 lb per min. The effluent from the Dorr thickener (ca- 
pacity of 300,000 gal) rose to 220 ppm of HCN in 7 hr. The at- 
mosphere was tested at all likely points. Concentration was less 
than 15 ppm. Industrial atmosphere limit is 20 ppm. High con- 
centrations existed in closed sumps. The cyanide was rapidly 
removed at the venturi, the total content dropped rapidly after 
addition of sodium cyanide stopped and was down to 50 ppm 
in 5 hr. The maximum removal rate was 1.5 lb HCN per min, 
whereas the addition was 3 lb per min so that maximum 
equilibrium concentration had not been reached. It was found 
that atmospheric dilution of the cyanide in the atmosphere was 
very rapid. For example, on opening a closed sump containing 
the effluent from the Dorr thickener, the cyanide concentration 
at foot level was 60 ppm. This became 15 at waist level and 5 
ppm at head level when standing in a direct steam flow rising 
from the sump. The Dorr thickener had been surrounded with 
an 8-ft wall of transite siding. The HCN concentration inside this 
wall on the bridge over the lake was approximately 3 ppm and 
was essentially undetectable at the roof of the wall. This ex- 
periment provided confidence in going ahead with the ferroman- 
ganese run so far as personnel danger and cyanide stripping was 
concerned. 

Also during this 3-month period of operation in which the two 
furnaces were on complete recycle for several-week periods, no 
plugging of the return lines, the sprays, nor the throats of the 
venturi scrubbers was noticed. 

Operation started on spiegeleisen which is very similar to 
that of ferromanganese and also gives high cyanides. At this 
time, B furnace was also operating. Cyanides rose very rapidly 
in the Dorr thickener and leveled out at about 200 ppm of cy- 
anide. Toxicity measurements in the general areas fitted very 
closely with the earlier experiment described. No trouble was 
noticed in the recycle and the pressure drop across the ven- 
turi showed no problem. The gas was cleaned out of the venturi 
scrubber to less than 0.01 grains per cu ft. At this time, the dis- 
integrator was used for the gas going to the stoves. Ferroman- 
ganese was started and again no problems resulted, but the 
cyanide in the Dorr thickener did rise to from 250 to 400 ppm. 
There was a slow increase in pressure drop across the venturi on 
A furnace which was on ferromanganese. There was none, how- 
ever, at B furnace showing that the build-up was not due to the 
recycle water since this was identical for each venturi, but rather 
to the dust-laden gas hitting the sprays. There was no particular 
difficulty in cleaning out the venturi when the pressure drop rose 
above 70 in. of water. The results on the ferromanganese gas are 
shown in Fig. 6. The soluble solids rose to 1.8 per cent. These 
salts were 94 per cent K and Na carbonates and chlorides. The 
solids content of the Dorr overflow kept in good control after the 
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Fig.6 Summary of blast-furnace venturi data 


start-up when again it was high. The pH maintained itself in a 
9.0 to 10.5 range so that cyanide stripping was maintained 
throughout without having to pump water to the gas washers. 
The system has been used continuously for the past year and a half 
though only part-time on ferromanganese. 

Prior to the installation of the present system, the blast-fur- 
nace stoves had to be cleaned out every six months. In the total 
of the past year and a half, a stove has not been cleaned and 
small checker openings have been installed to get some 30 per 
cent more efficient hot-blast heat for the furnace operations and 
to allow higher temperatures resulting in substantial savings in 
coke and production costs. The boilerhouse blast-furnace gas 
burners formerly required monthly maintenance cleaning. These 
burners have not been cleaned for over a year and a half. All 
maintenance charges have certainly been at a minimum. 

The heat balance came out as predicted and the water out of 
the venturi is at 150 F, the gas slightly above this, and roughly 
15-F cooling is obtained in the cycle in the Dorr thickener. The 
water fed back to the venturi is at 135 F. There is one disad- 
vantage in that rather wet gas has to be taken into the gas mains. 
Actually there is no trouble in handling this gas either for stove 
underfiring or for boilers. A fair amount of condensation is ob- 
tained in the drip legs in the lines to the burners. This water 
contains some HCN since it is in equilibrium with the gas and 
the water can be recirculated to the Dorr thickener. The disin- 
tegrator worked equally well on recycle water or Ohio River 
water. As a matter of fact, it worked better on recycle water be- 
cause it was unnecessary to filter. Actually the use of disintegra- 
tors has been stopped and again there was no trouble with the 
ferromanganese gas. During the two periods of eight-weeks 
operation there was no plugging in recycle lines. There was never 
any trouble on the B-furnace system in the venturi. Some 
build-up did occur in the venturi on A furnace. We are not sure 
that this might not have been due to a few periods when recycle 
water was lost for short periods. 

In the second long run, it was decided to operate on one furnace 
alone. This was partly necessitated by the shutdown of the other 
furnace. In this case, the pH of the overflow water gradually 
rose to 12 and the cyanide level rose to above 1500 ppm. While 
this would be a very dangerous level of cyanide at lower pH’s it 
caused no problems at all in this case. The reason is explained in 
Fig. 1 which shows that the potassium cyanide was not decom- 
posed and there would be essentially no free HCN in the solution. 
Thus the vapor pressure would be essentially zero. However, it 
was not desired to continue this and the alternative mentioned 
earlier was adopted. Water was pumped to the gas washer so 
that carbon-dioxide absorption would occur. In a matter of a 
few hours, the pH droppéd to 10.5 and the cyanide dropped to its 
original level of 600 ppm. 

The material balance for the one furnace system is shown in 
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Table 1 Material-flow single-furnace operations, furnace A, in Ib per min 


Dorr Thickener 


— Venturi gas — ———- Water ——-~. ———discharge———. 


Out 


Wate 
Total solids. . 
CN total. 


soluble. . 


Table 1. Unfortunately the cyanide content was not obtained 
either in the entering or outlet gas so the production rate is not 
available at this time. This was run on hourly composites over 
an 8-hr period. The soluble salts had been at equilibrium for 
some time. The balance is not very good in places but it does 
give a reasonable picture of what is going on. 

In a balance run on the two-furnace system, the gas from A 
furnace contained 60 lb of dust per min and B furnace 20 lb per 
min. Solution to B venturi had 250 ppm of cyanide and came out 
at 100 ppm. The drip legs have contained from 60 to 80 ppm of 
cyanide showing that good equilibrium is attained. 

In the two-furnace system, water never had to be pumped to the 


gas washers. All the data will be gathered together and full de- 


tails published on solid flows, rates, material balances, solids 
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In Out Underflow 


build-up, and so forth, in a later publication. This would be out 
of place in a discussion of air pollution. 

Thus it is known that the process can operate with very low at- 
mospheric pollution, no river pollution, essentially no personnel 
hazard, and very low maintenance. It is also known that this 
can be done with the lowest investment commensurate with any 
gas-treatment system. It is not known what problems would re- 
sult from having gas coolers so that variable temperatures of a 
liquid system would exist. It is not known how long the process 
could run without a cleanout. Certainly, observation has shown 
that it could be considerably over the eight-weeks period and there 
is no point at which it would be difficult to remove even the hard- 
est solids. 

United States Patent No. 2,877,086 has been issued covering 
the process as described. 
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Influence Coefficients for 
Open-Crown Hemispheres 


G. D. GALLETLY 


Mechanical and Electrical Engineering 
Department, Emeryville Research Center, 
Shell Development Compeny, Emeryville, 

Calif. Assoc. Mem. ASME 


Stresses and deflections of hemispherical shells with nozzle openings lying on the axis 
of symmetry may be determined accurately and rapidly with the aid of the tables given 
in the paper. The tables give the influence coefficients (i.e., stresses and deflections ob- 
tained with unit applied loads) for the edge and interior points for various types of edge 


bending load. The radius-ihickness ratios of the hemispherical shells considered varied 
from 15 to 250; the total central angle of the opening varied from twenty to one hun- 


dred degrees. 


An illustration of the utility of the coefficients is given by considering a hemisphere- 
cylinder intersection problem. 


Introduction 


— Sa SHELLS with crown openings (see 
Fig. 1) are frequently used in a number of industries as portions 
of a complete pressure vessel. In order to design these vessels 
properly at their point of attachment with the cylindrical nozzle 
it is necessary to undertake a bending analysis of both the open- 
crown hemisphere and the cylinder. The reason for this is, of 
course, because the pressure loading will cause the two parts of the 
vessel to undergo different elastic displacements (and rotations 
if the pressure is not uniform). Since both cylinder and open- 
crown hemisphere must fit together eventually, this means that 
axisymmetric edge moments and edge transverse shears will be 
produced at their common junction. These edge moments and 
shears then cause bending in the cylinder and hemisphere. The 
analysis of the bending of the cylinder is well-known and offers 
little difficulty (Timoshenko [8]).!. The bending analysis of the 
open-crown hemisphere, however, is not so easy. 

In order to determine the unknown edge moments and shears 
acting at a junction the designer would like to have available the 
edge influence coefficients (by edge influence coefficient is meant 
the displacement or rotation at an edge due to a unit moment or 
shear applied at the edge or the far edge). If these are available 
for the various shells meeting at a junction, then the formulation 
of the compatibility equations is relatively simple. Solving the 
compatibility equations then yields the unknown edge moments 
and shears. The final stress distributions can then be quickly 
determined if the interior influence coefficients are also available. 

The problem of determining the edge influence coefficients for 
hemispheres with crown openings has been investigated before by 
Galletly [3]. In that work he selected a specific hemispherical 
shell (radius-thickness ratio of 90.5, total central opening of 21 
deg) and compared the results obtained by different methods for 
calculating the edge influence coefficients. The criterion of ac- 
curacy used in that investigation was a solution due to Love [6] 
which involves associated Legendre functions of complex degree. 
It was found that an approximate Geckeler-type solution gave 
errors in the edge influence coefficients of about 15 per cent; the 
use of the approximation cot @ ~ 1/¢ (which yields a solution in 
terms of Kelvin functions of zero order) was better in that it gave 
errors of only 5 per cent. 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Power Division and presented at the Semi- 
Annual Meeting, St. Louis, Mo., June 14-18, 1959, of Tar American 
or MEecHANICAL ENGINEERS. 

Notre: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, February 
5, 1959. Paper No. 59—SA-4. 
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Fig. 1 Open-crown hemisphere used as a pressure vessel head 

While the work cited is of value, it is of limited application since 
only one shell geometry and size of opening was investigated. 
Also, interior influence coefficients were not considered. It would, 
of course, be possible to consider various shell geometries and de- 
termine the ranges over which the cot ¢ ~ 1/@ approximation 
gave good results. However, it was felt to be more profitable to 
calculate the influence coefficients by numerical integration of the 
differential equations of equilibrium. In this way no approxima- 
tions need be made and the influence coefficients are immediately 
available without having to resort to consulting other tabulations 
to calculate them. 

Purpose and Scope of the Present Work. The aim of the present 
investigation was to determine the edge and interior influence co- 
efficients for open-crown hemispheres subjected to edge bending 
moments and edge shearing forces. The coefficients were obtained 
by numerically integrating the classical small-deflection differen- 
tial equations of equilibrium (given in the Appendix) on a high- 
speed digital computer. The ranges of the parameters covered 
were as follows: 

r/h = 15, 25, 50, 75, 100, 125, 150, 200, 250 
gd: = 50, 40, 30, 20, 10 deg 
where 
r = radius of hemisphere 
h = thickness of hemisphere 

¢: = semiangle subtended by the opening (see Fig. 2) 

The governing differential equations were rendered dimension- 
less as were the edge bending loads and the influence coefficients. 
Thus the results given apply to many more shells than just those 
investigated. 
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Methods Used to Determine the Influence Coefficients 


From equations (37), (38), and (39) in the 4 ppendix it can be 
shown that the meridional bending moment My and the hori- 


zontal shearing force Hg (see Fig. 2) may be expressed as follows: 


1 


Me = — 
{Alvf(2:) — + Blof(n) + 
+ Clvf(a2) — + + (1) 
1 
rH, = aan [Am + By, + Cre + Dy] (2) 
where 


A, B, C, D = constants chosen to satisfy the boundary 
conditions 
vy = Poisson’s ratio 
wir = [12(1 — — (3) 
= co-latitude angle 
Z:....¥Y¥: = linearly independent solutions of the 
homogeneous differential equations 
fl...) = d/d&...) + v cot o(...) 


The boundary conditions of the problem are shown in Fig. 2. 
Using these boundary conditions in conjunction with the equa- 
tions given one obtains: 

7H = + a2B + + 
My = anA + axB + axC + anD 
= anA + + + ayD 
= + + aC + ayD 


(4) 


where the a;; represent numerical values. The system of equa- 
tions (4) can be written more conveniently in the matrix form 


{x} = [al{y} (5) 


where the vector z represents the edge moments and shears, y the 
constants A ....D, and the square matrix a represents the a;;. 


Solving equation (5) in the usual way one finds for the constants 


{y} = la} {z} (6) 


Equation (6) gives the constants A .... D in terms of the edge 
moments and shears acting at both edges, i.e., @ = 90 deg and 
¢ = dx. To determine the influence coefficients for any particu- 
lar loading, say a unit moment at @ = @», the following procedure 
is adopted: 

(i) Mo, Ho, and Hg» are set equal to zero in equation (6). This 
then determines the values of the constants A .... D which are 
consistent with the applied loading of Mg: = 1. 

(ii) The results of step (i) are then used in conjunction with 


Fig. 2 Boundary conditions of problem 
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equations (37) to (40) and the numerical integration results. This 
then determines the required influence coefficients. 


The coefficients for the other edge loadings are found in a similar 
manner. 


Edge Loadings Considered and the Quantities Tabulated 


In order to avoid using a specific value of r the quantities 
tabulated for various values of r/h were: 


Er?V, Eré, rNe, and Mg 
where 
V = angle of rotation of a tangent to a meridian (+ when | 
the positive semitangent turns inward) 
Ne = circumferential direct stress resultant (++ when tensile) 
Mg = meridional bending stress resultant (+ when producing 
tension on inner surface) 
6 = horizontal displacement (+ outward) 
r = radius of hemisphere 
Young’s modulus 


& 


The other stress resultants Ng and Mg are usually not the 
controlling ones and, hence, have not been given. 

A not unexpected result was also found for the influence co- 
efficients. This was that the influence coefficients at the edge 
away from the edge where the loading was applied were insig- 
nificant. That is to say, an edge effect existed. Even in the 
worst case (r/h = 15, d2 = 50 deg) the influence coefficients at 
the distant edge were less than 10 per cent of those where the 
loading was applied. Thus, in the tabulation of the influence co- 
efficients (Tables 1 to 3), it was assumed that there was no interac- 
tion between the two edges. 

The tables list first the influence coefficients at @ = 90, 85, and 
80 deg due to edge loadings at @ = 90 deg. The two edge loadings 
considered are shown in Fig. 3(a). The first lines in the influence 
coefficient tabulation correspond to the rHg-loading, the second 
lines to the M g-loading. 

The coefficients due to edge loadings at the various ¢,-angles 
are tabulated next. In each case the influence coefficients are 
given for @ = dx, bd: + 5, d2 + 10 deg due to edge loadings at dy. 
The two edge loadings considered for the crown opening are shown 
in Fig. 3(b). As was the case for @ = 90 deg, the first lines in the 
tabulation correspond to rH -loadings, the second lines to Mg- 
loadings.” 


2 An illustration of the accuracy of the numerical results is afforded 
by computing the edge influence coefficients for a shell with r/h = 
200, ¢2 = 10 deg. For this case, the cot ¢ = 1/¢ approximate method 
(Galletly, [3]) should be quite good. The results obtained were 


For rHg: = 1, Eré = 237.04 and Er?V = 20,735 
For Mg: = 1, Eré = 20,735 and Er?V = 4,381,600 


The discrepancy between these values and those in Table 3 (for r/h 
= 200, ¢: = 10 deg) is less than 2 per cent. 


J 
LINE 1: * 1 


(a) = 90° 


LINE 3: rg, LINE 4° Mg, 


(b) 


Fig. 3 Axisymmetric edge loadings considered in influence coefficient 
tabulations 
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Influence Coefficients Due to Uniform Internal Pressure 


From equation (34) in the Appendix it can be seen that the dif- 
ferential equations of equilibrium for hemispherical shells sub- 
jected to edge bending loads and uniform pressure are homo- 
geneous, i.e., there is no term on the right-hand side involving 
the pressure. This means that the membrane equations give the 
stress resultants and deformations due to uniform pressure. 
Thus, due to uniform internal pressure, the influence coefficients 
are: 


V=0 
rNg = pr*/2 
rNe = pr?/2 (7) 
Eré = pr?/2(r/h) (1 — v) sing 
My, = 0 
Ms = 0 
where 
Ng = meridional direct stress resultant (+ when tensile) 
Ms = circumferential bending stress resultant (+ when pro- 
ducing tension on inner surface) 
p = uniform internal pressure 


and the other symbols are as defined previously. 


Effect of Vertical Line Loadings 


The influence coefficients given in Tables 1 to 3 are for uni- 
formly distributed bending moments and horizontal shearing 
forces applied at the edges. The coefficients due to uniform pres- 
sure have just been discussed in the preceding section. Another 
loading case of interest is depicted in Fig. 4(a), i.e., uniformly dis- 
tributed concentrated vertical loads applied around the edges. 
In this section, one method of solving this problem is given. 

The vertical load V at @ = a may be regarded as being the re- 
sultant of two other loadings. One of these, V cot a, is assumed 
to act horizontally. The other loading, V/sin a, is assumed to 
act along the tangent at @ = a. It is easily verified that these 
two loadings have a resultant equal to V acting vertically. The 
problem shown in Fig. 4 can then be solved if means are available 
for solving the two partial problems shown in Figs. 4(b and c). 

The problem shown in Fig. 4(b) can easily be solved with the aid 
of Tables 1 to 3. That in Fig. 4(c) has the following solution: 
V sin @ 8) 
sin? (8: 


Ny = = 


Ehé (9) 


V (10) 


Ms = Ms =0 


The fact that V and dV /d®@ are zero, and hence Mg and Mg, can 
be checked by inserting the relevant quantities in the following 
compatibility equation relating the middle surface strains, €, and 
€, and the rotation V: 


Fig. 4 Resolution of axisymmetric vertical edge loadings into auxiliary 
loadings 
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phere-cy 
open-crown hemisphere. 


V = — &) cot — r2/ri d/dg(e) (11) 
where 
= 1/Eh[Ng — vNo] 
€ = — 


re = mr, = r for a hemisphere 


(12) 


lilustrative Example: Hemisphere-Cylinder Intersection 


As an illustration of the use of the influence coefficients the 
problem of the intersection of a cylinder with a hemisphere will 
be considered (see Fig. 5). It will be assumed that the opening in 
the hemisphere subtends an angle of 2a at the center of the 
hemisphere and the unknown stress resultants at @ = @ will be 
denoted by: 


Qa = transverse shearing force/unit length 
M, = meridional bending moment/unit length 
Na = meridional direct force/unit length 


The following relation between Q, and Ng also holds: 


Ne = —Q,ccta + (13) 


To determine the edge forces acting on the ends of the cylinder 
Na and Qq are resolved horizontally and vertically. Thus, 


Horizontally outward: Qa sina — Ng cosa (14) 


= sina + cosa cota — cos [from equation (13)} 


Qa Pp 
— — cosa 
sin @ 2 
R 
-H.- cot (15) 
where Hg is the horizontal shearing force. 
Vertically upward: On cos a + Ng sina (16) 


w 


intersection. Edge forces and moments on 
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| 
V sina V sing Vsina V sing 
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= Qa cosa — Qa cosa + 5 sin a [from equation (13)] 


pk 
2 


(17) 


Thus the vertical forces on the cylinder and hemisphere balance 
each other. A horizontal force of magnitude Hg — pR/2 cot a 
has, however, to be applied inward at the cylinder ends to satisfy 
horizontal force equilibrium requirements. To satisfy moment 
equilibrium requirements a moment Mg must also be applied to 
the cylinder ends. The final force system (Hq and M, still un- 
known, however) acting on the cylinder and hemisphere is as 
shown in Fig. 6. 

The magnitude of M, and Hq are determined by demanding 
that the final displacements and rotations (due to internal pres- 
sure and edge bending loads) of both cylinder and open-crown 
hemisphere be alike. 

The edge influence coefficients for a cylinder are (Roark, 1943): 
Due to uniform internal pressure p: 


2 
= (2—yv); =0 (18) 


Due to uniformly distributed unit shears perpendicular to the 
axis of symmetry at the end of the cylinder: 


1 


2DB6”’ 
Due to uniformly distributed unit moments at the end of the 
cylinder: 


1 1 


= = BD (20) 


mean radius of cylinder | 


thickness of cylinder 


12(1 — v?) 


and the superscripts p, H, and M denote internal pressure, shears, 
and moments, respectively. 

Using these relations the total displacement and rotation of the 
cylinder at its junction with the hemisphere is given by: 


pk? pR ) 1 M. 
6, 2Eh, (2 — v) (1. cota opp: + 


PR 


BD 
In these equations 6 was taken positive outward and V positive 
when the tangent rotated counterclockwise on the left. 

The displacement and rotation of the open-crown hemisphere 
due to the Hg and M, can be obtained from Tables1to3. Thus, 
relations are obtained of the form (assuming no interaction be- 
tween the edges of the hemisphere at ¢ = a and @ = 90 deg): 


Er°V, = ( )rHa + ( )Ma (24) 
Erb, = ( + ( )Ma (25) 


where the parentheses indicate numerical values of the edge in- 
fluence coefficients and the subscript h denotes the hemisphere at 


(22) 


(23) 


Journal of Engineering for Power 


Shears, forces, and moments acting at cylind Pp 


Fig. 6 
section 


@ = a. The displacement and rotation of the hemisphere due 
to uniform internal pressure are obtained from equation (7). 
Thus the total diysfcement and rotation of the hemisphere at 
@ = @ are given by 


Er°V, = ( )rHa + ( )Ma (26) 


pre f(r 
Eré, = 2 (<) (1 — »p)sina +( + ( )Ma (27) 


As the cylinder and hemisphere must eventually fit together at 
their common junction, the following compatibility equations 


must be satisfied: 
Vi 


From equations (22), (23), (26), (27), and (28) two simultane- 
ous equations for Hg and Mg are obtained. Having solved these 
equations it is then a simple matter to obtain the stress distribu- 
tion in the vicinity of the edge of the hemisphere by using the in- 
fluence coefficients given in Tables 1 to3. The stress distribution 
throughout the cylinder can also easily be found using the tables 
given in Timoshenko. 

Discussion of the Differential Equations Used in the Analysis. The 
differential equations used to calculate the influence coefficients, 
i.e., equation (35) is based on classical small-deflection theory. 
The assumptions made in that theory and its limitations should 
be remembered when using the influence coefficients. In particu- 
lar, the following points should be checked at the end of an 
analysis: 


(a) That the shell deflections are small in comparison to the 
shell thickness. 

(b) That no portions of the shell are inelastic. 

(c) That the direct forces Ng and Ny are small in comparison 
with their critical compressive values. 


For thick shells, i.e., those in which the shell thickness is not 
small in comparison to the principal radii of curvature, it is 
probable that shear deformations are important. These latter 


(28) 
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have not been considered herein. However, they will be investi- 
gated in the near future. 

With regard to item (b), if the final stress distribution in the 
shell should be substantially inelastic, then the limit load of 
the shell should be investigated. In this connection, the work 
of Drucker and Shield [2] should be of value. 

Creep effects have also not been considered in the present work, 
although there will be many instances in which they will be im- 
portant. 


APPENDIX 


Hemispherical Shell Subjected to Edge Bending Loads and 
Uniform Internal Pressure 

Differential Equations of Static Equilibrium. The differential equa- 
tions of equilibrium for a spherical segment subjected to edge 
bending loads can be obtained by specializing the equations given 
in Timoshenko [8]. 1t will be found that if a uniform pressure is 
also taken into account the differential equations remain un- 
changed. Thus one obtains the following equations: 


L(U) + »/r U = EhV 
L(V) — v/r V = —U/D 


(29) 
(30) 


rQ¢ 
transverse shearing force/unit length 
angle of rotation of a tangent to a meridian 
flexural rigidity of the shell = Eh*/12(1 — v?) 
Young’s modulus 
thickness of the hemisphere a 

= radius of the hemisphere 

= Poisson’s ratio 


(31) 


and 
L(...) = 1/r{(.. .)” + cot Gf... .)’ — cot? Gf. . .)] 


and the primes denote differentiation with respect to the colati- 
tude angle @. 
Introducing now the dependent variable Z defined by 


Z=U+KkV (32) 
where 
k => 


wer = [1211 — v®)(r/h)? — 


—D/r(v + ip?*r) 
(33) 


the two equations (29) and (30) can be reduced to the following 
single complex, second-order, ordinary linear differential equa- 
tion: 
L(Z) iwZ = 0 (34) 

where 

i=vV-1 
The minus sign in equation (34) is obtained if the positive sign 
is taken for k in equation (33), and vice versa. 

If the two linearly independent solutions of 


L(Z) — iwZ = 0 (35) 


ZA = 11 + ip 
(36) 
= X2 + tyr 


then it can be shown (Burrows [1], Galletly [4] ) that the solutions 
for U and V (where Z = U + kV) may be written as: 
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U = Ax + By + Cre + Dyz (37) 


EhrV = A(vx, — p*ry:) + + 
+ C(vr2 — + D(vye + w*rae) (38) 
Knowing U and V the stress resultants can be found from the 
relations: 
—U cot o/r + pr/2 
—U'/r + pr/2 
—D/r(V’ + vV cot >) 
—D/r(V cot + 
U/r 
Q¢/sin 


meridional direct stress resultant (+ when tensile)* 
circumferential direct stress resultant (+ when tensile) 
meridional bending stress resultant (+ when producing 
tension on inner surface) 
circumferential bending stress resultant (+- when pro- 
ducing tension on inner surface } 
transverse shearing force (+ when directed toward the 
center of the hemisphere for the larger value of $) 
horizontal shearing force (+ in the same sense as Qg). 
With no load acting on the shell, the relation Ng = 
—Qz cot @ holds. Hg is the horizontal resultant of 
Ngand Q¢. The vertical resultant is, of course, zero. 
The horizontal displacement of a point on the shell 6 (+ out- 
ward) can be found from the relation: 


Ehéb = rsin — 


Solution of the Differential Equations of Equilibrium. [quation (35) 
is a linear, complex, homogeneous, ordinary differential equation 
of the second order. . As the equation is complex the complete 
solution will contain [as shown by equation (37)] four linearly in- 
dependent solutions. These solutions were obtained by substitut- 
ing Z = x + ty in equation (35) and equating real and imaginary 
parts. One obtains: 


(40) 


= cot? — pry | 
= + ycot? d + wre | 


(41) 


By substituting 


(42) 


= Us 
the two second-order equations (42) can be reduced to the follow- 
ing four simultaneous first-order differential equations: 
= —u, cot d + uw cot? d — wrug 
= —u; cot d + us cot? @ + wru 


= U3 


(43) 


The four equations (43) were integrated on a high-speed digital 
computer using the Runge-Kutta integration method. The step 
length used was !/. degree. The integration process was com- 
menced at @ = 90 deg using the following initial values: 


Ist solution 
2nd solution 


3The N’s, M’s, H’s, etc., are per unit length. 
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As the shell has two edges the problem is really a two-point 
boundary value problem. However, as the differential equations 
are linear, it is only necessary to assign initial values at one point 
(Milne, [5]). 

Once the numerical solutions have been obtained, the influence 
coefficients are determined by using equations (37) to (40). 
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ERRATA 


Journal of Engineering for Power, Trans. ASME, series A, vol. 81, July, 1959 


Three errors have been brought to our attention by Escher Wyss Ltd. They 
occur in the special issue on gas turbines. They are as follows: 

Page 321, Table 9 Gas turbines in naval propulsion applications. The manu- 
facturer of JN 10, 10,000 hp should be Mitsui (not Escher Wyss). 

Page 322, paragraph entitled 10,000-shp Prototyge Closed-Cycle Plant. The first 
sentence should read: ‘A closed-cycle propu yn plant, utilizing air as the 
working fluid, is under construction by Mitsui’ (not Escher Wyss). 

Page 326, Table 13 Characteristics of marine gas turbines installed or building. 
The manufactiet 'y “"W JN 10, 10,000 hp should be Mitsui (not Mitsubishi). 
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